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ABSTRACT 


The  120-foot  antenna  of  the  Haystack  Microwave  Facility  and  the  60-foot  antenna 
of  the  Westford  Communications  Terminal,  both  operated  by  M.l.T.  Lincoln 
Laboratory,  were  coupled  to  form  a  planetary  radar  interferometer  operating  at 
X-band  and  were  used  to  observe  Venus  at  a  wavelength  of  3.8  cm  during  the 
1967  inferior  conjunction.  The  antennas  are  separated  by  approximately  4000  feet 
along  a  line  22°  east  of  north  At  maximum  projection  in  the  direction  of  the 
planet,  this  baseline  gives  a  fringe  spacing  of  5  seconds  of  arc,  or  a  maximum 
of  about  10  fringes  across  the  planetary  disk  at  inferior  conjunction. 

By  transmitting  a  CW  signal  from  the  120-foot  antenna  and  frequency  analyzing 
the  received  echo,  it  was  possible  to  resolve  the  planetary  surface  scattering 
into  strips  parallel  to  the  apparent  axis  of  rotation.  Crosscorrelation  of  the 
complex  frequency  components  obtained  at  the  two  sites  yielded  corresponding 
spatial  Fourier  components  which  resolved  the  scattering  along  the  strips.  With 
1-Hz  frequency  resolution  and  a  maximum  of  10  fringes  along  the  rotation  axis, 
the  planetary  hemisphere  visible  to  the  radar  during  inferior  conjunction  was 
mapped  with  approximately  100  resolution  intervals  along  a  direction  perpendic¬ 
ular  to  the  apparent  rotation  axis,  with  10  resolution  intervals  in  the  orthogonal 
direction. 

For  a  limited  region  on  the  planet,  surrounding  the  center  of  the  visible  disk, 
higher  resolution  was  obtained  by  transmitting  pulses  of  500-psec  effective 
length  The  pulse  resolution  enabled  the  planet  to  be  resolved  in  echo  delay, 
leaving  only  a  twofold  hemispheric  ambiguity  to  be  resolved  by  the  interferom¬ 
eter.  In  addition,  in  the  range-gated  observations  the  effects  of  significant  in¬ 
terferometer  sidelobes  (arising  from  the  limited  range  of  projected  baselines 
available)  were  avoided. 

Maps  obtained  from  the  observations  show  Venus  to  be  smoother  on  the  average 
than  the  moon  at  3.8  cm,  although  some  regions  of  the  planet  exhibit  strong  local 
radar-scattering  enhancement.  The  positions  of  these  regions  agree  well  with 
those  previously  reported  if  the  rotation  period  of  Venus  is  assumed  to  be  earth- 
synchronous  at  243. 16  days  retrograde. 


Accepted  for  the  Air  Force 

Franklin  C.  Hudson 

Chief,  Lincoln  Laboratory  Office 
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A  RADAR  INTERFEROMETER  STUDY  OF  VENUS  AT  3.8  cm 


I.  INTRODUCTION 

The  planet  Venus  is  surrounded  by  an  optically  opaque  atmosphere  such  that  neither  the 
nature  of  its  surface  nor  its  period  of  rotation  can  be  determined  with  eertainty  from  telescopic 
observation.  However,  its  atmosphere  is  reasonably  transparent  at  radio  wavelengths,  and  much 
information  has  been  gained  from  earth-based  radar  observations.  For  example,  radar  obser¬ 
vations  made  during  both  the  1962  and  1964  inferior  eonjunetions  at  the  Jet  Propulsion  Laboratory 
(JPL)*  have  shown  anomalous  peaks  in  the  speetrum  of  the  reeeived  echo  which  appear  to  be 
caused  by  regions  of  loeally  enhaneed  radar  baekseatter  on  the  planetary  surfaee.  The  locations 
of  these  regions,  as  well  as  the  surface  rotation,  were  determined  through  a  least-squares  fitting 
procedure.  From  these  observations,  a  sidereal  rotation  period  of  250  (+4,  —7)  days  retrograde 
and  a  north  polar  direction  of  255°  (  +  10°,  —4°)  in  right  ascension  and  68°  (±4°)  in  declination 
were  obtained. 

Measurements  in  which  echo  delay,  as  well  as  frequency  information,  was  used  were  re- 
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ported  from  the  Arecibo  Ionospheric  Observatory  (AlO)  for  the  1964  inferior  conjunction  of 
Venus.  These  data,  which  are  intrinsically  more  powerful  as  a  means  of  determining  pole  posi¬ 
tion,  yielded  a  rotational  period  of  245.1  (±2)  days  retrograde  at  270.3°  (±1°)  right  ascension 
and  66.7°  (±1°)  declination. 

Measurements  involving  both  a  simple  CW  echo  analysis,  as  well  as  a  more  complex  com¬ 
bination  of  delay  and  frequency  data,  can  benefit  from  the  inclusion  of  additional  information 
concerning  the  spatial  distribution  of  the  received  echo.  Since  the  disk  diameter  of  Venus,  even 
at  inferior  conjunction,  does  not  exeeed  roughly  one  minute  of  arc,  this  information  to  be  useful 
must  have  a  resolution  of  better  than  a  few  tens  of  seconds  of  arc.  With  this  in  mind,  therefore, 
Lincoln  Laboratory  in  late  summer  of  1967  undertook  an  observing  program  in  which  the  existing 
120-foot-diameter  Haystack  planetary  radar  system  and  the  Westford  60-foot-diameter  antenna 
were  interconnected  at  a  frequency  of  7840MHz  as  a  phase-coherent  radar  interferometer. 
Measurements  using  both  simple  CW  transmissions,  as  well  as  more  complex  time-eoded  signals, 
were  made. 

In  this  report,  the  theory  of  radar  interferometry  is  first  developed,  followed  by  sections 
dealing  with  surface-mapping  techniques  and  descriptions  of  the  actual  equipment  and  data- 
reduction  procedures  employed.  Only  a  modest  amount  of  interpretation  of  the  results  has  been 
attempted,  sinee  it  is  anticipated  that  several  journal  articles  will  treat  this  aspeet  of  the  re¬ 
search  in  greater  detail. 
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II.  THEORY  OF  RADAR  INTERFEROMETRY 


A.  Interferometer  Geometry 


The  interferometer  geometry  is  shown  in  Fig.  1.  The  coordinate  system  used  is  centered 
at  Haystaek.  r^  is  a  fixed  veetor  from  the  intersection  of  the  azimuth  and  elevation  axes  of 
Haystack,  to  the  interseetion  of  the  azimuth  axis  and  a  horizontal  plane  through  the  elevation 
axis  of  Westford;  r?  is  a  vector  in  that  plane  representing  the  offset  of  the  elevation  axis  of  the 
Westford  antenna;  and  r^  is  a  unit  vector4  toward  the  apparent  position  of  the  subradar  point  on 
the  planet  being  observed  (parallax  between  the  sites  can  be  neglected).  These  vectors  can  be 
expressed  as 


r.  =  D  sinE„  +  D  eos  E„,  sinA„,  i~A  +  D  eos  E„r  cosA„.l,^ 

1  W  ZE  W  W  EA  W  W  NO 

(differences  in  zenith  neglected) 


/\  A 

r»  -  d  sin  A  iT„  A  +  d  cos  A  E7~ 
2  p  EA  p  NO 


A  A  A  A 

r0  =  sinE  i~T^  +  eos  E  sin  A  i^  A  +  eos  E  eos  A  L7~ 
3  p  ZE  p  p  EA  p  p  NO 


(1) 

(2) 

(3) 


where  D  (  j  r^  |  )  is  the  distance  between  Haystack  and  Westford,  d  (  |  r^|  )  is  the  elevation  axis 
offset  of  Westford,  A^  and  E^  are  the  apparent  azimuth  and  elevation  of  the  planet,  and  A^  and 
E^r  are  the  azimuth  and  elevation  of  Westford  as  seen  from  Haystack.  The  portion  of  the  differ¬ 
ential  delay  which  is  azimuth  and  elevation  dependent  is  given  by 

a  i) 

r  ,  =  —  (r,  +  r„)  •  r^  =  —  [sinE„,  sinE  +  eosE...  eos  E 
d  e  1  2  3  e  1  W  p  W  p 


x  cos  (A  —  Ad]  +  — 
'  p  W  1  c 


cos  E 


(4) 


(  ZENITH) 


WEAST> 


ELEVATION  AXIS 
OF  WESTFORD 


Fig.  1  .  Interferometer  geometry. 


t  Throughout  this  report,  the  symbol  (A)  is  used  for  o  unit  vector. 
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where  cg  is  the  velocity  of  propagation  along  the  baseline.  Results  of  a  baseline  survey  gave 
the  following  values  for  the  baseline  parameters: 

D  =  4066.16  feet 
Aw  =  S21  °53'47"  W  ±  6" 

E,  =  -  1  °22'53"  ±  4" 

w 

d  -  12.5  inches. 


is  a  funetion  of  time  as  a  result  of  the  motion  of  the  planet  and  gives  rise  to  a  differential 
Doppler  effect  or  fringe  rate  of  (co/2n)  [dr^(t)/dtl.  The  change  of  with  position  in  the  sky 
produces  an  interferometer  fringe  pattern  for  which  it  is  more  convenient  to  use  celestial  co¬ 
ordinates  wherein 


r  ,  «  —  [sin<5_,  sin  6  +  cos  6^.  eos6  eos(L  -  L„)l  (5) 

d  c  1  B  p  B  p  p  B  1 

where  and  6^  are  the  hour  angle  and  deelination  of  the  baseline.  The  expression  is  only  ap¬ 
proximate,  as  refraetion  ean  only  be  taken  into  aeeount  to  first  order  (plane  parallel  atmosphere) 
in  this  coordinate  system.  However,  the  expression  is  quite  preeise  enough  to  eompute  the 
fringe  spacing.  Expressing  r^  in  eelestial  coordinates  ( |  r ^  I  is  mueh  smaller  than  |  r^  I  and  is 
neglected),  the  projections  of  r ^  in  a  plane  normal  to  the  direction  of  the  planet  are 

drd  — 

e  -gg-  =  (rd)N  =  D  [cos<5p  sindR  -  cosSR  sinSp  cos(L  -  LR)|  (6) 


cos  6 


drd  — 

dT7  =  (lVw  =  D  eos6B  sin(Ls-LB) 


(7) 


The  fringe  spaeing  in  radians  is  the  reciprocal  of  the  projected  baseline  component  in 
wavelengths.  Figure  2  shows  a  plot  of  the  projected  baseline  in  wavelengths  for  Venus  on  the 
day  of  the  1967  inferior  conjunction. 


Fig.  2.  Projected  baseline  (i)  for  Venus. 
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B.  Earth-Venus  Geometry 


Since,  in  this  discussion,  we  are  primarily  interested  in  the  mapping  of  the  planetary  sur¬ 
face,  the  orbital  and  rotation  parameters  will  be  assumed  known.  Any  inconsistencies  which 
might  be  attributed  to  errors  in  these  parameters  will  be  discussed  in  See.  VI -C.  The  orbit  of 
Venus  will  be  described  by  its  geocentric  right  ascension  RA^,  declination  6  ,  and  distance. 


HAYSTACK 


Fig.  3.  Earth-Venus  geometry. 

It  will  also  be  assumed  that  Venus  is  a  sphere  of  radius  r  which  rotates  with  an  angular  velocity 
fip  relative  to  the  local  star  system.  The  direction  of  is  given  by  its  right  ascension  RA^ 
and  declination  6^.  Figure  3  shows  the  Earth-Venus  geometry.  Three  Venus-centered  co¬ 
ordinate  systems  are  used  in  the  complete  description  of  the  radar-mapping  process.  The  first 

A  A 

is  a  system  (X,  Y,  Z)  where  X  is  directed  toward  the  transmitter,  and  Z  is  in  the  direction  of 

—  a 

the  projection  of  the  apparent  rotation  axis  ona  plane  normal  to  X.  In  this  system,  the  sub- 

A 

radar  point  is  (r^,  0,  0).  The  second  coordinate  system  (X,  E,N)  is  the  first  rotated  so  that  N 
is  in  the  direction  of  increasing  declination.  The  third  coordinate  system  (X  ,  Y  ,  Z  )  is  fixed 
to  the  planet  so  that  Z^  is  the  pole  or  the  direction  of  ft  ,  and  X^  goes  through  the  circle  of  zero 
longitude,  —40°  being  defined  as  the  longitude  of  the  subradar  point  on  20  June  1964  at  0*1  UT 
(Ref.  1).  Thus, 

Xp  -  r  cos  (Eat)  cos  (Long) 

Yp  =  r  cos  (Eat)  sin(Eong) 


Z  =  r  sin  (Eat) 
P  P 


(8) 


where  Eat  and  Eong  are  the  latitude  and  longitude  of  a  point  on  the  surface  of  Venus.  Conversion 

A 

from  (X,  E,N)  to  (Xp,  3  ,  Zp)  requires  four  rotations.  The  first  rotation  about  E  by  the  angle 
6p  makes  NT  parallel  to  the  earth's  axis 

X'  =  X  cos  6  +  N  sin  6 

P  P 

E'  =  E 


N1  =  N  cos  6  —  X  sin  6 

P  P 


(9) 
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The  second  rotation  about  NT  by  the  angle  A  makes  EM  in  the  opposite  direction  to  the  right 
ascension  of  the  rotation  axis,  so  that 


A  =  RA  -  RA  +  90 c 
P  a 


(10) 


as  illustrated  in  Fig.  4.  Thus, 

X"  =  X'  cosA  —  E’  sinA 
EM  =  E'  cosA  +  X'  sinA 
NM  =  N' 


(11) 


Fig.  4.  View  of  Eorth-Venus  geometry  from  celestial  pole. 


The  third  rotation  about  XM  by  O  makes  N,M  parallel  to  Z^,  where 
O  90°  -  6 


(12) 


so  that 


Xm  =  XM 

E,M  =  EM  cos  G  4  NM  sin  O 


Z  -  N' 1  ’  =  NM  cos  0  -  E"  sin  G 
P 


(13) 


From  Eqs.(9),  (11),  and  (13), 
X’ 


X  cosA  cos  6  —  E  sinA  +  N  sin 6  cosA 

P  P 


E,M  =  X(sinA  cos  6  cos  G  —  sin  G  sin 6  )  4  E(cosA  cos  G) 

P  P 

+  N(sinG  cos  6  +  cos  G  sinA  sin  6  ) 

P  P 

Z  -  X(— cos  G  sin  6  —  sin  G  sinA  cos  6  )  —  E(cos  A  sin  G) 

p  p  P 


+  N(cosG  cos  6  —  sin  G  sin 6  sinA) 
P  P 


The  final  rotation  about  Z^  by  angle  P  makes  X^  go  through  the  definition  of  zero  longitude. 


(14) 

1 


p  =  p  4  n  t 

o  p 


(15) 


5 


Thus, 


where  P  is  the  angle  between  X  and  the  direction  (HA  +  90°)  when  t  =  0. 
o  p  a 

X  =  X"'  cos  P  -  IT"  sin  P 
P 


Y  =  KMI  cos  P  +  X,M  sinP 
P 


(16) 


Since  the  subradar  point  is  given  by 

X  r 
P  P 

E  =  0 

N  -  0 

the  latitude  of  the  subradar  point  is 


_  i 

LatT^  =  sin  (—  cos  0  sin  6  —  sin  0  sin  A  cos  6  ) 

H  p  p 


(17) 


from  Eq.(14)  by  setting  E  N  0.  Similarly,  the  longitude  of  the  subradar  point  is 


Long^  P  +  tan 


■( 


sin  A  cos  6  cos  0  —  sin  0  sin  6  ■ 
_ P _ P 

cos  A  cos  6 


(18) 


from  Hqs.(14)  and  (16). 

The  transformation  from  (X,  Y,  Z),  in  which  Z  is  the  apparent  rotation  axis,  to  (X,  E,N) 
requires  one  rotation  by  the  angle  D  which  is  computed  from  the  components  of  center-to-limb 
Doppler.  The  apparent  rotation  consists  of  the  projection  of  the  planet's  rotation  on  the  plane 

/N 

normal  to  the  direction  X  plus  the  rotation  due  to  the  change  in  position  of  the  planet, 


H  =  12  (sin  6  cos  6  —  sin  6  cos  6  sinA) 

N  p  a  p  pa 


dRA  r  .12  cos  ( LatTI)  cos  L 
pee  H  p 

dt  X„ 


(19) 


and 


d6 

12...  —  12  cos  6  cos  A 
W  p  a  dt 


(20) 


where  LatTT  is  the  latitude  of  Haystack,  r  is  the  radius  of  the  earth,  12  is  the  rotation  of  the 

II  J  e  e 

earth,  and  X  is  the  distance  to  Venus.  The  last  term  is  a  correction  for  the  earth's  rotation 
o 

which  has  to  be  included  if  the  right  ascension  and  declination  of  the  planet  are  computed  for  the 
geoccnter.  From  Eqs.(19)  and  (20), 


D  -  tan  1  (fiw/«N)  • 


(21) 


The  values  of  X  ,  HA  ,  and  6  can  be  obtained  from  ephemeris  tables.  The  following  rota- 
2  °  P  P 

tion  constants  were  assumed  for  Venus: 

6  =  66.7° 

a 

RA  =  270.3° 
a 

12^  =  1  rotation  per  243.16  days  —  retrograde. 


Values  of  the  latitude  and  longitude  of  the  subradar  point  were  checked  with  computations  by 
3 

1. 1.  Shapiro. 
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C.  Relation  of  Received  Signal  to  Planet's  Surface 

Using  the  coordinate  system  (X,  Y,  Z)  defined  in  Sec.  B  above,  the  wave  incident  on  the  sur¬ 
face  of  the  planet  is  related  to  the  transmitted  signal  xfj  (t)  by 


x.(w,  V,  Z)  -  oxT(tc')  exp  {—  [Xo(t)  -  X(t)|} 


where 


x,rM  = 


j  x,j, ( t )  c”1^  dt 


(22) 


(23) 


X^(t)  is  the  distance  to  the  center  of  Venus,  and  o  is  an  attenuation  constant.  The  plane  wave 
assumed  in  Eq.(22)  will  be  justified  in  Eq.(40b).  If  it  is  now  assumed  that  each  element  of  unit 
surface  area  on  the  planet  inclined  at  an  angle  G  to  the  incident  beam  returns- a  signal 


x^(w,  Y,  Z)  -  x.(w,  Y,  Z)  s(Y,  Z) 


then  the  integrated  return  will  be 

Xj^(w)  =  a  ^  xT(w)  s(Y,Z)  exp{-  ~  [XQ(t)  -  X(t)|}  ^ 


dYdZ 
sO 


(24) 


(25a) 


where 


Q  x  (i-yZ-z2)1/2 

cos  9  -  —  -  - - - - 


(25b) 


Expanding  XQ(t)  and  X(t)  into  a  constant  and  time-changing  component  due  to  rotation  and  radial 
velocity  V  , 


X  (t)  -  X(t)  =  X  (t  )  -  X(t  )  +  V  t  +  YQt 
o  o  o  o  r 


(26) 


so  that  Eq.  (  25)  becomes 


xr(cj)  =  a 


ss 


/  i\  7)  -iwT(Y,  Z)  dYdZ 
xT<aj  )  S<Y’Z>  e  ^FO 


(27) 


where  t  is  the  delay  to  the  element  at  (X,  Y,  Z) 

2  [X  (t  )  -  X(t  )| 

1  o  o  o’ 


(28) 


and  u;?  is  the  Doppler-shifted  frequency 


w  -  —  (2Yfl  +  2  V  ) 
c  r 


If  the  planet's  orbit  is  known,  it  is  possible  to  refer  the  signal  to  the  center  of  the  planet  by 
continuously  shifting  one  of  the  local  oscillators  by  the  amount 


^  2V  (t) 
c  r 

so  that 

x(w)  =  xR(w)  exp[i£  2Xo(t)l  =  a  j  (  xT(w”)  s(Y,Z)  e"la;T(Y'Z)  (29) 

where  con  =  w  —  (2coYfi/c).  The  simplest  signal  to  transmit  is  a  sine  wave  (to  within  the  fre¬ 
quency  stability  of  the  standard),  in  which  case 
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(30) 


I  XT(w)|  =  1 ±o(0)  -  L>o) 

For  a  single  antenna  whose  antenna  beam  is  much  larger  than  the  planet,  the  receiver  spectral 
power  is 


|x(")|2=°2M 


xT(w")  x*(w"t)  s(Y,  Z)  s*(Yt,  zt) 


v  -iwr(Y.Z)  +iwT(Yt,  Zt)  dYdZdYtdZt 
x  e  e  - 2 - 


COS  0 


n  rc(u)  —  lo  )  i 


dZ 


(31) 


cos  0 

where  the  power-scattering  function  S(Y,  Z)  is  a  real  function 
S,Y,  Z,  - 

The  function  within  the  integral  is  a  delta  function  in  y  and  7  so  that  the  power  at  frequency  w 
originates  from  a  line 


v v  ry  v  -io)T(Y,  Z)  iwT(Y-y,  Z-z)  ,  , 

s(Y,  Z)  s*(Y  —  y,  Z  —  z)  e  e  J  9  '  dydz 


(32) 


Y  = 


(wQ  -  w)  c 
2  ft  Cl) 


(33) 


parallel  to  the  apparent  rotation  axis,  a  is  the  two-way  attenuation  which  is  related  to  the  antenna 
gain  G  by 


2  GZ\2 
a  - 


64tt3X4 

o 


provided  S  is  defined  as  unity  for  a  perfect  isotropic  scatterer.  For  simplicity,  it  is  convenient 
to  change  the  scale  of  the  coordinates  so  that 

X  =  x/r 

n  p 


Y  =  Y/r 
n  /  p 

Z  =  Z/r 
n  p 

and  to  introduce  the  center-to-limb  Doppler  shift  F  so  that 

i  /  \  1 2  2  2  f  0,  FAti)  ^  ,  dZn 
1  '  M  p  J  '  27r  '  n  cos  0 


(34) 


(35) 


where 


Agj  =  cj  —  c o 

o 


If  the  spectral  analysis  of  the  received  signal  is  performed  with  a  frequency  resolution  function 
f ( Cl) ) ,  which  for  Fourier  analysis  of  a  sample  of  length  T  is 


sin  [(w  -  wc)  -|] 


then  the  output  of  each  "frequency  channel  cM  is 


(36) 
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I  xc  I  2  =  J  lx(“)|2  *<.(">  ^ 


do; 

2n 


(37) 


since 


=  C  x(t)  e  ia;t  dt 
Jr* 


The  signal  y(t)  received  by  the  second  antenna  from  an  element  on  the  planet  is  shifted  in  phase 
because  of  the  different  distances  of  the  element  from  the  two  stations.  From  Fig.  5, 


y{w)  =  0a  x^w")  S(Y.Z)  e'l<JT(Y’  Z) 

2niz  try  2niY  Lz  Ary 
^  n  Z  n  Y  dYdZ 

x  e  e  - 

cos  0 


-10JT 


(38) 


where  is  the  delay  differential  to  the  subradar  point,  given  by  Eq.  (4). 


Fig.  5.  Relative  distance  between  twa  antennas  and  paint  an  planet. 


and  l  are  baseline  projections  in  the  Y  and  Z  directions  in  units  of  wavelengths  per 
planetary  radius.  Thus, 

,  (r  )Y 

Y  2?rcX 

o 

(r  )ry  Wr 


l  _  z,  p 
Z  2ttcX 


(39) 


where  the  Y  and  Z  components  of  the  baseline  can  be  obtained  from  the  N  and  W  components 
given  in  Eqs.  (6)  and  (7). 


(r  )z  =  (r  )N  cos  D  +  (r  )w  sinD 


(r  )y  = -(r  )w  cosD  +  (r  )N  sinD  .  (40a) 

The  phase  differences  expressed  in  Eq.  (34)  assume  the  received  signal  components  to  be  plane 
waves,  but  are  in  error  at  most  by  only 
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D*.  / 

2^  ^Z7i  ~  10  radian  (o?/27r  =  7840  MHz,  and  X  -  0.3  astronomical  unit)  (40b) 
o 


as  a  result  of  the  wavefront  curvature. 

From  Fqs.(30)  and  (38),  the  cross-spectral  power  of  received  signals  is 


y(o;)  x*(w)  =  P 


-2  M 


xT(w')  x*  (W't)  S(Y,  Z)  s*(Yt,  Zt) 


-io)r(Y,Z)  iwT(Yt,  Zt)  'la>Td  2lrlZ nZ  dYdZdYtdZt 

a.  C  C  0  0  0  - - — — ^ - — — - 

cos1'  0 


2„  2 
“  rpSo 


fs,- 


-1WT 


AwF 


Z  )  e 
n 


,  2 7riZ  i_  - 

d  n  Z 

e  e 


27riAo;Fi^  dZ 


n 


cos  G 


(41) 


for  a  CW  transmission.  If  y(t)  is  rotated  by  exp[io)r^),  then  the  lines  of  constant  phase  become 
stationary  on  the  planet.  Further  rotation  by  exp  [27riAu;Fiy]  makes  the  line  of  zero  phase  per¬ 
pendicular  to  the  Doppler  axis.  If  the  frequency  resolution  <5co  used  is  sufficiently  small  that 

<5a;F27riY  «  1  (42) 


then  the  power  and  cross-power  outputs  for  each  channel  are 


2  2 
r  a 
P 


I 


dZ 

S  (Z  )  - 

c  n  cos  0 


(43) 


and 


2  2„ 

y*  =  rp«  P 


C1'  C 


S  (Z  )  e 
c  n 


2 TriZ  i_ 
n  Z 


dZ 
_ n_ 

cos  0 


(44) 


where  is  the  scattering  function  integrated  over  the  Doppler  strip. 

If  the  transmitted  signal  is  phase-reversal  coded  by  multiplication  with  a  code  c(t)  which 
has  an  element  length  r ^  and  repeats  every  seconds,  then 


xT(w)  =  KiQ( Cl?  -  WQ)  ©  c*(w) 


(45) 


The  decoding  of  the  received  signal  is  achieved  by  multiplication  with  the  transmitted  code  de¬ 
layed  by  different  amounts.  There  are  Tr/Te  delay  channels  formed  which  can  then  be  spectrum 
analyzed  in  the  same  fashion  as  the  CW  signal  previously  described. 

The  code  used  in  the  coding  process  is  controlled  by  a  drifting  clock  (see  Sec.IV-D)  so  that, 
if  c(ca)  is  the  frequency  domain  description  of  the  coding  signal,  c(a;n)  exp[— icur^l  is  the  fre¬ 
quency  domain  description  for  the  n**1  delay  of  decoding  signal.  After  decoding,  the  frequency 
domain  description  of  the  signal  for  the  n^  delay  is 


(«)  =  Of  JJ  [no(«"-wo)  ®  c*(w" )]  ©  c*(w")s(Y,  Z) 


x  1WTn  -iwr(Y.Z)  dYdZ 
cos  0 


I s(- 


dZ 


AwF,  Z  )  R  [t  —  r(— AcuF,  Z  )] 
'  n  c  1  n  *  1 


n 


n  1  cos  0 


where  R  (r)  is  the  autocorrelation  function  of  the  code, 
c 


(46) 

(47) 
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This  expression  assumes  the  planet's  Doppler  spread  is  less  than  the  code  repetition  rate, 
as  it  neglects  the  frequency  foldover  effects  that  arc  produced  when  this  is  not  the  case.  Thus, 
the  output  of  the  delay  channel  and  c^  frequency  channel  is 


[  I  2  2  2/cU  ^  c  L  , 

x  =  r  a  (S  „+S 

c, n  p  c,  n  c,  n  ’ 


and  the  output  of  the  "cross  channel”  is 


Z  2 

x  y*  =  r  a  (3 
c,nJc,n  p 


\  c,n 


2  7ri  Z  ( 

e  n  Z  +  SL  e 
c,n 


—  2 7r i  |  zn  | 


( 48) 


(49) 


where  S c  is  the  scattering  function  integrated  over  the  frequency  (resolution  along  Y  )  and 

delay  (resolution  along  X  )  resolutions.  Since  the  delay  resolution  is  ambiguous  in  the  sign  of 
Z  ,  the  superscripts  U  and  L  signify  the  upper  and  lower  hemisphere  points.  If  the  frequency 
channel  number  is  measured  from  the  zero  frequency  or  frequency  of  the  return  from  the  sub¬ 
radar  point,  and  the  delay  channel  number  is  measured  from  the  delay  that  just  grazes  the  sub¬ 
radar  point. 


where 


U 

x,n 

L 

x,n 


Y  2±Z  2 
n  n 


if 


S(Y  ,  Z  )  dY  dZ 
_ n  n  n  n 

cos  0 


(50) 


Y  1±Z  1 
n  n 


,,  .  ,  1  .  6oj  ^ 

Y  1  =  (c  —  -  F 

n  2  27 r 


Y  2=  (c  +  i)  ^  F 
n  2  2  7r 


(51) 


and 


X  1  =  1  - 
n 


(n  —  1 )  6t2c 


X  2  =  1  - 
n 


n6i-2c 


(52) 


Zn1  =Jl  -[(Xn1)2  +  (Yn1,2l 
Zn2  =  71  -  «V)2  +  (Yn2,2l 


(53) 


ED.  METHODS  OF  MAPPING  SURFACE  OF  PLANET 

A.  One-Dimensional  Interferometry  to  Reconstruct  Received  Power 
Distribution  Along  a  "Doppler  Strip" 

It  was  shown  in  Sec.  II  that,  when  a  CW  signal  is  transmitted,  the  echo  received  in  a 
given  frequency  channel,  after  correction  for  the  Doppler  shift  to  the  subradar  point  (Doppler 
tracking),  arises  from  a  strip  on  the  face  of  the  planet  toward  the  radar.  While  the  power  in  a 
given  frequency  channel  for  a  single  antenna  represents  an  integration  along  the  entire  Doppler 
strip,  the  crosscorrelation  power  between  the  two  antennas  of  the  interferometer  yields  a  Fourier 
component  of  the  distribution  along  the  strip.  After  we  correct  for  the  planet's  motion  relative 
to  the  interferometer,  the  cross  power  is  given  by  Eq.(44)  which  can  be  rewritten 
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(54) 


Ac^Z^ “ 


I 


S(c,  Z  )  e 

\  ,  n/ 


2?riZ  t 
n  Z 


dZ 

_ n_ 

cos  0 


A  (i7)  is  the  complex  fringe  amplitude  (x  y*)  for  projected  baseline  i  omitting  the  attenu- 
e  z  c  c  z 

ation  constants.  Equation  (54)  can  be  inverted  by  one-dimensional  Fourier  theory 


S|c'Zn> 

eos  0 


I 


Ac(iz)  e 


-27riZ  !„ 
n  Z 


dl 


z 

27T 


(55) 


where  the  limits  of  integration  would  have  to  be  infinite  for  a  perfect  inversion.  If  the 


maximum 


projected  baseline  is 

J  Z  max 


,  then  the  resolution  pattern  is 


f 


Z  max  "27rl/',niZ  ^fz 
2tt 

Z  max 


sin  27rZ 

_ n  Z  max 

2ttZ  I- 

n  Z  max 


(56) 


or  a  half-power  width  of 


AZ 

n 


1 

21  „ 

Z  max 


(57) 


in  practice,  it  is  difficult  to  cover  a  wide  range  of  i ^  without  a  fortuitous  baseline  orientation 
to  make  full  use  of  changes  with  hour  angle  that  result  from  a  combination  of  those  ehanges  I 
and  I  illustrated  in  Fig.  2.  is  the  projection  onto  the  dashed  line  in  Fig.  2.)  However,  since 

S  is  real, 

A*(iz)  =  Ae(_iZ)  (58) 


so  that  only  the  magnitude  of  the  projected  baseline  is  important;  also,  it  ean  be  assumed  the 
refleeted  power  is  limited  to  the  planet  surface  so  that 

S(c,Zn)=0  (59) 

for  values  of  I  Z  ^  Z  ,  where 

1  n  n  max 


1  -  Y  (e) 
n 


(60) 


This  restriction  of  S(c,  Z^)  allows  Ae(^)  to  be  reconstructed  from  sampled  values.  From  the 
sampling  theorem,  (f ^ )  need  only  be  sampled  at  intervals  of 


AiZ  2Z 


(61a) 


If  only  the  Fourier  components  of  the  spatial  distribution  corresponding  to  baseline  projections 

between  and  are  measured,  then  the  equivalent  resolution  pattern  for  a  uniformly 

Z  min  Z  max  J 

weighted  transform  is 


cos  [7rZ  [Hr?  .  +  S-  )]sin[7rZ  (^v  —  .  )] 

1  nv  Z  min  Z  max  1  1  n  Z  max  Z  min 

ttZ  (!_  l7  .  ) 

n  Z  max  Z  min 


(61b) 


B.  Resolution  of  ” Range -Doppler  Ambiguity"  with  Two-Element  Interferometer 

When  the  transmitted  signal  is  eoded  so  that  single -antenna  measurements  resolve  the 
planet  into  range-Doppler  cells,  there  is  a  twofold  ambiguity,  that  is,  two  points  on  the  planet 
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have  the  same  range  and  Doppler  shift.  However,  the  cross-power  range-Doppler  cells  are 
the  vector  sum  of  contributions  from  the  ambiguous  points  as  shown  in  Eqs.(48)  and  (49). 


P(c,  n) 


su  +  s  L 

c,n  c,n 


(62) 


and 


A(c,  n) 


SU  e 
c,n 


icp(c,n,l7) 


+  S  L  e 
c,n 


-iip(c,n,l  ) 


(63) 


where  P(c,  n)  is  the  single-antenna  power,  and  A(c,  n)  is  the  cross  power  for  the  c'  frequency 
channel  and  n^  delay.  The  interferometer  phase  for  the  upper  or  northern-range  Doppler  cell 
is  <p{  c,  n,i^).  Equation  (63)  assumes  that  range-Doppler  cell  is  much  smaller  than  the  fringe 
spacing.  The  range-Doppler  mapping  geometry  is  illustrated  in  Fig.6(a-b).  Equations  (62)  and 
(63)  can  be  solved  for  and  S*J  provided  <p  is  not  a  multiple  of  2tt.  However,  if  measurements 
are  made  over  some  change  of  baseline  projection  a  least-squares  fit  to  the  data  can  be  per¬ 
formed,  thereby  eliminating  the  case  of  ip  being  a  multiple  of  2ir.  The  least-squares  solution 
to  Eq.  (63)  is 


u 

Sc,  n 

SL 
c,  n 


2  2 
Z  Re  A(c,  n)  cos  <^(c,  n,  ly )  Z  sin  Z  Im  A(c,  n)  sin  </?(c,  n,  l7 )  Z  cos  <p(c,n,t7) 


i  - 


2  2 

Z  sin  <p(c,n,l7)  Z  cos  <p(c,n, /„) 

t  ^  t  ^ 

fZ  Z 


(64) 


which  provides  the  best  estimate  of  the  ratio  /S^  while  the  single -antenna  measurements 
are  used  to  estimate  the  sum  of  and  S^. 


before  multiplicolion  by 


exp  [2wi  AujF/yJ  ) 

(°)  (b) 


Fig.  6.  (a)  Twa  regions  with  same  range  and  Dappler  shift;  (b)  their  locations 

in  interferometer  fringe  pattern. 

C.  Signal-to-Noise  Ratio  Analysis 

It  is  convenient  to  convert  received  power  into  temperature  and  to  assign  a  system  temper¬ 
ature  Tg  to  the  receiver  and  sky  background  noise.  From  Eq.(48), 
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(65) 


2  2  11  T 

P  a  PTrlSl  +S1j) 
c,n  I  p  c,n  c,n 

where  P^  is  transmitted  power.  If  the  planet  were  uniformly  rough,  then  it  could  be  described 
by  a  scattering  function  which  is  only  a  function  of  the  angle  of  incidence  O,  in  which  case 


„  2  2  S(0) 

P  a  P  r  '  a 

c,n  T  p  cos  O  c,n 


(66) 


where  a  is  the  area  of  the  range-Doppler  cells  in  the  Y  ,  plane.  Converting  to  temperature. 


Tjj(c,  n) 


G2\2P_r  2S(0)  a 

T  p  x  c,n 

64tt  X  4K  cos  0 
o 


(67) 


where  K  is  Boltzmann's  constant.  The  fractional  cross  section  relative  to  a  perfect  isotropic 
reflector  is 


f  S(0) 


- ^  dY  dZ 

cos  0  n  n 


(68) 


The  Westford  antenna  temperature  is 


T 


W 


riIGW 

gh 


(69) 


and  effective  system  temperature  for  the  cross  power  is 

ft 


Ts- 


Vv 


(70) 


From  the  theory  of  a  Gaussian  process,  the  rms  deviation  in  the  measured  temperature  is 


At 


(Ts„*T'>) 


(71) 


^rms  \  AfT 

where  Af  is  the  frequency  resolution,  and  T  is  the  integration  time.  The  ’’noise  threshold"  of 
the  cross  power  is 


'AT 


(Ts„  *  T„)  (Ts,„  +  Tw) 


T  +  T 
S  W I 

W  w/ 


WH  rm  s  nr&z 

V  At  I 

and  the  rms  deviation  of  the  measured  phase  is 


(72) 


Ae 


rms 


K  [  T»)  (Tsw  +  Tu  ) 

^  |twh|  -/AfT 


(73) 


IV.  DESIGN  AND  CONSTRUCTION 

A.  Radio  Frequency  and  Local  Oscillator  Systems 

The  radio  frequency  (RF")  portions  of  the  interferometer  are  shown  in  Fig.  7.  The  trans¬ 
mitter  is  coupled  to  the  right-handed  circular  polarization  port  of  Haystack,  while  both  receivers 
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are  coupled  to  the  orthogonal  polarization  (left-handed)  ports.  A  system  temperature  of  about 
60° K  was  achieved  at  Haystack  using  a  maser,  and  about  150°K  at  Westford  using  a  parametric 
amplifier. 

The  local  oscillator  (LO)  system  is  the  critical  part  of  the  interferometer,  since  phase  co¬ 
herence  and  stability  depend  on  the  LO1  s  relative  phase  noise  and  stability.  The  first  and  most 
critical  LOs  are  klystron  oscillators  (at  7710  MHz)  which  are  phase  locked  to  reference  signals 
derived  from  a  station  standard.  Figure  8  shows  a  block  diagram  of  the  LO  phase-lock  system 
in  which  the  klystron  is  servo  controlled  to  the  64  harmonic  of  120  +  30  MHz. 

B.  Intersite  Coupling 

In  order  to  maintain  constant  phase  difference  between  the  klystrons  at  the  stations,  it  was 
necessary  to  compensate  for  changes  in  the  electrical  length  of  the  cables  carrying  the  refer¬ 
ence  signals  between  the  two  stations.  The  transmission  line  servo  system  shown  in  Fig.  9  main¬ 
tains  constant  electrical  length  by  reflecting  some  of  the  120-MHz  reference  so  that  a  phase 
comparison  with  the  transmitted  signal  can  be  performed  and  used  to  provide  an  error  signal. 

In  order  to  identify  the  reflected  signal,  a  100-kHz  phase-reversal  modulation  is  applied  to  the 
reflected  signal.  Figure  10(a)  shows  the  mechanical  line  "stretchers"  used  to  correct  the  line 
length.  Figures  10(b),  (c),  and  (d)  show  the  Westford  receiver  and  parametric  amplifier. 

C.  Data  Interface  to  CDC  3300  Data- Processing  Computer 

The  final  intermediate  frequency  (IF)  of  both  the  Haystack  and  Westford  receiver  channels 
was  2  MHz.  The  CDC  3300  Direct  Data  Interface  converted  the  signals  into  a  form  suitable  for 
digital  processing  in  the  CDC  3300  data-processing  computer. 

The  IF  signals  were  first  band  limited  to  reduce  the  dynamic  range  requirements  of  the  signal 
processing.  Then,  each  channel  was  frequency  translated  to  a  0-Hz  IF  by  two  identical  balanced 
mixers  fed  by  quadrature  reference  signals  at  2  MHz.  The  resulting  pair  of  quadrature  video 
signals  represented  orthogonal  components  of  a  complex  signal.  These  video  signals  were  then 
filtered  by  an  identical  pair  of  low-pass  filters  that  produced  the  equivalent  effect  of  a  single 
bandpass  filter  at  IF. 

The  filters  used  for  the  CW  data  were  5-pole  low-pass  filters  with  cutoff  frequencies  at 
256  Hz,  which  defined  the  total  processing  bandwidth  of  512  Hz  (±256).  The  filters  used  for  the 
coded-pulse  data  were  rectangular  500-psec  pulse-matched  filters. 

In  the  coded-pulse  case,  the  planetary  Doppler  frequency  spread  was  used  to  limit  the  spec¬ 
tral  width  of  the  signal  data  processed.  Noise  data  were  unavoidably  aliased  by  the  subsequent 
sampling.  The  CW  signals  were  deliberately  offset  in  frequency  by  +100  Hz  so  that  a  symmetrical 
noise  comparison  bandpass  was  available  and  to  avoid  centering  the  signal  on  the  troublesome 
zero-frequency  point.  (DC  offsets  in  the  analog  system  show  up  as  false  zero-frequency  signals 
in  the  subsequent  processing.) 

Filtered  signals  were  amplified  and  simultaneously  sampled  by  a  4-input  analog  sample-and- 
hold  multiplexer  (these  four  signals  are  the  two  pairs  of  complex  signals  for  Haystack  and  West¬ 
ford  receivers).  The  "held"  signals  were  then  sequentially  encoded  by  an  8-bit  A/D  converter. 
Encoder  output  was  plaeed  in  the  proper  format  and  transmitted  directly  to  the  lower  16K  core 
memory  of  the  CDC  3300  data-processing  computer  by  the  Direct  Data  Interface. 
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Fig.  7.  RF  black  diagram. 


16 


HAYSTACK 


CLOCK  100->*MC  QUANTA 


Fig.  7.  Continued. 
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Fig.  8.  LO  phase-lack  system. 
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Fig.  9.  Transmission  line  servo  system. 
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Fig.  10(o).  Line  stretchers. 


Fig.  10(b).  Westford  receiver  at  Haystock. 


Fig.  10(c).  Westford  receiver  at  Westford. 


Fig.  10(d).  Close-up  view  of  Westfard 
porometric  omplifier. 
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The  Interface  is  under  timing  control  of  the  Radar  Sequencer,  which  generates  start  and 
stop  functions  as  well  as  supplying  range  rate  offset  sampling  pulses.  The  Interface  directly 
controls  the  location  in  core  memory  for  each  data  word,  and  periodically  communicates  control 
and  time  information  to  the  computer  program  by  a  standard  computer  communication  channel, 
to  synchronize  the  computer  program  cycling  to  the  radar  timing. 

The  CW  data  were  sampled  at  a  rate  of  512  samples  per  second  for  each  complex  pair  of 
samples,  covering  an  unambiguous  frequency  range  of  51 2  Hz.  The  coded-pulse  data  were 
sampled  for  31  points  at  500-psec  spacing  for  eaeh  eode  interval  of  15.5  msee.  Eaeh  range  box 
then  effectively  eovered  an  unambiguous  frequency  range  of  64.5  Hz,  which  produces  some  fre¬ 
quency  aliasing  of  the  noise  in  the  analog  signal  bandwidth  being  eneoded. 

D.  Doppler  and  Range  Rate  Tracking 

During  the  experiments,  the  radar  receiver  was  continuously  tuned  to  remove  the  Doppler 
frequency  shift  of  the  received  signal  arising  from  the  motion  of  the  planet's  eenter  with  respeet 
to  the  observer.  The  ranging  experiments  also  required  a  eloek  signal  that  was  offset  in  fre¬ 
quency  to  compensate  for  the  range  rate  of  the  moving  planet.  Equipment  used  for  tracking  both 
the  Doppler  shift  and  range  rate  was  under  control  of  the  Univae  490  pointing  computer.  The 
basic  frequency  control  information  was  derived  from  the  planetary  ephemeris. 

Doppler  frequency  was  removed  by  offsetting  the  second  receiver  LO  at  100  MHz  by  the 
amount  of  the  predicted  Doppler  shift.  This  LO  signal  was  generated  using  a  digitally  controlled 
frequency  synthesizer  (MP5100A).  The  frequency  was  incremented  in  0.1 -IIz  steps  by  command 
of  the  U490  pointing  computer  at  a  rate  of  20  commands  per  second.  Since  the  typical  rate  of 
change  of  frequency  was  less  than  2  Ilz/sec,  the  Doppler  shift  was  compensated  with  a  precision 
of  ±0.1  Hz  and  with  an  accuracy  determined  by  the  basic  ephemeris  and  knowledge  of  time.  Sinee 
the  coherent  data  processing  did  not  attempt  an  equivalent  frequency  resolution  finer  than  1  Hz, 
the  frequency  tuning  system  itself  should  not  have  produced  significant  speetral  smearing  or  off¬ 
set.  Although  the  frequency  of  the  synthesizer  was  stepped  20  times  per  second,  phase  discon¬ 
tinuities  of  the  LO  were  produced  only  by  1-kIIz  eontrol  digit  changes  which  occurred  typically 
at  15-minute  intervals. 

The  equipment  used  for  generating  the  range  rate  offset  eloek  signal  took  advantage  of  the 
fact  that  the  fractional  frequency  offset  required  is  the  same  as  the  ratio  of  Doppler  frequency 
shift  to  carrier  frequency.  The  carrier  frequency  was  7840  MHz  and  the  basic  eloek  was  derived 
from  the  1-MHz  station  standard  frequency.  The  generated  Doppler  frequency  offset  present  in 
the  seeond  LO  signal  conceptually  was  divided  by  7840  and  applied  as  a  single-sideband  modula¬ 
tion  to  the  standard  1-MHz  eloek  signal,  preserving  sign  as  well  as  magnitude.  The  aetual  hard¬ 
ware  implementation  was,  in  fact,  more  complex  but  did  produce  an  offset  clock  signal  that  was 
used  to  sample  the  received  signal  for  digital  processing.  In  effeet,  a  moving  time  base  that 
tracked  the  planetary  motion  was  used  to  produec  a  "stationary"  target  with  respect  to  the  com¬ 
puter  data  processing. 

E.  System  Sequencer 

The  timing  signals  for  the  Rayford  experiments  are  controlled  by  the  Haystack  Radar  Se¬ 
quencer.  All  experiments  start  on  a  selected  minute.  Transmission  is  then  enabled  for  an 
interval  equal  to  the  round-trip  signal  flight  time  (received  from  the  target  ephemeris  and  quan¬ 
tized  in  1-fisec  increments). 
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After  the  flight  time  has  elapsed,  control  lines  are  activated  which  signal  the  changeover 
from  transmit  to  receive.  The  signal  A/D  converter  and  computer  interface  are  enabled  in  prep¬ 
aration  for  the  data-taking  phase. 

On  the  occurrence  of  a  preset  time  pulse  (about  3  0  sec  later),  sample  commands  are  gen¬ 
erated  and  transferred  to  the  A/D  control  unit.  The  incoming  return  is  sampled,  the  resulting 
8-bit  word  is  formatted,  and  the  characters  are  transferred  into  the  computer's  storage  bank. 
Samples  are  taken  continuously  until  another  flight-time  interval  has  elapsed,  at  which  time 
stop  commands  are  generated,  control  is  transferred  to  the  computer  program,  and  the  system 
is  readied  for  another  run. 

All  timing  during  the  receive  interval  is  based  on  a  Doppler-corrected  1-MHz  clock  which 
compensates  for  the  time  compression  (expansion)  of  the  target  return.  The  sample  pulse  gen¬ 
erators  and  interval  counters  thereby  track  the  return  in  time.  System  timing  functions  are 
generated  relative  to  the  initial  start  with  a  precision  of  ±125  nsec  and  to  an  accuracy  of  1  psec. 
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F.  System  Tests 

Assuming  that  both  systems  already  operate  as  s ingle -antenna  planetary  radars,  the  major 
tests  are  those  concerned  with  testing  the  opera¬ 
tion  of  the  system  as  an  interferometer.  The 
first  test  consisted  of  running  both  receivers  ad-  -g 

jacent  to  each  other  (at  Haystack),  with  the  in-  jj{ 

< 

ter  site  coupling  cables  making  a  loop  through  the  £ 

Millstone  radar  (approximately  half  way  between 
Haystack  and  Westford).  By  feeding  a  test  sig¬ 
nal  into  both  receivers,  the  relative  phase  sta¬ 
bility  of  the  interferometer  receivers  could  be 

checked  directly.  Figure  11  shows  the  long-term  phase  stability,  while  Figs.  12(a)  through  (c) 
indicate  short-term  stability  or  phase  noise  of  the  system. 

Further  phase  stability  tests  were  made  with  the  Westford  receiver  installed  on  the  60-foot 
antenna.  These  tests  were  first  made  with  a  test  transmitter  located  on  a  pattern  truck  50  km 
away.  Large  relative  phase  fluctuations  were  observed,  as  shown  in  Fig.  13;  however,  these 
fluctuations  were  shown  to  be  largely  due  to  refractive  index  fluctuations  by  repeating  the  test 
with  a  test  transmitter  at  Haystack  (Fig.  14). 


Fig.  11.  Instrumental  phase  stability. 


V.  DATA- REDUCTION  TECHNIQUES 

A.  Atmospheric  Fluctuations  and  Refraction 

When  computing  the  fringe  pattern  or  the  position  of  the  lines  of  constant  phase  relative  to 
the  celestial  sphere,  it  is  necessary  to  take  careful  account  of  the  refraction  in  the  earth's 
atmosphere  and  ionosphere.  If  the  earth's  atmosphere  and  ionosphere  were  plane  parallel  and 
uniform  along  the  plane,  then  the  ray  paths  reaching  two  antennas  at  equal  elevations  are  equal, 
so  that  the  interferometer  phase  ojt^  could  be  computed  from  Eq.(5).  In  other  words,  in  this 
case  the  baseline  can  be  considered  to  be  outside  the  atmosphere  since  the  atmosphere  does  not 
affect  the  fringe  pattern  location.  Because  of  the  curvature  of  the  earth  and  the  nonuniformity 
of  the  atmosphere,  a  precise  estimate  of  the  location  of  the  fringe  pattern  can  only  be  made  by 
integrating  the  electrical  length  along  the  ray  path  to  each  antenna.  However,  the  effect  of  earth 
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Fig.  12(o).  Phosor  of  Haystack  (top  trace) 
and  relotive  phase  (lower  trace)  with  bright 
spot  showing  zero-frequency  offset. 
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Fig.  12(b).  Phose  naise  af  Haystock  and  Westfard  signals  relative  ta  station  stondard 
(10-sec  exposure  of  2-MHz  signals  with  the  station  standard  1  MHz  as  scope  sync). 


Fig.  12(c).  Relotive  phase  naise  (10-sec  exposure  of  Hoystock  2-MHz  IF 
with  Westford  IF  as  scope  sync). 


HOURS 

Fig.  13.  Differential  phose  between  Westfard 
ond  Haystack  measured  using  test  signol  trans¬ 
mitted  from  pottern  truck  50  km  owoy. 
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Fig.  14.  Phase  between  Westford  and  Hoystack 
far  one-woy  transmission  of  test  signol  from 
Hoystock  to  Westford. 
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curvature  can  be  taken  into  account,  when  the  antennas  are  spaced  only  a  few  miles  apart  or  less, 
by  using  Eq.(4).  The  apparent  elevation  E^  differs  from  the  elevation  without  the  atmosphere 
by  the  refractive  bending  e.  When  E^.  and  d  are  zero,  Eqs.  (4)  and  (5)  give  the  same  result  if 
e  is  approximated  by 


e  »(N  x  10"6)  cotE 

P 


(74) 


where  N  is  the  refractive  index  on  the  ground  in  N -units.  More  precise  bending  values,  together 

4 

with  refractive  index  values,  have  been  tabulated  by  Bean  and  Dutton.  At  X-band,  the  ionospheric 
refractivity  is  only  a  few  N-units  (owing  to  a  decrease  proportional  to  frequency  squared),  while 
the  major  contributions  to  the  tropospheric  refraction  are  water  vapor  ( — 100  N-units  for  90- 
perccnt  relative  humidity)  and  the  dry  atmosphere  (  —  220  N-units  at  25 °C).  Unfortunately,  the 
contribution  of  the  water  vapor  to  the  refraction  is  liable  to  be  very  nonuniform  and  can  result 
in  large  relative  path  differences  between  the  two  antennas.  The  relative  electrical  path  differ¬ 
ence  which  produces  a  phase  shift  in  the  fringe  pattern  of  27tANAD/a  radians  is  time  varying  and 
results  in  "bad  seeing."  For  example,  a  cloud  of  15  N- 
units,  3000  feet  deep,  will  produce  a  phase  shift  of  2ir 
radians  over  a  time  span  of  a  few  minutes  as  it  drifts 
through  the  beam  of  one  and  then  the  other  antenna. 

Clouds  are  not  the  only  form  of  atmospheric  irregular¬ 
ities  —  weather  fronts  and  turbulence  also  produce  bad 
seeing  conditions.  Figure  13  shows  the  interferometer 
phase  fluctuations  produced  by  relative  path  differences 
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Fig.  15.  Interferometer  phase  of  region 
around  subrodar  point. 


of  a  signal  from  a  pattern  track,  while  Fig.  15  shows  the  phase  fluctuation  of  the  subradar  point 
on  Venus  as  measured  with  the  planetary  radar  interferometer.  Since  the  subradar  region  on 
Venus  is  constrained  by  the  range-Doppler  gate,  it  is  difficult  to  see  how  those  fluctuations  could 
arise  from  Venusian  atmosphere  rather  than  the  earth's  atmosphere. 


B.  Use  of  Subradar  Point  as  Phase  Calibrator 

As  a  result  of  the  refractive  fluctuation  discussed  in  Sec.  A  above,  it  was  desirable  to  stabi¬ 
lize  the  fringe  pattern  relative  to  the  planetary  surface.  This  was  done  by  referring  all  the 
measured  phases  to  that  of  the  first  range  box  and  zero  Doppler  on  the  cell  which  surrounds  the 
subradar  point.  In  the  "coded"  mode,  the  maximum  error  of  the  phase  calibration  due  to  the 
energy  within  the  range-Doppler  cell  being  nonuniformly  distributed  is  just  the  phase  difference 
between  the  subradar  point  and  the  edge  of  the  cell.  Thus,  it  is  most  desirable  to  make  the  first 
range-Doppler  cell  as  small  as  possible  by  making  the  first  range  box  just  graze  the  leading 
edge  of  the  planet.  Figure  16  illustrates  the  calibration  box  and  the  fringe  pattern  on  the  planet. 
In  the  CW  mode,  the  phase  calibration  was  performed  on  the  zero-Doppler  strip.  This  technique 
could  have  resulted  in  a  large  error,  as  the  phase  of  this  strip  could  be  far  removed  from  that 
of  the  subradar  point.  However,  because  of  the  strong  highlight  at  the  subradar  point  or  the 
sharpness  of  S(0)  about  0=0,  the  phase  is  essentially  determined  by  the  return  from  the  sub¬ 
radar  point. 

C.  Removal  of  "Mean  Planet" 

The  planet's  scattering  law  S(Yn,  Z^),  which  is  defined  as  the  power  returned  from  unit  area 
of  the  planet  for  an  incident  plane  wave  of  unit  power  per  unit  cross  section,  can  be  decomposed 
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Fig.  16.  Calibration  cell. 


into  a  "mean  planet"  scattering  law  S(G)^  and  deviation  from  the  mean  planet  S^( Y^,  Z^).  The 
S(O)  is  only  a  function  of  the  angle  of  incidence  G,  while  the  deviation  results  from  surface 
"features."  From  Eq.(50),  the  power  in  a  particular  delay-Doppler  cell  in  the  projected  Y  , 
plane  is 
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S(«) 

cos  0 


dY  dZ  + 
n  n 


I 


S.(Y  ,  Z  ) 
n'  n 


dY  dZ 


n  n 
cos  G 


(75) 


where  the  integration  is  carried  over  the  area  of  the  cell.  In  the  CW  mode,  the  power  in  each 
Doppler  strip  from  the  mean  planet  scattering  law  is 


S(Q) 
cos  G 


dZ 

n 


(76) 


if  the  Doppler  strips  are  very  thin,  while  the  power  in  each  range  cell  before  frequency  analysis 
is 


S(O) 

cos  O 


X  (area  of  Doppler  cell) 


(77) 


or  just  S(G)  since  the  area  is  proportional  to  cos  G.  Since 

r  —  2tc 

cos  G  =  X  =  _ E - 

n  r 

P 


(78) 


0  is  determined  by  the  delay. 

The  above  decomposition  into  mean  planet  and  features  is  especially  useful  when  the  mean 
planet  dominates  the  signal  returned.  This  is  the  case  for  Venus,  which  has  a  scattering  func¬ 
tion  S(G)  (for  the  polarized  echo)  which  is  highly  peaked  at  G  =  0  and  drops  by  more  than  30  dB 


$  S(0)  includes  effect  of  the  planet's  atmospheric  attenuation. 
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at  the  limb  or  as  0  -  90°.  The  power  in  each  range-Doppler  cell  at  angle  of  incidence  0  is 
approximately 


S(0) 
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dY  dZ 


n  n 

cos  0 


(79a) 


where 
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and 


cos  0  =  X  =  1  - 

n  2r 

P 


(79c) 


If  power  in  each  cell  after  subtraction  of  the  mean  noise  is  divided  by  the  mean  scattering  law, 
then  the  power  in  each  cell  will  be  proportional  to 


1  + 


SHY  ,  Z  ) 
<3V  n  n 

S(0) 


cell  area 
cos  O 


(80) 


or  a  planet  with  no  features  will  appear  uniformly  bright.  A  similar  scattering  law  compensation 
was  performed  on  the  CW  data,  where  the  power  in  each  Doppler  strip  was  divided  by  the  in- 
tegrated  mean  [Eq.  (76)).  Also,  the  correlation  function  of  the  mean  planet 


dZ 

f  S(O)  COS  27rZ 

n  Z  cos  0 
dZ 

/  S( 0)  - 7T 

J  '  cos  0 


(81) 


was  subtracted  from  the  power  so  that  only  the  features  should  appear  in  the  transformation. 
This  was  done  to  eliminate  the  sidelobes  of  the  mean  planet  that  result  from  transformation 
[ Eq.  (55)]. 


D.  Real-Time  Computer  Programs 

The  purpose  of  the  real-time  radar  signal-processing  program  is  to  collect  the  sampled 
data  of  the  received  radar  signal  and  to  perform  as  much  processing  of  these  samples  as  the 
computer  speed  and  rate  of  signal  input  will  allow.  Finally,  the  results  of  this  processing  must 
be  saved  (in  this  case  digitally  recorded)  for  any  subsequent  ” off-line"  analysis. 

For  the  Ilayford  interferometer  experiment,  two  real-time  processing  programs  were 
written  for  the  two  basically  different  experiments  that  were  attempted:  the  continuous  wave, 
and  the  coded-pulse  transmission,  respectively.  These  programs  were  written  for  the  CDC  3300 
computer  at  the  Haystack  Microwave  Facility.  The  CW  and  codcd-pulsc  programs  are  described 
separately  below. 
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Fig.  17.  Display  of  spectra  from  real-time  CW  program. 


1.  C\Y  Program 


Prior  to  a  day's  operation,  the  relative  phase  cot^  and  the  rate  of  change  of  relative  phase 
of  the  two  received  signals  were  computed  from  elements  of  the  planetary  epheme rides  and  the 
known  interferometer  geometry.  Provisions  for  inputs  relative  to  the  expected  atmospheric 
conditions  expected  were  also  provided  for  this  calculation.  These  elements  of  phase  rotation, 
in  intervals  of  1  minute,  were  stored  on  magnetic  tape  and  used  as  an  input  to  the  real-time 
program.  During  the  transmit  portion  of  a  run,  the  starting  time  of  the  receive  period  is  input 
to  the  computer,  the  tape  is  searched,  and  the  phase  correction  components  are  interpolated 
and  stored  in  memory  for  every  coherent  interval  (approximately  1  sec).  Since  512  samples 
are  taken  within  each  coherent  interval,  the  phase  corrections  for  each  sample  are  computed 
by  interpolation  in  real  time  by  means  of  the  following  recursion  formulas:  . 

cose.,  ,  cos  </>  cosAcp  —  sin  sin  Ac7  (82) 

K  + 1  K  r\ 


sin  -  sin  </>  cosAc 

K+l  k 


+  cos  sin  A(p 


(83) 


where  c p  is  the  phase  rotation,  and  A <r  is  the  phase  rotation  increment.  This  phase  correction 
K 

is  then  applied  to  the  proper  signal.  Subsequent  to  the  phase  correction,  the  frequency  spectra 

5 

of  the  signals  are  computed  by  means  of  a  fast  Fourier  transform  routine  which  produces  com- 
plex  spectra  components  for  512  1-Ilz  filters  which  are  then  recorded  on  magnetic  tape.  Suf¬ 
ficient  time  remains  after  this  operation  to  produce  incoherent  sums  of  the  two  frequency  spectra 
which  are  displayed  at  the  end  of  each  run  as  shown  in  b  ig.  17. 


2.  Coded-Pulse  Program 

Since  the  modulation  applied  to  the  transmitted  radar  waveform  is  simply  a  phase -reversal 
type,  the  decoding  process  employed  by  the  program  consists  of  a  pattern  of  additions  and  sub¬ 
tractions  of  the  data  samples  identical  to  the  pattern  of  phase  shifts  impressed  on  the  trans¬ 
mission.  The  coding  selected  for  this  experiment  was  a  31 -element  shift  register  sequence 
with  a  baud  length  of  500jjisec.  The  decoding  process  yields  31  range  sample  components  for 
each  signal,  with  a  sample  interval  equal  to  the  baud  length.  At  the  start  of  each  receive  period, 
the  time  to  lOOgsec  is  read  into  the  computer  from  the  station  digital  clock.  This  time  is  re¬ 
corded  on  the  output  tape  along  with  the  decoded  range  samples,  since  it  is  necessary  for  com¬ 
puting  the  phase  corrections  in  later  analysis.  Sufficient  time  remains  after  the  decoding  and 
recording  to  incoherently  sum  the  power -vs-delay  functions  for  each  signal  for  display  at  the 
conclusion  of  the  run  as  shown  in  Figs.  18(a)  and  (b). 

Listings  of  the  real-time  programs  appear  in  the  Appendix. 


E.  Data- Averaging  Programs 

In  the  C\V  mode,  the  complex  samples  for  each  second  are  used  to  obtain  the  power  and 
cross  powers 


PH<W)  = 

Xjj(ct;)  xj^(co) 

(84) 

pW<«>  = 

xw(w)  x\iv(a,) 

(85) 

'iiw((i)  = 

Xj  j(Cc’  )  X*y(0D 

(86) 
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Fig.  18.  Real-time  display  of  power  vs  delay:  (a)  Haystack,  (b)  Westford 
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Fig.  18.  Continued 


No.  of  Runs 

(~5  minutes  each) 

15 

23 

18 

9 

13 

10 

20 

13 

10 

13 

Longitude 
of  Subrodar 
Point  Longp 
(deg) 

—  80.5 

-68.3 

-60.5 

-59.3 

-53.0 

-47.6 

-46.8 

-40.5 

-39.6 

-33.1 

O 

a  o 

3  c:  — J  cn 
+-  _q  .  a) 

3.8 

5.5 

6.8 

7.0 

8.0 

8.7 

8.7 

8.7 

8.7 

8.0 

F 

(center-to-limb 

Doppler) 

(Hz) 

107 

88 

73 

70 

58 

52 

52 

59 

62 

75 

D 

(angle  between 
apparent  rototion 
axis  ond 
celestial  north) 
(deg) 

30 

32 

32 

31 

28 

25 

24 

17 

16 

13 

Range  of  Baseline 
(wavelengths  per 
rodius  along  the 
apparent  rotation  oxis) 

2.5  to  3.2 

2.4  to  3.7 

2.2  to  4.0 

2.3  to  3.9 

1.85  to  3.3 

1.95  to  3.75 

2.3  to  4.1 

2.7  to  3.8 

2.55  to  2.9 

2.45  to  3.45 

Day  No. 
(1967) 

214 

222 

228 

229 

235 

241 

242 

249 

250 

256 
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which  are  then  summed  for  every  15  sec.  The  15-sec  averages  of  the  cross  powers  are  rotated 
by  the  phase  of  the  zero-Doppler  channel  (the  strip  through  the  subradar  point)  to  remove  the 
atmospheric  fluctuations  as  described  in  Sec.  V-B.  Further,  the  cross  power  is  rotated  by 
exp  [ZTiAcoFiy]  to  make  the  fringes  perpendicular  to  the  Doppler  axis  [see  Eq.  (41)],  and  the 
averages  are  taken  for  the  whole  run  —  a  run  being  a  complete  round-trip  transmit-receive  se¬ 
quence  of  about  5  minutes  for  the  conjunction  period. 

The  coded  data  are  rotated  and  averaged  in  the  same  manner  after  the  range  boxes  have 
been  spectrum  analyzed;  64  time  samples  are  spectrum  analyzed,  yielding  1.008-Hz  resolution. 

F.  Coordinate  Transformation  and  Data  Display 

For  each  day’s  observations,  the  projected  baseline  changes  in  a  manner  similar  to  that 
illustrated  in  Fig.  2.  In  the  CW  mode,  the  cross  power  is  transformed  to  give  the  distribution 
along  a  Doppler  strip,  at  which  point  the  data  are  in  the  (X,  Y,  Z)  coordinate  system  which  has 
subradar  point  at  (r^,  0,  0)  and  Z-axis  along  Doppler  axis.  The  data  in  each  cell  (1  Hz  in  Doppler 
and  ~l/l0  maximum  fringe  spacing)  are  transformed  to  a  system  (X^,  Y^,  Z^)  which  is  fixed  to 
the  planet's  surface.  Details  of  the  transformation  were  presented  in  Sec.  111. 

The  coded  data  are  first  transformed  from  the  range-Doppler  cells  of  the  upper  (northern) 
and  lower  (southern)  hemispheres  to  a  raster  in  (Y,  Z)  coordinates  and  then  to  (Y^,  Z^)  coordinates. 

Details  of  the  transformation  were  given  in  Sec.  II.  The  coded  data,  which  now  exist  as 
functions  of  X^  and  Y  ,  are  also  transformed  into  functions  of  planetary  latitude  and  longitude 
by  relations  which  again  are  the  inverse  of  Eqs.  (8)  through  (16).  In  addition,  the  ambiguity 
resolving  function  is  applied  to  the  data  by  multiplying  each  matrix  point  by  S^  q  or  n  of 
Eq.  (63)  according  to  whether  the  value  of  Z^  is  positive  or  negative.  In  this  way,  an  81  x  81  ma¬ 
trix  of  points  is  generated  in  a  square  planetographic  grid,  with  the  interval  between  each  matrix 
point  being  equal  to  1.25  planetary  degrees.  Spatial  smoothing  is  obtained  by  actually  computing 
the  power  every  5/12  of  a  degree  and  setting  the  matrix  point  equal  to  the  average  of  nine  of 
these  values  taken  from  a  3  X  3  array  centered  on  the  matrix  point. 

The  data  are  now  quantized  to  45  discrete  levels  and  transferred  to  an  intensity  plot  by  gen¬ 
erating  on  a  very  accurate  CRT  display  a  raster,  each  point  of  which  has  a  duration  directly 
proportional  to  the  value  of  the  matrix  point.  This  intensity  plot  is  photographed  by  a  special 
CRT  camera. 

VI.  RESULTS  OF  VENUS  MAPPING  FOR  1967  CONJUNCTION 

Venus  was  observed  during  the  period  2  August  through  13  September  1967  for  several  hours 
a  few  days  a  week.  Most  of  the  data  were  taken  in  the  CW  mode  because  of  the  larger  baseline 
coverage  needed  to  resolve  the  Doppler  strips.  However,  during  the  period  around  conjunction 
(30  August),  about  half  the  observing  time  was  spent  in  the  coded  mode. 

A.  Map  from  CW  Transmissions 

Because  of  limitations  in  the  operating  time  and  in  the  orientation  of  the  Haystack-Westford 
baseline,  it  was  not  possible  to  obtain  data  over  a  wide  range  of  projected  baseline.  Con¬ 
sequently,  determination  of  the  distribution  of  the  power  along  a  Doppler  strip  suffers  rather 
seriously  from  sidelobe  effects.  Table  1  lists  the  days  on  which  Venus  was  observed,  as 
well  as  the  range  of  the  component  of  the  projected  baseline  in  the  direction  of  the  apparent 
rotation  axis  covered  in  each  case.  On  many  days,  the  baseline  change  is  so  small  that  the 
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Fig.  19.  Resolution  patterns  for  doys  242  ond  256. 

equivalent  spatial  resolution  pattern  differs  little  from  a  cosine  function.  The  equivalent  res¬ 
olution  patterns  along  a  Doppler  strip  for  days  242  and  256  are  shown  in  Fig.  19.  Due  to  the 
inability  of  the  interferometer  uniquely  to  resolve  a  planetary  feature  because  of  the  large  side- 
lobe  effect,  it  is  advantageous  to  use  the  change  in  the  subradar  point  on  the  planet  as  well  as 
the  change  in  direction  of  the  rotation  axis  as  a  means  of  distinguishing  the  mainlobe  from  the 
sidelobes.  As  the  planet  rotates,  the  angle  between  the  Doppler  axis  and  the  planetary  axis 
changes  so  that  the  locations  of  the  interferometer  sidelobes  move  with  respect  to  the  mainlobe 
as  referred  to  the  planetary  surface.  Thus,  if  a  feature  is  inadvertently  associated  with  the 

position  of  a  sidelobe,  it  will  appear  to  move  from 
day  to  day,  while  correctly  placed  features  remain 
fixed.  However,  rotation  of  the  planet  is  suffi¬ 
ciently  slow  (~1°  per  day)  that  this  effect  will  not 
be  seen  in  one  day's  observations.  Figure  20 

shows  a  map  of  the  features  on  Venus  in  the  co- 

1 

ordinate  system  defined  by  Carpenter.  The  map 
represents  the  average  of  all  the  data  for  the  entire 
observing  period  placed  on  a  common  coordinate 
system;  the  amount  of  data  for  the  other  side  of 
the  planet  is  insignificant,  and  consequently  is  not 
displayed.  Figure  21  shows  the  same  map,  but 
indicates  (by  means  of  enclosures)  those  strong 
features  which  can  be  associated  with  the  main¬ 
lobe  of  the  resolution  pattern.  The  test  for  these 
features  was  made  by  a  search  through  the  maps 
obtained  for  each  individual  day.  The  fact  that 
some  of  the  scattering  centers  do  not  appear  to  move  during  the  observing  period  (days  214  to 
256)  lends  confidence  to  the  planetary  rotation  parameters  assumed. 

Figure  22  shows  the  lines  of  constant  Doppler  shift  which  pass  through  the  enclosed  features. 
The  circles  indicate  the  sidelobe  positions  for  these  features  for  different  days  on  which  ob¬ 
servations  were  made.  However,  some  days  have  been  omitted  for  clarity.  The  dotted  line 
indicates  the  position  of  the  subradar  point  during  the  observation  period. 
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Fig.  21 .  CW  map  with  features  that  remain  stationary  enclosed. 


Fig.  22.  Lines  af  constant  Dappler 
and  sidelabe  patterns  for  enclosed 
features. 
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TABLE  II 

CODED  OBSERVATION  DATES 

Mean  Angulor  Distance 

from  Subradar  Point 

for  S/N  Rotio 

b  O) 

29 

24 

29 

35.5 

35.5 

38 

38 

39 

35.5 

32 

b  '~5> 
2  3 

20 

18 

20 

26 

26 

27 

26 

34 

24 

22 

No.  of  Runs 
(~5  minutes  each) 

CO< —  CNCNCN'O'^CO'O'O 

Longitude 
of  Subrodor 
Point 
(deg) 

-78.4 

-60.2 

-59.1 

-53.1 

-51.2 

-47.8 

—46.9 

-46.1 

-40.6 

-33.0 

Latitude 
of  Subrodar 
Point 
(deg) 

3.9 

6.6 

6.8 

7.8 

8.0 

8.3 

8.4 

8.4 

8.2 

7.5 

D 

(angle  between 
apparent  rotation 
axis  ond 
plonetory  north) 
(deg) 

7.5 

10.4 

10.2 

9.0 

8.5 

4.8 

3.8 

2.8 

-3.1 

—6.4 

F 

(center-to-limb 

Doppler) 

(Hz) 

107 

74 

72 

59 

56 

53 

53 

53 

59 

74 

Day  No. 
(1967) 

215 

228 

229 

235 

237 

241 

242 

243 

249 

256 

B.  Maps  from  Coded  Transmissions 

Maps  made  from  the  results  of  the  coded  mode  of  operation  do  not  suffer  from  the  sidelobe 
effects  seen  in  the  CW  data,  but  are  not  entirely  free  from  ambiguities.  First,  the  upper-lower 
hemisphere  ambiguity  is  well  resolved  by  the  interferometer  only  when  the  signal-to-noise  ratio 
is  sufficient  to  produce  small  errors  in  the  measured  phase  of  the  cross  power.  Second,  the 
repetitive  31 -element  code  is  ambiguous  in  delay,  since  the  planet  is  81  code  elements  deep 
(40.5  msec).  Also,  there  is  a  certain  amount  of  frequency  foldover,  since  the  Doppler  spread 
at  the  later  delays  was  generally  larger  than  the  unambiguous  frequency  window  of  64.5  Hz.  The 
regions  of  frequency  foldover  have  been  eliminated  from  the  maps.  However,  the  regions  which 
are  aliased  in  delay  cannot  be  removed  and  some  features  may,  in  fact,  be  associated  with  regions 
which  lie  15.5  msec  deeper  into  the  planet.  In  general,  however,  the  rapid  falloff  of  power  in 
the  scattering  law  renders  this  possibility  unlikely. 

Figures  23(a)  through  (j)  are  intensity  plots  made  from  the  power  received  at  Haystack  as 
a  function  of  delay  and  Doppler.  In  these  plots,  the  two  ambiguous  hemispheres  are  given  equal 

weight.  Thus,  there  is  symmetry  about  the  line  Z  =0.  Figures  24(a)  through  (j)  are  intensity 

U  l  n 

plots  with  the  cross  powers  n  and  Sc  n  used  to  redistribute  the  power  between  the  two  hemi¬ 
spheres.  Figure  23(k)  shows  the  average  of  data  from  the  whole  experiment  without  hemispheric 
resolution,  while  Fig.  24(k)  shows  the  corresponding  average  using  the  ambiguity  resolution 
functions.  Table  II  lists  the  number  of  runs  for  each  day  and  the  approximate  distances  from  the 
subradar  point  for  which  the  Haystack  signal-to-noise  ratio  is  both  10  and  3  times  the  standard 
deviation  of  the  associated  noise  fluctuations  (a).  In  viewing  the  maps  in  which  the  cross  power 
has  been  used,  it  should  be  noted  that,  owing  to  the  smaller  antenna  size  and  higher  system 
temperature  at  Westford,  the  weighting  functions  used  to  redistribute  the  Haystack  power  become 
very  noisy  for  regions  where  the  Haystack  signal  alone  is  less  than  about  10  a. 

The  basic  spatial  resolution  is  usually  set  by  the  delay-Doppler  cell  size,  as  given  in 
Eq.(79b).  Reference  to  Fig.  16  shows  that  these  cells  become  larger  near  the  Doppler  equator. 

As  an  example,  on  Day  228,  F  =  74  so  that  at  delay  box  31  (about  49  planetary  degrees  from  the 
subradar  point)  the  width  of  the  cell  is  82  km  and  the  length  is  about  99  km.  In  comparison  with 
this,  the  cell  bordering  the  Doppler  equator  is  approximately  83  by  900  km. 


C.  Interpretation  of  Results 

Mean  Planetary  Scattering  Law:—  A  theoretical  scattering  law  derived  by  Muhleman^ 
was  found  to  fit  the  data  extremely  well  with  the  addition  of  3-dB  one-way  attenuation  for  a  vertical 
path  through  the  Venusian  atmosphere.  The  law  used  to  fit  to  the  CW  data  was 


S(0)  = 


sec  0  ^3 

H  cos  0Q- 

3 

3in  0  +  a  cos  0) 


(87) 


where  the  exponential  includes  the  effects  of  atmospheric  attenuation.  The  best  fit  to  the  data 
was  obtained  using 


P  -  1.3  ±  0.3 


a  =  0.13  ±  0.2 


From  Muhleman*  s  theory,  the  value  of  a  yields  an  effective  mean  slope  of  about  7°  at  3. 8 -cm 
wavelength. 
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Fig.  23(a-k).  Maps  from  Haystack  coded  data. 
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Fig.  23(a-k).  Continued. 
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Fig.  23(a-k).  Continued. 
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Fig.  24(a-k).  Maps  with  coded  interferometer  data. 
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Fig.  24(a-k).  Continued. 
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Fig.  24(a-k).  Continued. 
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Surface  Features:—  From  Figs.  20  through  24,  a  number  of  the  planetary  surface  regions 
are  seen  to  exhibit  a  local  increase  in  the  scattered  power  as  compared  with  the  average  scatter¬ 
ing  behavior.  Local  reductions  also  exist,  of  course,  but  these  are  far  less  obvious  in  the  data. 
The  common  characteristic  of  the  enhancements  is  that  they  display  relatively  little  dependence 
on  angle  of  incidence,  i.e.,  they  contrast  far  more  prominently  with  their  mean-scattering  sur¬ 
roundings  when  located  far  from  the  subradar  point  than  when  near  it.  From  this  behavior,  it 
is  deduced  that  they  are  unlikely  to  be  composed  primarily  of  relatively  smooth,  preferentially 
oriented  inclined  facets,  but  rather  possess  substantial  roughness  at  the  scale  of  the  wavelength. 
Observation  of  their  depolarized  scattering  properties  would  obviously  be  desirable  to  confirm 
this  hypothesis,  but  such  observations  were  not  possible  with  the  instrumentation  available. 
Depolarized  observations  of  these  features  have  been  made  at  longer  wavelengths,  however, 

and  these  have  in  fact  behaved  in  a  way  which  is  characteristic  of  diffuse  scattering  from  a  locally 
rough  surface. 

The  question  of  the  actual  composition  of  the  anomalously  scattering  features  is  unresolved. 
As  discussed  in  the  next  section,  it  seems  clear  that  these  features  are  quite  firmly  attached  to 
a  rigid  surface.  But,  whether  they  merely  represent  patches  of  locally  rough  horizontal  surface 
or  correspond  to  rough  mountainous  regions  is  impossible  to  tell  from  the  present  data.  Long 
continued  observation  using  extremely  short  observing  pulsewidths  might  resolve  the  topography 
along  the  track  of  the  subradar  point  to  the  order  of  several  hundred  meters.  And  if,  as  seems 
likely,  there  exists  substantial  atmospheric  attenuation  at  centimeter  observing  wavelengths, 
it  is  possible  that  variations  in  the  height  of  the  subradar  point  could  affect  the  observed  radar 
cross  section. 

Location  of  Features  and  Rotation  Period  of  Venus:—  The  system  of  planet-centered 

coordinates  used  in  the  reductions  reported  here  is  based  on  a  sidereal  rotational  period  of 

2  1 
243.16  days  (retrograde),  a  north  pole  direction  pf  RA&  =  270.3°  and  6^  =  66.7°,  and  an  epoch 

which  defines  the  longitude  of  the  subradar  point  as  —40°  at  0hUT,  20  June  1964.  It  is  also  re¬ 
quired  that  the  longitude  of  the  subradar  point  increase  in  the  positive  sense  with  time,  and  that 
northerly  latitudes  be  positive.  These  assumptions  are  based  upon:  (1)  the  growing  suspicion 
that  the  sidereal  rotational  period  of  Venus  is  exactly  equal  to  a  resonance  value  in  which  Venus 
makes  precisely  four  rotations,  as  seen  by  an  earth-based  observer,  between  successive  infe¬ 
rior  conjunctions;  (2)  the  direction  of  the  rotational  pole  of  Venus,  as  determined  from  the 
delay-Doppler  observations  of  Ref.  2;  and  (3)  the  epoch  defined  by  Carpenter.  The  last  assump¬ 
tion  is  a  convenience  to  facilitate  comparison  between  the  present  results  and  those  obtained 
earlier  at  JPL. 

Based  on  this  coordinate  system,  the  features  listed  in  Table  111  have  been  identified.  The 
ordering  is  according  to  decreasing  reliability:  the  first  four  features  listed  (1  to  IV)  are  con¬ 
sidered  to  be  quite  reliably  established;  the  last  four  features  (A  to  D)  are  less  firm.  Also 
listed  are  the  locations  of  matching  features  previously  identified  at  the  Jet  Propulsion  Laboratory 
during  the  1964  inferior  conjunction.  From  the  agreement  in  position  to  about  2°  in  longitude, 
and  from  the  fact  that  Venus  is  known  to  have  made  precisely  eight  rotations  between  1964  and 
the  present  observations,  the  rotational  period  of  assumption  (1),  above,  can  be  verified  as 
243.16  ±  0.Z  days,  retrograde.  The  systematic  agreement  in  position  between  five  of  the  Haystack 
features  and  those  identified  earlier  at  JPL  also  requires  that  these  features  be  firmly  attached 
to  the  rigid  surface  of  the  planet  and  that  they  have  a  lifetime  of  at  least  three  years.  The  failure 
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TABLE  III 

LOCATION  OF  RADAR  FEATURES  ON  VENUS 

Feature 

Locotion  fram 

Our  Observations 

Probable  JPL 
Identification 

Location  fram 

JPL  Observations 

Latitude 

(deg) 

Longitude 

(deg) 

Latitude 

(deg) 

Longitude 

(deg) 

Goldstein^ 

a 

-29  ±  2 

0  ±  ? 

Haystock  1 

-26  ±  3 

-1  ±2 

Corpenter^ 

F 

-26.7  ± 1.8 

0  ±  0. 7 

Haystack  II 

+24  ±2 

-81  ±  2 

Goldstein 

P 

+23  ±4 

-78  ±  6 

Hoystack  III 

+31  ±2 

-78  ±2 

None 

Haystack  IV 

-7  ±  3 

-65  ±  3 

Carpenter 

C 

-6.8  ±5.8 

-68.9  ± 1.3 

Haystack  A 

+23  ±  2 

-68  ±  2 

Carpenter 

D2 

+22.7  ±  1.7 

-70.0  ±0.7 

Haystack  B 

-12  ±  3 

-81  ±  3 

Carpenter 

B1 

—  11.9  ±  4.4 

-75.8  ±0.6 

Haystack  C 

— 13  ±  2 

-36  ±  3 

None 

Haystack  D 

+10  ±4 

-39  ±  4 

Nane 

to  observe  other  1964  JPL  features  in  some  cases  may  be  due  to  their  lying  outside  the  planetary 
region  observed  by  us,  and  in  others  possibly  by  the  confusion  associated  with  JPL' s  use  of  a 
single  parameter  (Doppler  frequency)  to  analyze  the  1964  observations. 

vn.  CONCLUSION 

Despite  the  fact  that  maps  of  the  Venusian  surface  presently  obtainable  by  radar  are  rela¬ 
tively  crude,  we  feel  that  the  techniques  and  basic  applicability  of  planetary  radar  interferometry 
have  been  adequately  demonstrated.  Because  of  the  opaque  cloud  cover,  radar  appears  to  be 
the  only  existing  means  of  mapping  the  planet  Venus.  We  hope  that  future  interferometric  studies 
can  be  made  with  a  more  complete  projected  baseline  coverage,  as  well  as  with  a  greater  system 
sensitivity.  At  the  present  time,  the  phase-stability  limitations  set  by  the  earth's  atmosphere 
have  not  been  reached,  and  the  corresponding  limit  to  the  resolution  of  the  technique  is  still  to 
be  realized.  Ultimately,  of  course,  a  limit  will  be  reached,  although  operation  at  longer  wave¬ 
lengths  or  the  use  of  an  extraterrestrial  phase  calibrator  may  extend  the  available  resolution 
even  further. 
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APPENDIX 


The  computer  program  listings  used  in  the  reduction  of  the  Venus  data  are  presented  in 
this  appendix. 

Real-Time  Programs:—  Figure  A-l  lists  the  eoded-pulse  real-time  program,  and 
Fig.  A-2  is  the  CW  real-time  program.  Both  these  programs  run  on  the  CDC  3300  computer 
at  Haystack.  Figure  A -3  consists  of  the  program  used  to  generate  an  ephemeris  tape  for  the 
real-time  fringe  rotation  of  the  CW  data. 

Coded  Post  Real-Time  Programs:—  The  post  real-time  phase  calibration  program 
used  to  refer  the  phases  of  the  coded  data  to  the  subradar  region  is  shown  in  Fig.  A -4.  Fig¬ 
ure  A-5  is  the  eoded-pulse  averaging  and  weighting  function  generation  program  which  averages 
the  power  for  a  complete  day  of  observing  and  generates  the  hemispheric  ambiguity  resolution 
functions.  The  programs  used  to  transform  the  data  to  the  planet's  coordinate  system  and 
display  the  data  on  a  CRT  are  given  in  Fig.  A -6. 

CW  Post  Real-Time  Programs:—  Figure  A -7  is  the  CW  phase  calibration  program 
which  refers  the  phases  to  the  subradar  region.  The  CW  averaging  program  (Fig.  A-8)  aver¬ 
ages  the  data  for  runs  with  the  same  projected  baseline  and  then  takes  the  transform  over  dif¬ 
ferent  projected  baselines.  Figure  A-9  lists  the  programs  used  to  display  the  CW  data  in  the 
planetary  coordinate  system. 
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Fig.  A-l.  Hayford  coded-pulse  real-time  processing  program. 
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Fig.  A-l.  Continued. 
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Fig.  A-l.  Continued. 
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Fig.  A-l.  Continued. 
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Fig.  A-l .  Continued. 
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Fig.  A-l.  Continued. 
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Fig.  A-l.  Continued. 
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Fig.  A-l.  Continued. 
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Fig.  A-l.  Continued. 
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Fig.  A-2.  Hayford  CW  real-time  processing  program. 
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Fig.  A-2.  Continued. 
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20  0  00  r 3 
*•1 

CONNECT  Chan  3 

1  R  3 

TEST  FOR  READY 

*♦2 

YES 

*-2 

5,3 

CLEAR  CHS  I  NT 

*-l 

5> 

CLEAR  CHAN  0  AND  2 

2,2 

connect  unit  2  to  channel 

*-l 

6,2 

SET  800  BPI 

♦-1 

16 1 1 

CH2 

C I  T ,  1 

12,1 

manual 

CIT.1 

SET  UP  CH2  INTERRUPT  IOC. 

17,1 

CH3 

C  I  T  ,  1 

SET  CH3  I  NT  ADDRESS 

3,3 

ENABLE/START  write  CORE 

*-l 

0 

REINIT 

ENABLE  INTERRUPTS 

4,3 

TEST  FOR  WRITE  CORE 

**2 

*.3 

2 , MO  I  SO , I T I TLE*17 


UJR  *-2 

EXS  2,2 

UJR  *-l 

UJR  * 

REINIT  SJ6  **3 

SEU  0,2 

UJP  *-l 

I OCL  R 

SJ3  **3 

CON  1.2 

UJP  *-l 

UJP  RESTART 


WAIT  TILL  FINI9HED 
WAP  FOR  INTERRUPT 


CONNECT  UNIT  1 


Fig.  A-2.  Continued. 
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PRQC  IS  THE  MAIN  PROCESSING  ROUTINE 


PROC 


A  « 1 


A. 2 


A, 3 


ENA  f  S 

0 

TAM 

220 

LDI 

ICNT,1 

ISO 

610,1 

UJP 

*♦2 

UJP 

TOQLQNG 

ENI 

2044,1 

ENI 

*11,2 

ldaq 

BUF0,1 

S9A 

BIAS 

STA 

ITR/2 

SHAO 

24 

SBA 

BIAS 

STA 

I T  1 , 2 

INI 

-4,1 

I  JD 

A, 1,2 

RTJ 

F0URIER2 

ENI 

2046,1 

ENI 

511,2 

ENI 

1022,3 

lda 

I  TR ,  2 

STA 

DBLOCK , 3 

MUA 

I  TR ,  2 

STA 

ZR  +  1 

STQ 

ZR 

lda 

ITI  ,2 

STA 

D8L0CK*1 , 3 

MUA 

ITI  >2 

SHAQ 

24 

ADAQ 

ZR 

SHAQ 

-5 

ADAQ 

spectra, 3 

STAQ 

SPECTRA, 3 

ldaq 

BuF  0 , 1 

SEA 

bias 

TAM 

700 

SHAQ 

24 

SBA 

BIAS 

TAM 

710 

MUA 

S  I  NO 

SHAQ 

1 

STQ 

tempi 

TMA 

7  08 

MUA 

COSO 

SHAQ 

25 

ADA 

TEMPI 

STA 

I  TR,  2 

TMA 

708 

MUA 

S  I  NO 

SHAQ 

1 

STQ 

TEMP2 

RESET  CLOCK 


THIS  LOOP  *  12  MS*  APPROX 
HAYSTACK  SPECTRUM  300  MS.  APPROX 


0*  A 

I  **2*R**2 


Q  p  A 


W*C*CQS0*W*S*S1N0*WC 
W  AC 


W , C*S I N  0 


Fig.  A-2.  Continued. 
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A,  4 


TMA 

71B 

MU  A 

COSO 

SHAG 

25 

SBA 

TEMP2 

sta 

ITI  ,2 

REM 

end  of  phase  c 

LDA 

COSO 

MUA 

COSDO 

SHAQ 

1 

STQ 

TEMP3 

LCA 

SINO 

MUA 

SINDO 

SHAO 

25 

ADA 

TEMP3 

STA 

COSO 

LDA 

SINO 

MUA 

CoSDO 

SHAQ 

1 

STQ 

TEMP4 

LDA 

COSO 

MUA 

SINDO 

SHAQ 

25 

ADA 

TEMP4 

sta 

SINO 

INI 

-4,1 

INI 

-2,3 

I  JD 

A. 2, 2 

RTJ 

/•'0URIER2 

ENJ 

1022,3 

ENI 

511.2 

LDA 

ITR.2 

STA 

DBLOCK*102«, 3 

MUA 

ITR.2 

STA 

ZR*1 

STQ 

ZR 

LDA 

ITI, 2 

STA 

DBLOCK-1029,3 

MUA 

ITI  ,2 

SHAQ 

24 

ADAQ 

ZR 

SHAQ 

-5 

ADAQ 

SpECTR A*1 024 , 3 

STAQ 

SPECTRA-1024,3 

INI 

-2,3 

I  JD 

A. 4, 2 

REM 

PICK  UP  NEXT  S 

LD  I 

ICNT.l 

LDA 

COSTBL.l 

STA 

coso 

LDA 

SINTBL, 1 

sta 

SINO 

LDA 

CoSDTBU  > 1 

STA 

COSDO 

LDA 

S I NDTBL » 1 

U*S*CQSO-W*C«SIN0*WS 
CORRECT  I ONi 8EQ I N  TRIO  RECURSION 


WESTFORD  SPECTRUM  277-369  MS 


DBL0CK*2047  INITIALLY 


PHASE  CORRECTIONS 


Fig.  A-2.  Continued. 
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RECORD 


ENDTEST 

END 


ST  A 

sindo 

IDI 

JCNTil 

TMA 

226 

ISO 

100  #  1 

STA 

I T ! ME1 # 1 

lda 

ICNT 

LOQ 

saved 

stag 

FIRST 

OUTWENT 

2,FIR5T,DHOCK*2Q48 

UJP 

*-  2 

ENA  ,  S 

0 

STA 

Saved 

EXS 

4,3 

UJP 

* 

I OCL 

10B 

EXS 

2,2 

UJP 

*-l 

set 

13,2 

UJP 

*•1 

SJ6 

*  +  3 

SSI 

V,2 

UJP 

♦-1 

EXS 

2,2 

UJP 

*-l 

I  OCL 

4 

CON 

1,2 

UJP 

*-l 

ENI 

21,1 

LDA 

C I T  S A VE ,  1 

STA 

C  I  T  ,  1 

I  JO 

*-2,1 

DINT 

0 

SCJM 

ioe 

SSJM 

7  B 

EINT 

0 

UJP*  I 

RUN 

YES, THEN  W*IT 

clear  interface 

WAIT  TILL  FINISHED  HRITINQ 
endfile  UNIT  2 


REWIND  AND  UNLOAD  IF  SSW6  ON 


CLEAR  CH2 

CONNECT  UNIT  1  TO  CH2 


RESTORE  C  1 T  TABLE 


Fig.  A-2.  Continued. 
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CH3  INTERFACE  INTERRUPT  PROCESSOR 


CH3 

UJP 

** 

COPY 

3 

COPY  STATUS  of 

CH  3 

TAM 

778 

ANA ,  S 

21008 

AZJ,NE 

fault 

EOBTEST 

TMA 

778 

AN  A ,  S 

608 

A  Z  J , EQ 

NOTEOI 

ENI 

8uF  0  ,  2 

ASE 

208 

ENI 

BuF 1  ,  2 

ST  | 

A, 1,2 

MODIFY  A.l  FOB 

BUFFER 

ST! 

A. 3, 2 

MODIFY  A  .  1  #  A  .  3 

for  buffer 

ENA,S 

1 

RAO 

ICNT 

I CNT ■  NO.  OF 

BUFFERS 

LDf 

I CNT  $  2 

TMA 

228 

ISO 

100,2 

STA 

I T I ME2 , 2 

SSL 

108,3 

UJP 

*•1 

EXITS 

ENI 

PROC , 1 

ST  I 

I NT1 , 1 

NOTfcOB 

ENA  ,  S 

Cm3 

ENI 

C  I  T  ,  1 

STA 

17,1 

UJP#  I 

Ch3 

RETURN 

r  AULT 

EN  A ,  S 

1 

RAD 

F  CNT 

UP  FONT 

SEU 

118,3 

clear  fault 

I  NT 

UJP 

♦-1 

UJP 

EqBTEST 

Fig.  A-2.  Continued. 
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CH2  TAPE  CHANNEL  INTERRUPT  PROCESSOR 


CH2 

UJP 

* 

ENAaS 

CH2 

EN  J 

C I  T  #  1 

STA 

16,1 

COPY 

2 

I  NCL 

0  0  04 

TAM 

76B 

IPA 

*02400 

azj,ne 

ERT 

UJP,  I 

CH2 

ERT 

ENA ,  S 

1 

RAD 

IERRTCNT 

UJP,  I 

CH2 

MANUAL 

UJP 

** 

SEL 

4,3 

UJP 

*-l 

ucs 

0 

SJl 

end 

IOCL 

10B 

UJP 

RESTART 

TOOLONG 

SBCD 

0 

UJP 

ENDTEST 

INTI 

EQU 

4 

BUFO 

EQU 

8192 

BUF1 

EOU 

12288 

BUS 

OCT 

177 

SAVED 

BSS 

1 

COSO 

BSS 

1 

SINO 

BSS 

1 

COSDO 

BSS 

1 

SINDO 

BSS 

1 

TEMPI 

BSS 

1 

TEMP2 

BSS 

1 

T6MP3 

BSS 

1 

TEMP4 

BSS 

1 

CITSAV6 

BSS 

25 

1H 

BSS 

DATA 

2 

I  tr 

BSS 

512 

I  T I 

BSS 

512 

SPECTRA 

BSS 

2048 

COSTBl 

BSS 

610 

SINTBL 

BSS 

610 

COSDTBL 

BSS 

610 

SINDTBL 

BSS 

COMMON 

610 

MOISO 

BSS 

6 

IRUN 

BSS 

1 

ititle 

BSS 

18 

IERRTCNT 

BSS 

1 

I  CNT 

BSS 

1 

FONT 

BSS 

1 

FIRST 

BSS 

2 

dblock 

BSS 

2048 

ITIM61 

BSS 

100 

ITIME2 

BSS 

100 

END 

MAIN 

RESTORE  C I T*16 
CLEAR  INTERUPT  CH2 
TEST  FOR  ERROR 
RETURN 

UP  TAPE  ERROR  COUNT 

MANUAL  I  NT  ROUTINE 
DISABLE  WRITE  CORE 

STOP 

SENSE  SWITCH  1  ON»  GO  TO  END 
off  clear  interface, 
restart 

set  bcd  fault  lite 


Fig.  A-2.  Continued. 
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SUBROUTINE  INPUT 

COMMON  MOlSO, JOURO. IHRO.FJMNO 

COMMON  I  RUN.  ITITLE(H).  IBR*  ICNT#  IFCNT,  I  BLOCK  (2050  ) 

COMMON  ITI ME 1(100). (TIME 2(100) 

TYPE  DOUBLE ( 2 )  SPECTRA 

COMMON/DATA/ I TR( 512). I T I ( 512 ), SPECTRA ( 1024 ) 

common/dat a/ i cospo( 610). isinpo (6io ) , icosdo(6io ) , isindo < 6io) 

DIMENSION  I  DENT ( 72 ) , I  DA  TES  C  72 ) 

RE AD ( 60 , 20 0  )  (ITI  TLE ( I ) » I *1* 18) 

200  F0RMAT(18A4) 

C  SENSE  SWITCH  3  On.  TEST  ROTATION 
GO  TO ( 1 , 13 )  SSWTCHF ( 3 ) 

13  WR I TE ( 5 V , 12 ) 

12  FORMAT ( 3lH  TYPE  MONTH,  DAY.  HOUR,  MINUTE) 
READ(58,11)M0IS1,j0UR1, IHRl,  IMN1 
11  FORMAT (412) 

READ(6,1000)IDENT 
READ(6»1000) IDATES 
1000  FORM A  T ( 72R1 ) 

I MN1  s  IMN1  *  2 
I  F  (  IMNlX.3,3 

4  IHRl  *  IHRl  *  1 

I MN1  *  IMNl  +  60 
I F ( I HR1 ) 5 , 3, 3 

5  JOUR1  *  JOUR1  -  1 
IHRl  *  IHRl  ♦  24 

3  READ(6,10  02)M0IS, JOUR, I  HR, IMN.R1.R2.D 
1002  F0RMAT(4l2,3El8,ll) 

IF(M0IS-M0IS1)3,6,1 

6  I F ( J0UR-J0UR1)3,7,1 

7  I F ( I  HR- I HR1 ) 3 «  8 ( 1 

8  IF(IMN-IMN1)3,1,1 
1  I RUNal 

RETURN 

END 


Fig.  A-2.  Continued. 
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SUBROUTINE  i n i t 

COMMON  MOlSO,jOURO, IHRO, FIMNO 

COMMON  IRUN, 1TITLE<18), IBR, 1CNT, IFCNT, I810CK(2050 ) 

COMMON  ITIME1U00)>ITIME2(100) 

TYPE  DOUBLE ( 2 )  SPECTRA 

COMMON/OATA/ITR(512) , I T I < 512 ), SPECTRA < 1024  ) 

COMMON/DATA/ I cOSPO (610 ) ,  ISINPO (610  > , ICOSD0  <610 ) , ISIND0(410 ) 
DIMENSION  I  HR (256), I W 1 (256), ITR2<256) 

DIMENSION  REPh‘15) ,  RPHD (15) , DOP (15) 

OO  TO  (20,21)  SSWTCHF ( 3 ) 

20  WR I TE ( 59 , 223 ) 

223  FORMAT ( 37h  TYPE  IN  ROTATION  IN  CPS.  FORMAT  F6.3) 

READ  <  98 , 222 )  FCPS 
222  FORMAT ( F6 , 3 ) 

ARGl'FCPS* 2, 0*3,1419926536 
'  ARG2*-ARG1/512.0 
DO  44  J*1 , 610 
AeARGl*( J-l) 

I COSPO ( j) *8388607. 0*COS( A) 

I SINPO ( J) *8380607. 0*8 IN( A) 

ICOSDO( J) *8380607. 0*COS( ARG2) 

44  I  S I NDO ( J)*83806O7.O*S!N(ARG2) 

00  TO  22 

21  WR I TE ( 59 , 112 ) 

112  FORMAT ( 55H  FORMAT  512  ,  TYPE  MONTH,  DAY,  HOUR,  MINUTE,  SECOND  ) 
READ (58, 110)  MOISO'JOURO, I HRO  *  MN  > ISEC 
110  F ORMAT (512) 

FlMN0«MN*ISEC/60.0 
TWOPI  *  2. *3. 1415926536 

INTERM  *  FIMNO  S  I F ( ( F I MNO  *  INTERM )  -  0.5)2,11,11 
2  I MNO  •  INTERM  S  GO  TO  13 
11  I MNO  *  INTERM  ♦  1 

13  I MN I  >  I MNO  -  2  S  I  HR  I  *  I HRO  S  JOURI  *  JOURO  S  MOISI  *  MOISO 
FIMNI  *  FIMNO  S  IF(  IMNI  )30 ,1,1 

30  I  HR  I  *  I  HR  I  -IS  IMNI  s  IMNI  *  60  S  |  F  <  I  HR  I  )3l,  1, 1 

31  JOURI  *  JOURI  -IS  I  HR  I  *  I  HR  I  ♦  24 

1  RE AD( 6, 10 02) MO  IS, JOUR, I  HR, I MN, REPH ( 1 ) , RpHD ( 1 ) , DOP < 1 ) 

1002  F0RMAT(4I2,3El8.H) 

IF(MOlS-MOISI >1,5,100 
5  IF< JOUR-JOURI >1,8,100 

8  IF( IHR-IHRI >1,9,100 

9  IF( IMN-IMNI >1.10.100 
100  DO  24  Ial,10 

24  BACKSPACE  6 
GO  TO  1 

10  DO  7  J*2» 15 

READ (6,  1002 )M0IS, JOUR, I  HR, I MN, REPH ( J ) , RPHD ( J ) , DOP ( J ) 

7  CONTINUE 

J  X  1  S  I  X  1 

15  LL  «  J  ♦  2 
AP  *  REPH  (  LL ) 

BP*(8.*(REPH(LL*1)-REPH(LL-1) )-<REPH<lL*2)-REPH<LL-2> ) ) / 12 . 
CP*(16.*(REPH(LL*1)*REPH(LL-1)-2.*AP)-(REPH(LL*2)*REPH(LL-2)-2.*AP 
1)  )/24. 

DPx( ( REPH (LL*2) -REPH (LL-2) ) -2 , * (REPH ( LL*1 ) -REPH < LL-1 > ) )/12. 

EPx( (REPH(LL*2>*REPH<LL-2)-2.*AP)-(l./3, ) * ( REPH ( LL*1 > *R0PH ( LL*1 > - 
12.*AP>  >/24. 

AD  =  RPHD ( LL  > 

BD= (8.* (RPHD ( LL*1 ) -RPHD (LL-1 ) > - ( RPHD < U*2 ) -RPHD (LL-2)  )  )/12. 
CD=(16.*<RPHD(LL*l>*RPMD(LL-l)-2.*AD>-(RPHD(LL-2)*RPHD(LL-2)-2.*AD 


Fig.  A-2.  Continued. 


64 


l>>/2*. 


DD*<  (RPHD(LL*2)-RPHDau-2) ) -2 , * < RPHD < LL*1 > -RPHD < LL-1 > )  >/12. 

ED*(  <RPHD(Ll*2>*RPHD(LL-2)-2.«AD)-U./3, ) * < RPHD ( LL*1 > *RBHD < LU - 1 ) - 
12.*AD>>/24, 

DT  s  FIMNI  •  IMNO  S  DELT  =  (1,  •  DOP < J*2 ) / . 78«E10 ) /60 . 

17  PH  I T I  *  ( ( <EP*DT  ♦  DP)*DT  *  CP)*DT  ♦  BP)*DT  *  AP 
ARG  *  PH l T I  *  TwOP l 

ICOSPOd)  *  COS(ARG)*8388607.  S  I  S  I NP  0  <  I  >  *  S  I  N  <  ARQ  )  *838860  7  . 
DPHT1  1  ( ( <  ED*DT  ♦  DD)*DT  ♦  CD>*DT  ♦  BD)*DT  ♦  AD 
ARG*-DPHTI*TWOPI/(60 .*512, ) 

ICOSDOU)  *  COS(ARG)*8388607  ,  $  ISINDOd)  <  S  I  N  (  ARG  )  *838860  7  . 

I  *  I  *  1  $  lF((FIMNI-FIHNO).Oe,lO,)l«#l9 
19  DT  «  DT  ♦  DELT  S  FIMNI  *  FIMNI  *  DELT 
I  F  < (FI MNI-IMNO) .LT. (0.9) >17.12 
12  J  *  J  ♦  1  S  IMNO  *  I MNO  ♦  1  S  GO  T3  15 

18  CONTINUE 

DO  6  1*1.20 
BACKSPACE  6 
6  CONTINUE 
22  NPTS=512 
IB!TS»12 

CALL  F0RINIT2( ITR, I T I , IWR.IWI.ITR2,NPTS.1.IBITS> 

CALL  CINIT 

RETURN 

END 


Fig.  A-2.  Continued. 
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SUBROUTINE  endrun 

COMMON  MOISO; JOURO. JHRO.riMNO 

COMMON  I  RUN. I  TITLE <18), J6R.1CNT, IFCNT, I  BLOCK (  2050  ) 

COMMON  ITIMEK100).  JTIME2<100) 

TYPE  DOUBLE ( 2 )  SPECTRA 

COMMON/DAT  A/ 1  TR<  512  ).m  (512),  SPECTRA  (1024) 

COMMON/DATA/ ICOSP0< 610), ISINP0(610),  ICOSD0( 610), I  SIND  0(610) 
DIMENSION  HAYSPEC(512),WESPEC(512) 

DIMENSION  P(200),X<200) 

EQUIVALENCE  (P,  I  BLOCK), ( X , I  BLOCK ( 50 0  ) ) 

c  convert  to  floating  point 
DO  2  1*1,256 

HA YSPEC ( I ) ‘SPECTRA ( 1*256) 

WESPEC< I ) sSPEcTRA  < I *768) 

H A YSPEC ( I*256)»SPECTRA( I ) 

2  WESPEC( I*256)=SPECTRA( 1*512) 

HR  I  TE  (  61, 20  0  )  (JTITLEd  ), 1*1,18) 

200  FOHMAT(IBAA) 

UR  I TE ( 61 , 20 0 0  )  IRUN 
2000  FORMAT ( 9X , 11H  RUN  NUMBER, 15) 

WR  I  TE ( 61 , 70  0  7  )  I CnT , I FCNT , I ER 

7007  FORMAT ( 6X.17H  No.  OF  BUFFERS  *,I8,16H  NO.  OF  FAULTS  *,I8.23H  NO.  0 
IF  PARITY  ERRORS  *,IB) 

DO  3  K*l,6« 

I FQ*  *256*8*  <  K-l ) 

J1*(K-1)*8*1 

J2*J1*7 

WRITE (61, 203)  < I FQ , ( HA YSPEC < J ) . J* J1 , J2 ) , ( WESPEC ( L ) . I* J1 , J2  )  ) 

203  FORMAT (IX, I  5 , 8F 16 . 0 , / , 6X . 8F16 . 0  ) 

3  CONTINUE 

WRITE (61, 101)  < ITIMEK  I  ),  ITIME2(  I  ),  1  *30,40  ) 

101  FORMAT ( IX, 2  I B ) 

C  TEMPORARY  TIME  PRINTOUT 

K*308 

DO  4  1*1,200,2 
P< I ) *H* YSpEC  <  K ) 

P< 1*1) *WESPEC ( K ) 

K*K*1 

X<  I  )*50*I/2 

4  X( I *1 ) »X < I )*.5 

CALL  SKETCH ( X , P ,  200  ) 

IRUN* IRUN*1 

RETURN 

END 


Fig.  A-2.  Continued. 
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HAYFORD  EPHEMERIS 
aFORTRAN  S/LS 

_ LL  -C_ -COMPUTATION  .Dfl  PARAMETERS  PERTINENT  JO  T-HE-JHAyF.ORp  1  -  INTERFEROMETER  _ 

experiment 

DIMENSION  MH( -1 ! 1 ) /M(-l ! 1 )/ APAZ( -1 ! 1 ) / APEL ( —  1 ! 1 ) / ALPHA ( *  1  1 1 ) /DELTA 

.  « 1 -Ull/REPHl  -11 11/I0£N-T1?2)  /  IDATES(72)/STtl2) _ 

NAME  LIST  FREQ/GCLA1/GCL01/GCDI1/GCLA21/GCL021/GCDI21/SMALLD/ 
*SMALAZ/OBJRRA/OBJRDE/O0JSP/OBJRAD/EPSILO/SMALLE/T/P/E/OMEGAI 
_T -j-283.l5_J_P_g_1013-»25  ;  E  ■»  10. _ 


1030 


1032 


ST ( 1 ) *23989*  709;ST(2)*3l322*982;ST (3) *37946 . 527; ST ( 4 ) *45279 *704 
ST (5) *52376. 353 ;ST( 6) *59709. 608; ST (7) *66806.344; ST( 8 ) *74139*627 
S.T  LSI  «81472_.863l-STllCl«.2l6i947iSTUU*9502.646;STU21  *16595*  356_ 
P I *3. 14 15926536; DET0RA.pl/18O.; FREQ-7840* ;GCLA1*42*62325*DET0RA 
C»  * 2997925E6  ;  GCD 1 1 *6368.4844944; GCLA21 *-0 *01049*DETORA 

GCL021*  -0 * 00672 *DETQRA;  GCl01»-71. 4&869*D£1QRA _ 

EPSILO«+0*3E-3;OBJRRA«90*3*DETORA; OB JR AD *608 9. 
SMALLD  •  1*239365568  ;  SMALAZ  *  201 *896389*DET0RA 

SMALLE-bI* 38138885*D£JQRAj  GCD 121 *SMALLD»S IN ( SMALLE ) _ 

0BJRDE*«66.7*DET0RA;0BJSP*2117664O*; OMEGA  I  * . 29670359E-6 
LICO  «  30 

READ17/103QJ  LAN  /MOISl/ JOURl/ IHR1/ IMN1/ IAN  / MO  I S2/ J0UR2/ IHR2/ IHN2 
FORMAT (1013) 

CALL  DAJU67 1  MOI SI / JOURl / JuLDl ) 

CALL  OAJU67!MOIS2j  J0UR2/-JULD21 _ 

WR ITE ( 5/ 1032 ) 

FORMAT(SlF  STANDARD 
»  INPUT  NAME  LIST*) 


PARAMETERS/  CR.  IF  NEW  VALUES/  TYPE  1/  CR  AND 


1  READ(5/ 1020) I 

1020  FORMAT ( I ) 
_ IF( I) 2/ 3/2 _ 

2  INPUT ( 5 ) 

3  CONTINUE 
WRI TE ( 4/ 1010 ) . 


1010 


FORMAT! 21H0EXPER I MENT 
READ ( 3/ 1000 ) IDENT 
READ! 3/ 1000 ) IDAT E£ 


'HAYFORD'  ) 


WRITE(6/1000> IDENT 
WRITE16/1000) IDATES 
FORMAT  L72RUL 


1000 

WRITE (4/ 1011) IDENT/  IDATES 

1011  FORMAT(40HOINTERFEROMETER  EPHEMERIS  DERIVED  FROM  .V72R1/72R1) 

WRJ  TE<4/1012)FREQ _ _ _ 

1012  FORMAT ( 32H0S ITE  1  !  HAYSTACK/ 


FREQUENCY  */F8*2/<  MHZ/  SITE  2  l  W 

•ESTFORD* ) 

WRITE<4/1013)AZ21/DETORA/EL21/DETORA/GCLA21/D£TORA/-GCL021/DETORA/ 


Fig.  A-3. 
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44!  *GCDI21 

_ 45.5 _ 1013  FORMAT (32H0DIFF.ERENCES  (SITE-2-#  SITE  1)  !/>  AZlMUlH  AXIS  DIFF.  jl 

46*.  #*/ F8 » 5/  *  DEGREES/  ELEVATION  AXIS  DIFF.  #*/F8.5/«  DEGREES*/*  GEOC 

47s  *ENTRIC  LATITUDE  DIFF.  »*/F8.5/*  DEGREES/  GFOCENTR I C  LONGITUDE  DIF 

49!  WRITE(4/1014)JULD1/IHR1/IMN1 

50!  1014  F0RMAT(30H00BJECT  !  VENUS/  JULIAN  DAY  «/I8/*/  GMT  ■$/2(lX/I2)) 

_ 51 !  _  WRITEIAz-1015) OB JRRA/DETORA/ OBJRDE/DETORA/ OBJSP/ OB JR AD  _ 

52!  1015  FORMAT! 38H0 COORD  I  NATES  OF  VENUS  ROTATION  AXIS  !  //*  R • A .  »*/F5.1/ 

53!  **  DEGREES/  DECL.  #*/F5.1/*  DEGREES/  SIDEREAL  PERIOD  «*/F8.0/*  SECO 

54  ! _ «NDPi  RADIUS  »»,F8.?.<  KM* ) _ 

55!  FACTOR  •  (77,6/T)*(P  ♦  4810. *E/T) 

56!  FREQ»FREQ*1 tE6; FC-FREQ/CJ  SLA1-SIN(GCLA1);CLA1»C0S(GCLA1) 

_ 571 _ DELTX-»--SMALLH»C0S(SMALLEUC0S(SMALAZ) _ _ 

58!  DELTY  »  SMALLD*C0S(SMALLE)*SIN(SMALAZ) 

59!  DELTZ  »  SMALLD*SI N ( SMALLE ) 

_ 601 _ IMN1*  IMM"1 .  i  IF  ( IMN1 19/19/19 _ _ _ _ _ 

61!  9  IMN1* IMN1+60  ;  lHRl.IHRl-1  j  IF( IHR1 ) 18> 19/ 19 

62!  18  IHR1#  IHR1*24  ;  J0UR1»J0UR1-1 

_ A3:  iq  T  MN?«  t  MN?> 1  ]  I F (  I MNP-6Q ) 1 0 . ?? . ?? _ 

64!  22  IMN2-IMN2-60  l  I HR2»  IHR2+ 1  ;  I F ( I HR2-24 ) 1 0/ 23/ 23 

65!  23  IHR2-IHR2-24  ;  J0UR2* J0UR2+1 

_ 66  I _ 10  _READ13a  1001 1  JDAYj  JEXPz  MHi  -1 1  /  M I  -  U  / 1  ST  z  ELOz  lELEXPz  AZD/I  AZEXPz  Ittl  / 

67!  »DOP/DTM/DDOP 

68!  IF(MH(-1).EQ. 24)00  TO  14 

_ 691  _ QO_TO_  15_ _ 

70!  14  READ! 3/ 1001 ) JDAY/ JEXP/ MH ( - 1 ) / M ( - 1 )/ I  ST/ ELD/ IELEXP/AZD/ IAZEXP/THRTM 

71!  */THRDOP/DTM/DDOP 

_ 72! _ 15  CONTINUE _ _ _ 

73!  1001  F0RMAT(3X/I7/2X/4I3/2(4X/I9/1X/I3)/F11.6/F12.3/2F11.4) 

74!  CALL  DA JU67 ( MO  I  SI / J0UR1 / JULD1 ) 

_ Z3J _ CALL  DAJU67_LMftLS2zJflUR2j.JuLD.2l _  _ _ 

76 :  IF! JDAY-JULDl J10/24/11 

77!  24  IF(MH(-l)*IHRl)10/25/ll 

_ 76) _ 25 _ LEAHl - 1 )  -  1MN1)  10z.lL/  11 _ 

79!  11  APAZ! "1 )»AZD*( 10***! IAZEXP-9) )*DETORA 

80!  CAZ  ■  COS! APAZ ( ■ 1 ) ) ;  SAZ  •  S IN ( APAZ ( -1 )  ) 

81! _ APEL(-l)  ■  ELD»10.«»(  IELEXP-9) _ _ _ _ _ _ 

82!  ROEL  •  APEL ( *1 ) *DETORA 

83!  CREL  •  COS(ROEL)  ;  SREL  -  SIN(ROEL) 

_ 841 _ TAU-"_180.»».)  «£;6*(  COS!  AP£l(  -1J  *Q£IQJ3A1/S  LNLAP£Li-ll*D£IORA,L_= _ 

85!  1  42 .5/ ( APEL ( »1 )  ♦  0 . 4 ) #*2,64 ) *FACTOR/PI  -  40* / ( APEL ( -1 )  +  2 • 7  )  *«4 • 

86!  APEL! »1 )  *  ( APEL  < - 1 )  +  TAu)*DETORA 

871 _ CEL  »  COS ( APEL (*11)  i  SEL  »  S I Ni APEL ( ■ 1 )  ) _ 


Fig.  A-3.  Continued. 
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8S5  CALL  JUDA67( JDAY/MOIS/ JOUR) 

89 ! _ SGMT*  ST ( MO  IS) t  ( JQUR*1  )  *236 >5 _ 

90t  SGMTU  »  SGHT  +  (MH(-1)*3600.  +  M< - 1 ) *60* ) * ( 1 •  +  236 . 5/ ( 24 . *3600 . ) ) 

91 :  STL*SGMTU'GCL01*( 12./PI )*360C. 

_ 921 _ SIND£L*SEL»SLAl*CEL*XAZ*CLAUCQSDEL«SQRTUjt.SINDEL»SINDEL) _ 

935  SNRDEL*SREL*SLA1+CREL*CAZ*CLA1)CSRDEL»SQRT( 1 t-SNRDEL*SNRDEL) 

94 S  SINTHE*»(CEL*SAZ)/C0SDELJC0STHE« (SEL*CLA1-CEL*CAZ*SLA1 j/cosdel 

_ 351  SNRTHE*»(CREL*SA2)/CSRDFL;CfiRTHF»(<;RFIl  »f.l  At.CRFI  *CA7*SI  A1  l/CSRDEl 

96:  THETA  *  ATAN2 ( SNRTHE/ CSRTHE ) 

97!  ALPHA(-l)  »  STL*PI/( 3600** 12 • )  -  THETA 

98: _ IF ( ALPHA  t  *1 ) »GE  » ( 2»  »P I ) ) Al  PHA( -1 )  >  AIPHA(-I)  -  2.*PI _ 

99:  DELTA ( -1 )  •  ATAN ( SNRDEL/CSRDEL ) 

loo:  DELZ  *  CEU*<SMALLD*COS(SMALLE)*COS( APAZ(-1 J-SHALAZ)  +  EPSILO)  ♦ 

101 : _ *SIN(  SMALLE)  *sel*smalld _ 

102:  DELZ  *  DELZ* ( 1 •  +  FACT0R*1 .E-6 ) 

1035  REPH(-l)  »  DELZ* ( FC  +  DOP/C) 

1041 _ _ 3EADl3.tlQQnjPAY^JEXPjMHlflnH10n-LS-IjELD/  IELEXP^AZD/  IAZEXP#  IHQ^  _ 

105:  *DOPO/ DTM/DDOP 

106:  IF ( MH ( 0 ) .EQ#24)G0  TO  20 

_ 1021 _ GO  TO.. 21 _ 

108:  20  READ(3/1001)JDAY/JEXP/MH(0)/M(0)/ IST/ELD/ IElEXP/AZD/ IAZEXP/ ThRTM/ 

109:  *THRDOP/DTM/DDOP 

lio: _ 21  CONTINUE _ 

111:  APAZ(O)  -AZD*(10.**( IAZEXP-9) }*DETORA 

112!  CAZO*  COS( APAZ ( 0 ) )  >  SAZO*  S I N ( APAZ ( 0 ) ) 

1138 _ APELtOl  ■  ELD*10.**( IELExP-9) _ 

114:  ROEL  *  APEL ( 0 ) *DETORA 

115:  CRELO  *  COS ( ROEL )  >  SRELO  *  SIN(ROEL) 

117:  1  42.5/CAPELC  0)  ♦  0 . 4 ) **2 . 64 ) *FACTOR/PI  -  40»/<APEL<  0)  +  2.7)**4. 

118:  APELt  0)  •  (APELt  0)  +  TAU)*DETORA 

119: _ CELO*  C0S(APEL(0)  )  l  SELO*  SIN(APEL(0)  ) _ 

120:  CALL  JUDA67I JDAY/MOIS/JOURO) 

1 2 1 :  SGMT  •  ST ( MO  IS )  +  ( JOURQ-1 ) *236 .5 

122! _ SGMTU  ■  SGMT  *  ( MH ( 0 )  *3600.  +  M(0)  *60.)»(1.  *  236. 5/(24. *3600/ ) ) 

123:  STL«SGMTU-GCL01*(12./PI )*360C. 

1241  SlNDEO«SELO*SLAl+CELO*CAZO*CLAliCOSDEO«SQRT( 1 . -SINDE0*S INDEO ) 

125: _ SNRDEO«SRELO»SLA1+-CRELO»CAZO*CLA1JCSRDFO*SQRT(1.-SNRDEO*SNRDEO) _ 

126:  SINTH0»-(CEL0*SAZ0)/C0SDE0 

127:  SNRTHO  »  -( CRELO*SAZO)/CSRDEC 

_ 128_: _ C0STH0«(SEL0*CLA1-CEL0»CAZ0»SLA1  )/COSDEO _ 

129:  CSRTHO  •  (SREL0*CLA1-CREL0*CAZ0*SLA1 J/CSRDEO 

130:  THETA  •  ATAN2CSNRTH0/CSRTH0) 

131 : _ ALPHA (0) _ »  STL*P I /  ( 3600 « *  12 « )  -  _IHETA _ _ 


Fig.  A-3.  Continued. 
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ALPHA ( 0 )  -  2 • *P I 

*C3S{APAZ(0>  -SYALAZ>  +  EPSILO) 


'AY>JEXPiMH(l 

s  pv  <»  p 


I  >!*(  1  ),  IST.»ELD>  IELEXPi  AZDi  I  AZEXPi  ThRTPU 


EPHEMFRIS 

IF  < ALpHA(01 

DELTA ( 0 )  «  ATAM ( SNRDEO/CSRDEC 
DELZ  =CELC*( S^ALLD*C5S(S^ALLE 
*SIM(S^ALLE) *SELC*SMALLD 
DELZ  =  DELZ* (  1  •  +  FACTOR* 1 • E -6 ) 

REPH(O)  =  DElZ*(EC  +  DOPO/C) 

12  READ(  3*  1001  )  JDAY,  JEXP^Mt  1  ) ,Y( 1  )/  I  ST i  ELD  <  I  ElEXP*  AZD/  I  AZExPj  TNiOOP, 
*DTMj  DDOP 

IE  (MH{1  )  »EQ»24)GP  TP  16 
00  TO  17 

16  READ { 3j 1001 ) J" 

♦  THRDOPj  DT P1  >  • 

17  CONTINUE 

APAZ1 1)  *AZP*<10.**( IAZEXP-3) )*DETORA 
CAZ  =  CCS( APAZ ( 1 ) )  ;  SAZ  *  SIN(APAZ(1)) 

APEL(l)  «  ELD*10.**( IELEXP-9) 

ROEL  *  APEL ( 1 ) *DETORA 

TAU  =>  180»*1*E-6*(C0S(APEL<  1 ) *DE TOR A  ) /S I N ( APEL (  1)*CETPRA) 

1  42 • 5/ ( APEL (  1)  +  0»4)**2»64)*FaCT0R/pI  -  40./MAPELC  1)  +  2« 

APEL (  1)  1  ( APEL (  1)  +  TAU)*CETORA 
CEL  «  CeS(APELtl)  )  ;  SEL  =  S I N ( APEL ( 1 )  ) 

CREL  *  COS ( ROEL )  j  SREL  =  SlN(ROEL) 

DELZ  =  CEL*(SMAlLD*C0S(SKALLE)*C0S(APAZ(1: 


>  7  )  **4 • 


•  SM.ALAZ  l  +  EPS  I LC  )  + 


*SIN(SNALLE)*SEL*SMALLD 
DELZ  *  DELZ* ( 1 •  +  EACT0R*1 »E-6 ) 

REPH(l)  *  DELZ* ( FC  +  DOP/C) 

REPHD  *  ( RE  PH  C 1 )  -  REPH ( - 1 ) ) /2 • 

SINDEL*SEL*SLA1+CEL*CAZ*CLA1 ;C0SDEL=S0RT( 1 . -SINPEL*SINDPL ) 
SNRDEL=SREL*SLA1+CREL*CAZ*CLA1; CSRDEL=SURT ( 1 • -S\RDEL*£NRDEL) 
SINTHE*-(CEL*SAZ )/CCSDEL;COSThE*( SEL*CLA1-CFL*CAZ*SLA1 J/COSDEL 
SNRThEk-(CREl*SAZ)/CSRDEL  j  CSRTHF*( S?EL*CLA1-CREL*CAZ*SLA1 j/csrdel 
CALL  JUDA67(UDAYjM0IS#JDuR) 

SGMT«ST(M0IS)+(o0UR-1 >*236*5 

SGMTU  =  SGMT  +  (  YH(  1  )  *3600.  +  1 )  *60O*(l*  +  236 . 5/ (  24  .  *3600 .  )  ) 

STLeSGMTU”GCL01* ( 12./PI )*36CC. 

THETA  -  ATAN2(S'RTHE>CSRThE) 

ALPHA(l)  =  STL.«PI/(3600**12.  )  -  THETA 
I F  (  ALPHA  (  1  )  *GF_.  {  2  •  *  P I  ))  ALPHA!  1 )  *  ALPHA  (  1 )  -  2**PI 

DELTA(l)  =  ATA'.fSNRDEL/CSRDEL) 

DX DP  =  DEL TX* (SIM  EC*SL A 1 *COSTHC  +  CL A 1 *COSDEC ) -DELT Y*S I NCEC *S I NTHO  + 

*  DELTZ*(CLA1*COSTHO*SINDEO-SLA1»COSDEO) 
DYDP«DELTX*SLA1*SINTHC+DELTY*C0STHC+DELTZ*ClA1*SINTH0 
GAMHa’»ATAN(  -DYDp/DXDP  ) 

ALPHAD  *  (ALPHA!  11  -  ALPHA  l  ■!).)/ lg.O_» 


Fig.  A-3.  Continued. 
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DELTAD  =  (DELTA! 1 )  -  DELTA ( -1 ))/ 120 • 

0MGXDP*CM£GA I  * (  COS ( OBJRRA- ALPHA ( 0 ) ) *C3S(23JRDE) *SNRDED-S IN ( DdJRDE) 
*«CSRDEO)+ALPHAD«CSRDEO 

OMGYDP* OMEGA  I* ( S IN ( OBJRRA- ALFHA ( 0 ) ) »COS(08JRDE) ) +  DELT  AD 
OMOPTH  =  SQRT<  DM3YDP*SMGYDP  +  3MGXCP*0M3XDP) 

CTOLD  *  2 • *  <  EC  +  D0PO/C)*OBJRAD*OMORTH 
GAMAXS=ATAN( -OMGYDP/OMGXDP) 

BL  =  (FC+C3PQ/C)*S2RT!DXCP*nXDP+Dy-DP#DYDP)*<-l.  +  FACTOR* 1 .£-61 
WRITE (6* 1002 ) MO  I S#  J5UR0#  Mh( 0 ) #  M ( 0 ) #  REPH ( 0 ) #  REPHD#  DOPG 
1002  F0RMAT(4l2#3E18.11 ) 

IF(LIC0-3O)3O#31#3O 
31  WRITE (4*  1052)JDAY#MCIS# jOuRC# I  AN 
1052  FORMAT ( 14H1JULIAN  DATE  =#I8#$#  GMT  DATE  =$313) 

WRITE!  4#1Q51 ) 

1051  FORMAT ( 106H0  GMT  APP.EL*  APP.4Z* 

1SELINE  BL  DIR.  C.L*DOP  COP. AX. DIR 

_ 2REES  _ DEGREES  WAVELENGTHS  wl/MIN 

3  DEGREES  SECONDS$//) 

LICO  =  0 
30  CONTINUE 

WRITE(4#1050)MH(0) #M(0># APEL ( 0 ) /DET3RA, APAZ ( 0 ) /DETORA# RFPH ( C ) * 

* REPHD #BL# GAMh W/DET ORA# CTOLD ^GAMAXS/DET ORA jT MO 
1050  FORMAT  (213#  IX#  2F10  •  3#  3X#  2F1C  .2#  F10 . 0#  3F10 . 1  #  4X#  F 10  •?.) 

LICO  =  LICO  +  1 

APAZ(O)  ;  APAZ(O)  =  APAZ(l) 

APEL(O)  ;  APEL(O)  =  APEL(l) 

ALPHA ( 0 )  ;  ALPHA(O)  =  ALPHA(l) 

DELTA(O)  ;  DELTA(O)  =  DELTA(l) 

REPH(O)  ;  REPH(O)  =  REPH(l) 

;  MH(O)  =  MH( 1) 

M  (  0  )  =  M  (  1  ) 


REL* PHASE  RFPhDOT  BA 
TI"E  DELAY/*  HR  M.\  DEG 

AVELE^GTHS  DEGREES  HZ 


APAZ! -1 ) 

APEL! -1) 
alpha ( -  1 ) 

DELTA(-1 ) 

.REPHC-l) 

MH(-l)  =  MH ! 0 ) 
M  (  —  1  )  =  M  (  0  )  J 


CAZO»CAZ;SAZO=SAZ; CFLC=CEl; SELO=SELlCRELO=CREL;SRELO=SRFL 
SINCEO=SINDEL;CCSDEC=COSDEL;SINTHO=3I\THE;C"'STHO=COSTHE 


SNRCEO=SNRDEL;CSRDEC=CSRDEL;SNRTHO= 
DOPO  «  COP  >  TM1  =  TMO  ;  TMQ  =  TM  j 
IF ( JULD2-JDAV) 13# 26# 12 

26  IF!  IHR2-MH! 1 ) ) 13#27# 12 

27  I F ( IMN2-M! 1 ) ) 13# 13# 12 
13  PAUSE  13 

IF ( SENSE  SWITCH  1 ) 7#  S 

8  STOP  _ 

END 

VARIABLES*** 


SNR  THE J  CSRTHP*CSRTHE 
J”' JRO  =  JP'JP 


Fig.  A-3.  Continued. 
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PROGRAM  HAYF0RD2 
CHARACTER  CH 

COMMON  RA(50) .DECS (50) >NDAYC50) .DISK  50) , I  DAT A ( 7940 ) fHA( 64,31), 
1CROTFR(64),W(2,64).CROT(64),8ROT<64> ,C (2,64. 31) , H (2, 64 ) rSAoTFR ( 64 > 
2,WF0UR3(64> 

DIMENSION  SCRaTA(64,31),SCRATBC2»64,31),WA(64»31  >,CH(84> 
EQUIVALENCE! I D*T A, SCRATA,SCRATI,CH >, (WA, HA) 

PI*3. 1415926536 
DO  2  1=1,50 

READ  1>  WRS, AMIN, SECS, DEGS, DM  I  NS , DSEC ,0 l ST < I) , NDA Y < l  ) 

1  FORMAT  (7F10.7.I3) 

RAC  I )=HRS*P J/12 ,0  *AM I N*P I /720 . 0 *SECS*P I /4320 0  .  0 

2  DECS ( I )*DEGS*Pl/180.0*SIGN( ( DM  I NS*P I /1 08 0 0 .0*DSEC«P J/6480  0  0 .0 
1),DEGS) 

3  CONTINUE 
REWIND  30 
I  NR  =  0 

IDA  =  0 
MO  I  S  =  0 
ICOUNT*0 

WESN=WESSsAWESS«AWESN»HAYN*HAYS*AHAYS*AHAYN»0 .0 
DO  4  1=1,64 
DO  4  K  =  1 »  31 

4  HA ( I , K ) =WA ( I «  K  >  =  C ( 1 » I ,K)*C(2, I , K  )  *  0 , 0 
BUFFER  IN(30,1>(IDATA(1),IDATA(21>) 

10  SO  TO  (10,11),UMTSTF(30) 

11  PRINT  12,  (CHC I >, 1*1,84) 

12  F0RMATC84R1) 

DO  19  1*1,72 
IFCIDA.EO.  0)111,19 

111  IFCCH(I)  .LE.  9)13,19 

13  IFCCH  (  I  )  ,  GE  •  0)14,19 

14  JFCCH  C  I  *1 )  ,  LE  .  9)15,17 

15  IFCCH  (1*1). GE.  0)16,17 

16  IDAY  =  CH  (p*10*CH  (1*1) 

GO  TO  18 

17  I  D  AY  =  CH  CD 

GO  TO  18 

18  I D A= I 

19  CONTINUE 

DO  130  I  *  I  DA , 72 

IFCCH  ( I ) . EQ . IRA ) 120 . 125 

120  IFCCH  C 1*1) .EC.1RUJ121.125 

121  MOis*e 

GO  TO  130 

125  IFCCH  (D.EQ.IRS)  126,130 

126  IFCCH  (I*1).EQ.1RE)127,130 

127  MO  I S*9 

130  CONTINUE 
IFCMOlS.EQ.  0)131,134 

131  WR I TE ( 59, 132 ) 

132  FORMAT ( 7HnO  DaTE) 

READ(58,133)M0IS, IDAY 

133  FORMAT ( 2  I  2  > 

134  T I ME= ( C  C C  CCH(80  )*0  )*64.*CN(8l) )*64.*CHC82) )*64.*CH(83) )*64.*CH(84 
1)  >/360CC000  .0 

NDAYTsMOIS*100*IDAY 
18  =  0 

DO  304  1=1,50 

I F ( NDAYT'NDAY ( I )  )  304,305,304 


Fig.  A-4.  Coded-pulse  phase  calibration  and  fringe  rotation  program. 
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309  1 8 *  I 
SO*  CONTINUE 

PRINT  306,NDAyT,NDAY( IB),RA< 1 0  ) , DECS U B ) , T I  ME 
30*  rORM*T(2I20,3F20 .5) 

CRA*RA< 1 8 ) *  C  R  A  < I C ) -BA ( ID) )*T I  ME/48.0 
CDEC*DECS( I  0 ) ♦ ( DECS ( I C ) -DECS ( ID) >  *T I M6/4  8 , 0 
ANO*CRA-179,7*PJ/180 .0 
NDAYN*(M0lS-8)*31  ♦  I  DAY  *212 

ST  I  ME* ( (TlME/24.0 ) ♦ ( NDA YN-212 > >*2.0*P»«!l. 002737909 
1  +(15.0/24.0*45,0/1440.0  *45 . 868/ ( 24 . 0*360 0  . 0  > ) *2 . 0*P I 
BN*32450 .0SBA*44.0*P 1/180 ,0*19.0 *P 1/(100.0*60.0) 

0HA«31.O*P1/18O,O  ♦  22. 0*PI/(180  .  0*60,0  ) 

BLN»0N*tCOS(CDeC>*SIN(0A>+COS(BA)*COS(BHA-STIME*CRA)*SIN(COSC) ) 
BLW'0N*(COS(0A)*S!N(BHA-STIME+CRA) ) 

BLA*ATAN(BLW/0LN) 

0LL*SORT(0LN*BlN*BLW*0LW) 

FLN*94.7*(COS(23.3*PI/180.0)*COS(CDEC)-SIN(CDEC)*SIN(23,3*PI/180.0 
1)*SIN(ANQ))*1.O60*(RA< IC)-RA( ID) )*180,0*30.0/PI-0.72*COS(CRA-ST 
2IMf )/D1ST( IB) 

FLU=94.7*SIN(23.3*PI/l«O.Oy*COS( ANO>- 1.068* (DECS( tC) -D8CS ( 10) ) 
1*180. 0*30. 0/P  I 
DOPA=ATAN(FLW/FLN) 

FI  =  SQRT(F|_N*F|_N*F|_W*FLw) 

ROT*2.O*PI*SlN<DOPA-BLA)*0LL*6.O55/(Ft*DlST( I B ) *14960 0 . 0 > 

THETA *360 .O*COS(DOPA-0LA)*BLL*6,O55/(DIST( JB) *149600 .0 ) 

DO  144  1*1,64 

SROTFR( I )*SJN( ( I-33)*R0T) 

144  CROTFR( I )*C0S( ( I-33)*R0T) 

PRINT  3O7,NDAyN,CRA,CDEC,ANG,STIME,BLA»0LL»DOPA«Fl,ROT»THSTA 
307  FORMATdX,  I3»l0F12.4) 

140  STIME*STIHE*15.*2.*PI/(60 ,*60 ,*24, ) 
PHLR*BN*CoS(8A>*'COS(CDiC)*itN(9HA-STlME*CRA>*2.0*P  1/86164,0 
PHL=BN*C0S(BA)*C0S(CD6C )*COS(BHA-STIME*CRA) 

PRINT  141 ,  PHLR »  PHL 

141  FORMAT(10X«2F20 ,6) 

142  DO  1*3  1*1,64 

SROT( I )*SIN(PHLR*2.0*PI*( (1-0/64. O  +  tCOUNT-l.O )*PHL*2,0*Pl  ) 

143  CROT ( I >«CO$<PhLR*2.0*PI*( ( 1-1 ) /64 , 0*1  COUNT-1 . 0 >*PHL«2 . 0 *P I ) 
IP*IPP*0 

20  BUFFER  I N ( 30 , 1 ) ( I D ATA ( 1 ) , I  DATA ( 3970  ) ) 

30  GO  TO(30,40,50,60),UNITSTF(30> 

60  PRINT  61, IP 

61  FORMATdX, 12HPARITY  ERROR, 110) 

IP* I P*1 

40  BUFFER  I N ( 3 0 , 1 ) ( ID*TA(3971), IDAT A ( 7940  ) ) 

41  GO  T0(4i, 70, 50,00). UNJTSTF<30> 

80  PRINT  81, IPP 

81  FORMATdX,  12HPARITY  ERROR, 110) 

IP*IP*1 

70  PRINT  62,  ID*TA(1),  IDATA(3971) 

82  FORMATdX, 6HREC0RD* 2120  ) 

I F ( I  P-1 ) 84 , 83 , 83 

83  DO  85  1*1,79*0 
85  I D A  T  A (  I  )  »  0 

84  DO  76  I R*l , 31 
DO  71  1*1,64 

K*4*IR+(I-l)*l24+lSIGN(l,-t*2*(I-32)) 
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H<1,I)*IDaTA<K>SH<2,I)*IDATA(K*1> 

W(l. I )  =  IDATA<K*2)*CR0T< I )*IDATA<K*3)*SR0T(  I ) 

71  W(2»I)=-IDATA(K*2)*SR0T(I)*IDATA(K*3)*CR0T(l) 

CALL  FOURIER3(H,WFOUR3»64,-1,0> 

CALL  FOURIER3(W,WFOUR3»64,-1,0> 

DO  72  1*1, 64 

L  =  I*ISIGN<32»1*2«<32-!  )  ) 

I  F  (  I R-3 ) 73  »  74 , 77 

73  WESN*WESN*(W<1,L>*W<1,L)*W(2,L)*W(2,L> )/(6A. 0*  2.0) 

hayN*HAYN»(H(1»L)*H(1,L)*H(2,L)*H(2,L))/<64,0*  2.0) 

GO  TO  77 

74  IF< 1-30)77,77,75 

75  I  F (  1  -36  )76.77,77 

76  WESS=WESS*(W<l.L)*W(l,L)*W(2,L)*W<2,L))/(5.0  ) 

WAYS*HAYS*(h(l»L)*H(l»L)*H(2,L)*H(2,L))/(5,0  ) 

77  MA< I » IR)*HA< I, IR)+H(1,L)*H(1»1>*H(2»L)*H<2»1) 
A=W(1,L)*CR0TFR(I)-W(2,L)*SR0TFR(I) 
B=W(1,L)*SR0TfR(I)*W(2,L)*CR0TFR(I) 

C  ( 1 , I , I R ) *C ( 1 , I , IR)*A*H(1,l)*0*H(2»l) 

72  C(2,I,IR)*C(2,I, IR)*B*H( 1,L)-A*H(2,L  ) 

78  CONTINUE 
ICOUNT* IC0UNT*1 

IF( I  COUNT >15 >142,20  0,20  0 

200  I  COUNT  *  0 

GO  TO  ( 908. 909  )  , SSWTCHF (5  ) 

908  M  =  4  $  GO  TO  9l0 

909  M*3 

910  CONTINUE 

A  A  =  B  B  *  0  *  0 
AWESS*AWESS*WESS 
awesn=awesn*wesn 
amays=ahays*hays 

AHAYNs AHAYN*HaYN 

WESN*WESS*HAYS*HAYN*0  .  0 
DO  201  1*31,35 
AA=AA*C(1, I ,M) 

201  BB=BB+C<  2, I , M ) 

PHASE* (180.0 /pl)*( AT AN(BB/AA»*(0.5*SIGN(0.5,AA>)*SIGN(PI.BB)) 

PRINT  202, phase 

202  FORMAK20X.9H  PHASE  *  ,F7.2) 

PRINT  203, AHAYN, AwESN,AHAYS, AHESS , C ( 1 , 33 . M ) , C < 2 , 33 , M ) 
1,HA(33.3),HA(33,2) 

203  FORMAT(lOX.8Fi5,2) 

*L  =  SORT  t AA*AA*BB*BB) 

DO  220  IR=l,3l 

DO  220  1*1,64 

AC*C(1« I* IR)*AA/AL*C(2»I , IR)*BB/AL 
AD=C(2,I,IR)*AA/AL-C(1»I,IR)*BB/AL 
ctl,  I. I R ) *  AC 
220  C ( 2 , I, I R ) *  AD 
DO  221  1*29,37 
A  A  s  C  ( 1 ,  I »  M  > 

BB*C  <2» I , H ) 

PHASE* <180 .0/Pl )*( ATAN(BB/AA)*<O.S-SIGN(0 ,5, AA) >  *S I QN  (  P 1 ,  BB) ) 
CC*SORT(C(l, I,M)*C<1»I»M)*C(2,I,M)*C<2, I , M ) ) /SORT <  A8S < (AH AYS- AHAYN 
1)*(AWESS-AWESN)) ) 

CC*CC*< INR*1) 

PRINT  230, PHASE, CC 
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23  0  FORMAT (10X.2HPH.F30 .  2 » 4HC0RR . F 3 0  .2) 

221  CONTINUE 

CALL  YYDISK(0,2,INR*15872*3968,WA,3968) 

CALL  Y Y W A  I  T  (  0  ) 

CALL  YYDISK(0,2, I NR*158 72*7936, C. 7936) 

CALL  YYWA  I  T ( 0  ) 

INR* INP*1 
no  222  I R  =  1 *  3i 
DO  222  1=1. 64 

222  H  A  ( I « I R ) *  W  A  < I , IR)=C(1< I, I R  )  *C ( 2 » I »  IR  > «  0 . 0 
GO  TO  140 

50  IAB=INR-1 

no  380  I R  =  1 »  3l 
DO  380  1=1.64 

380  MAC  I . I R ) =W A ( I , I R ) =C ( 1 » I » I R ) =C ( 2  »  I  .  IR)  =  0  ,0 
DO  40  0  I NR  =  0 , I A0 

CALL  Y  YD  ISM  0,1, I NR*158 72 *39 68 , SCR  A T A , 3968  ) 

CALL  YYWA  I  T ( 0  ) 

DO  391  IR=l,3l 
DO  39l  1=1.64 

391  WA(I»IR)*WA(I,Ir>*SCRATA(I,IP) 

CALL  Y  YD  ISM  0,1, I NR*1 5 872*7936. SCR ATB, 7936) 

CALL  YYWAITC0) 

DO  392  I Rs 1 , 3l 
DO  392  1=1,64 
DO  392  J=1 . 2 

392  C(J»I,IR)=C(J,I,IR)*SCRATB(J,I,IR) 

400  CONTINUE 

DO  601  K*1 »  3 
M=(K-1)*10*1 
N  =  M*1 0 

PRINT  602.  M.N 

602  FORMAT ClX,21HpOWER  AND  CORREL A T I  ON . 2  I  3 0  ) 

PRINT  603 

603  F0RMAT(1X,5HP0WER) 

DO  700  1=1.64 

700  PRINT  600, C  W A ( I , I R ) , I R «M , N ) 

PRINT  604 

604  F0RMATC1X.4HREAL) 

DO  701  1*1.64 

701  PRINT  600, (  C<1. I . IR) . IR=M,N) 

PRINT  605 

605  F0RMAT(1X,4HIM*G) 

DO  702  1=1.64 

702  PRINT  600, t  C < 2 . I » I R > . I R*M » N ) 

601  CONTINUE 

600  FORMATtlX,  HF12.0) 

GO  TOC900.901)  SSwTCHF(l) 

901  CALL  SEFFC20) 

BACKSPACE  20 
900  CONTINUE 

WRITE  (20)  NDAYN.CRA.CDEC.ANG, TIME, 0LA.BLL.DOPA.FL. ROT. THETA 

WRITE  (20)  WA.C.AHAYS.AHAYN.AWESS. AWESN 

END  FILE  20 

BACKSPACE  20 

PAUSE 

GO  TO  3 

FND 
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PROGRAM  HFP10T2 
— CHARACTER -MAPtMA - - 

COMMON  ICR<4,64,28) ,P<64,31)  ,CR< 2 , 64 , 31 ) »  INFO<2oO) 

DIMENSION  MAPU21.121>.MA(l01.101),  PP(  64,31) 

cOtJ  I VALENCF  <M*»CR)  I  IP, MAP)  ,  (  1CR»  PP) 

DIMENSION  IAA<l)$iAD«77777p  +  16  55B-L0CUAA<l))iTAA<  140*146000109 

P  I  =3 . 1415976536 

00  10  Ig1,64 

DO  10  IR  =  1» 28 

no  10  Ksi ,  4 

TCR<K, I » I R  >  ■  0 

10  CONTINUE 

11  N DAYA*ACRA-AC0EC*AD0PA*AFL=AThETA*0 . 0 
A AHAYN=0 ,0 

REWIND  30 
c*0  .  0 

20  PEAD<30)  n04Vn»CRA,CDEC,ANG»  TIME.  El  A » ELL » EOPA , FL» ROT .  THETA 
C,0  TO  <61,21>,  ECFCKF  <  30  ) 

23  RF A D  (30)  p,CR,AHAYS»AHAYN,AWESS. AW6SN 
FNCODl=<8M»60,  INFO)  ND  A  YN,  T  I  ME  •  F  L »  TMET  A 
60  FORMAT ( 4HD AY  ,  I7,1X»5MTIMR  , F5 , 2 » IX , 2 WFL » FS . 2 . IX . 6HTWF T A  ,F7.?) 
CAUL  L I  MI TS CO . 0 , 100 , 0 »0 . 0 , 120 . 0 ) 

CAI L  PD  I NTS <33 i0,AHAYS/<10.0**8.0)»1RH»1) 

CALL  POINTS(33.0,aWESS/(10.0**8.0),1RW,1) 

TAI  L  POINTS (1 ,0»  AHAYN/< 10 .0**8 . 0 ) *1RM»1> 

— CAtL-ROLNTSil. OtAWESN/ 110.0**8 ,0 ) » 1RW. 1 ) 

DO  500  I*1»6A 
**! 

PHASE*!  60-.37P  i  )*<  aTAN<CR<2»  1 ,4  )  /Cfi  <1.  !  .4)  >*<0 .5-SICMO ,3.CR(1.  I  »4) 
1) )*SIQN<P I ,CR<2. I .4) ) )*60 . 0 
C0RR*CR<1. I '4)*CR<1* I»4)*CR(2. I »4 )*CR<2. I .4) 

-CArfei.— P04NTSt  A#  PHASE.  1R*>1 ) 

CALL  POINTS(A,SQRT<CORR)/(lO,0**8.0),lRC»l) 

500  CALL  pOfNfS  <  A, P< I.4)/(10.Q**8.)»lRP»l) 

CAt  t  "L'ABEl5<4hFRF0  » lr3HP0W, 1 ) 

CALL  QRIDS<0. 0.10.0, 0.0. 20.0) 

CALL  GRAPHS <  I  NFO, 15, 200, 1 ) 

BAUSE 

no  TO  <20.504)  SSWTCHF ( 2 ) 

504  GO  T0(41,533).SSWTCHF<6) 

41  no  42  1*1,64 
no  42  IR*1»31 

42  pp<I.IR)spP(I,IR)+P<I»IR)-AMAYN 
AAwa¥N*AAhAYN*ARAYN 

«G  TO  43 

503  DO  40  IR=4»31 
CO  40  1*1,64 
IK= IR-3 

vs< 1-337*1. 008/FL 

-ftl^S8R44ri^-0“4-LS 4-. 94-1 R  )  /81 . 0  )  * *2  ) 

R2  =  SQRT(l,0-<  <83,94-IR)/81.0)**2) 

JF(R2*R2-Y*Y)  40,40,30 
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30  72  =  S1RT(R2*R2-Y*V  ) 

IF(R1*R1-Y*Y)  31,31*32 

31  71=0.0 
90  TO  33 

32  Zl=SQRT(Rl*Ri-Y*Y) 

33  7=(2.0*Zl+Z2>/3.0 
R  =  Ri 

<fCAs(0 .11**3  >/( (R+.ll*90RT(i, 0-R*R ) >**3) 

PHT=-THETA*Z*Pl/180 .0  :£=j= 

TCR(l.I.I<>3ICR(l,I»IK>*CR<l,I.IR>*1000.0*rOS<PWT>/(SCA*SQRT 
1 ( ARS (  (AHAYS"AwAYN)*(A1hESS*AHESN)  )  )  ) 
tcR(2.  I  1»IK>*CR(2»  ! . I R >*1000 . 0*S !N( PH l  >/<SCA«SQRT 

1  f  ARS  <  <  AHAY’?"  ARAyv)*  t  A  wESS-  AW  ESN) )  )  ) 
!CR(3,I.IK>BlCR(3»I»IK>+COS(Pt-I>*CO$(PHI>*lOOO.O 
TCR(4,  I,  U)*ICR(4,  I,  rK>*SIN(PHI  )  *S  I N  <  PH  I  >*l  000  .0 
40  CONTINUE 
43  0*0*1. 0 

-KDAYA=NDAYA*ND*YN 

acra=acra+cra 

acdec=acdec*cdEc 

ADOPA=ADOP'A*n0PA 
AFI.  =  AFI*Fl 
ATwFTA3ATHETA*TH5TA 
61  PAUSE 

00  TO  (900.20)  SSWTCHP ( 3 ) 

900  «0  TO  (90l*62)»SSWTCHF(6) 

901  no  902  tr*iTi<n 
no  902  N*1 » lOl 

902  M A ( L • N I  3 0 

50  903  1*1*64 
no  903  I»=1.3l 

«03  PP<  I»  IR>  =  PP(I, !R)/(C*1000.0  ) 

AAHAYNBAAHAYV/fC*lOOO ,0 ) 

HO  904  1=1*64 

904  PRINT  910,  (PP< I, IR), !R*1,14> 

PRINT  “911,  AAHAYN 

911  FORMAT  (10X.F20.fr) 
no  905  tsl «  64 

905  PRINT  9l9,(PPTl.fR>,fR*16r31> 

910  FORMAT ( IX , 16Ffr  •  ) 

fl  =  AF|_7C 

— nO-970-t* 1*1  di¬ 
ce  970  N*i » 101 
I=(L“51>*FL/60.0*33,5 

R  =  5ORT((<t-61>"(t-5n*(N-5l)*(N-51))/3600 .0) 

T  F  ( 1 . 0*R  >  9  70,970,974 

974  !R=(1,0-SQRT<1.0-R*R> >*81.0*3.94 

IF ( ABS( f-33,0)”ABS( 64.0 "R*l . 00  8/FL  >  >  973.973,970 
973  IF ( I '641  966.966.970 

966  IF(l-l)  970.967,967 

967  fFtIR  *3i>  968.968.970- 
96P  IF  (  IP  —  4  )  970.969.969 
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969  !SCA=(0.11**3>/<<R*.ll*SQRT<i.0-R*R>>**3) 
— «A(t.M>3PP( I « iR)/(SCA*PP(33»4> ) 

970  CONTINUE 
50  TO  971 

62  ”0  63  U ■ 1 i 1 0 1 
"0  63  N*1 « 1 0}. 

63  wA(L*N>=0 

50  603  K*l<4 
PRINT  601, K 

601  FORM*T<10X* 120  > 

50  602  1*1-64 

PRINT  600 , ( ICR (K. I . I R > . I R*l* 14 ) 

600  rORMAT<lX,14l8> 

602  CONTINUE 

nO  603  1*1 »  64 

PRINT  600, ( I  Or  Ik- I  * IR> « IR*15. 28) 


603  CONTINUE 


SNR*AMAYS*2.0/AwAYN 

no  64  1*1,64 

50"64  f-R*l » 20 - 

Al=ICR< 1. I . IR) 

A?* I  OR 02* I ,  I  R  ) 

A  3  =  I CR<  3 » I . IR) 

A4» I C  R  C 4. I » IR  ) 

TCP  (1* I* IR)*(A1*A4*A2*A3)*1000.0/(SNR*A3*a4) 

64  ICR  <2t3^IR>»<A1*a4-A2*A3)*10C0.0/(SNR*a3*a4) 
50  665  K*1 »  2 

PRINT  601,  K 
no  65  T*l,64 

PRINT  600, (ICR  <K,I,lR),IRil,14> 

65  CONTINUE 

50  665  3«t»64 

PRINT  600, (ICR  ( K ,  I, IR), !R*15, 28  > 

665  CONTINUE 
FL-AF^C 
no  70  L  =  1 , 1 0  1 

CO  70  N* 1 « 1 0 1 

)*rt/60 ,0*33,5 

RpSQRT ( ( (L'51) ♦<L-51)*(N-51)*(N-51  )  )/J600,0> 
IF( 1 . 0 -R) 70, 70,74 

74  T  R  =  ( 1 ,0-SQRT(l . 0 -R  «R ) )  *81 . 0*0 .94 

tF(A8SU*33,0)-A9S(64,0-R*1.0  0e/FL>  )  73,73,70 
73  IFU-64  )  66,66*  70 

66  IF ( I“1 )  70,67,67 

67  I F ( I R-28 >  68-68,70 
60  IF(IR-l)  70,69,69 
69  I F ( N'51 )  71.71,72 

72  ma(L,N>*ICR  (1*1.  IR)  /1000 
00  TO  70 
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71  MA(L’N>»ICR  <2» If IRJ/1000 
70  CONTINUE 
071  NDAY=NDAYa/C 

PECsACDEC*180,0/(PI*C> 

ANG  =  ACRA*180 . 0/(P  1*0 

THFTA*ATHETA/C 

FL=AFl'C 

DOPA=ADOPA*ieo.O/(PI*C)  < 

ENCODE (80 0.141, INFO  INCA Y, THETA. FL#C,SNR 
1  A!  FORMAT (4WD AY  ,17, 5Y»5HTHET A»5X,F7.2»5X»2HFL»F7.2»5X,1HC,5X,F7.2»5X 
1,3HSNR»5X,F7.2) 

CALL  LlMITSd. 0,101. 0,51. 0,111.0) 

*0  140  N*51 i 1 0 1 
CO  140  L*1 ' 10  1 

Al*N 

A2«L 

“QsMACL, N) 

140  CALL  POINTS! A2. Al, MQ, 1 ) 

CALL  GRAPHS  ( INFO, 15, 200,1) 

CALL  LIMITS  (1.0,101.0. >9, 0,51.0) 

CO  150  N*1 .51 
CO  150  L=1'101 
A1  =  N 
A2  =  L 

mQ=MA(L«N) 

150  CAI L  POINTS! *2,A1.M0,1) 

.CALL  GRAPHS  ( 0,0. 0,1) 

155  r>0  156  K  =  1 » 12l 
ro  156  M31 , 1 Zl 

156  wap(K,M)*0 
ANG  =  ANG-179 , 7 

CDFC  =  COS(D?C*P  1/180,0) 

SDFC=$IN(DSC»PI/180,0) 

CAnG=COS ( ANG*Pl/180.0) 

SAKG*SIN(ANG*Pl/180,0) 

CTH»COS(23.3*PI/180,0) 

STHsSIN(23.3*pI/180.0) 

ALATR=(180.0/pI  >*ASINF(-CTH*SCEC-ST).-*SANG*CDFC) 
ALCNR=(180.0/pI)*ATAN<(SANG*CCEC*CTH'SrEC*STH)/CANG*CDEO)- 
ALCNRsALONR*l 6. 0*(NDAY-241)*36 0.0/245.0 

pla*atani<cang*sth/  (Ctm*cdec-sth*scec*sang)  ) 

APA=PL  A*l80 ,0/PI 
CDCPsCOS(DOPA*P 1/180.0) 

SDCPSSIN(DOPA*PI/180.0) 

nrun*nday*iooqo»time*ioo 

print  160,NRUn»NDAY,DEC» ANG.TI-ETAiFL.COPA, ALATR, ALCNfi,  APA 
160  FOPMAT(2UO,8plO,2) 

PE=PI/180,0*(16,0  +  (NBAY -^2  41  7*360. 0/245.0  ) 

CO  180  L* 1 , 1 0 1 
CO  180  N* 1 . 1 0 1 

Y»(L-517*CBOR/60 .0-<N-517*SEOF/6OfO 
Z*(U-51)*SDOP/60 . 0*<N*51)*CL OP/6  0  »  0 
!F(l-Y*Y-Z*Z>l80.l61,l61 
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161  *sSQRT<l-Y*Y+Z*ZT 

yA=X*CANG*CCEC-Y*SANG*Z*SCEC*CANG 

YA*X*(CTW*SANg*CDEC-STH  +  SCEC  >  +  Y*CANG*eTM*Z*(CTH*SCEC  +  SANG*STH  +  rDFC 

%AsX+<'CTH*SDEC-«TH*SANG*COEC>-Y*CANG*STH*Z*CCTH*CCEe-S?M«SDEC*SAN 

1G)- 

X8=XA*C0S<PE)*Ya*SIN(PE) 

Y8=XA*SIN(PE)+YA*C0S(PE) 

K  *  YB  +  6  0  .0  +  61.5 
*sZA*6fr.O  +61.5 
t  F ( X8 )  170*175  *175 
170  MAP(K»h)*l8S- 
GO  TO  180 

175  MP(K,M)=MAP<K»M>+  MA(L.N) 

180  CONTINUE 

~M  190  J»1^19 
AA*( J-10>*PI/18,0 
CO  190  1*1 »36l 
PB*U-181)*PI/'560.0 
K»COS<AA)*SIN(BB >*60.0  +  61. 5 
P*SIN(AA>*60. 0+61.5 
If  (MAP<K*mU-190*181»1»0 

181  w AP ( K»  M ) *1R* 

190  CONTINUE 

CO  200  J*1 *  19 
AA=(j-io>*pi/ie.o 
CO  200  I  E 1 »  3 6 1 
PB=<i»181>*P 1/360,0 
K*CCS(BB>*SIMAA)*60,0  +  61.5 
m*S I N ( BB > *60 , 0  +61.5 
IF(MAp(K»M>)  200,191*200 

191  *AP(K,M)slR* 

200  CONTINUE 

CAi.  L~ ~L  I MITSI1 ,0  ,61.0*61 . 0* 1*1 «  Q ) 

CO  210  K* 1 » ^1 
f0_210  M=fclil21  - 
Y APKBMAP ( K  »H  > 


AK  =  K 
A  M  =  M 

lftMAfK )209i 210,209 

209  CALL  POINTS* AK» AM.MAPK* 1 > 

210  CONTINUE 

CALL  TJPAPHS  40r6-*6*l) 

CALL  LIMITS! 61 . 0,121. 0,61, 0*121.0) 

CO  220  K*61 » 121 

CO  226  N»61?421 

AK*K 

AMBM 

VAPKbMAP<K*M> 

IP(MAPK)2l9,220,2l9 
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219  CAIL  POIMS(Ak»  AM,MAP«»1) 

220  CONTINUE 

CAI l  GRAPHS  (0/0, 0.1) 

CALL  LIMITS  (1,0  .61.0,1.0.61.0) 

TO  230  K *  1 ,61 
TO  230  M*1  61 
AK*K 
AH  =  M 

yapk*map ( K  >  M ) 

IF(MAPK)229,230,229 

229  CAM  PCINTS(Ak,AM.MAPK,1> 

230  CONTINUE 

CAIL  GRAPHS  (  0  ,  0 , 0  > 1 > 

CALL  LlMITS(6l. 0,121, 0.1. 0,61.0) 

TO  240  K*61 , 121 
CO  240  M* 1 >  61 
AK*K 
ah  =  M 

„APK  =  M AP ( K  «  M  > 

IF(MAPK)239,240,239 

239  CAIL  POINTS(Ak,AM.MAPK,1> 

240  CONTINUE 

CAIL  GRAPHS  ( 0 , 0  ,  0 , 1 ) 

PAUSE 

CO  T0(242.241)  SSWTCHF (1 > 

241  CALL  SEFF  <  20  ) 

BACKSPACE  20 
CALL  SEFF  (10) 

BACKSPACE  10 

242  CONTINUE 
COPAA*COPA-APa 

GO  TO(800,e01),SSWTCHF(6) 

801  WRITE(20>  NDAY.THETA.FL.CCPAA, ALATR,ALCNR,nOPA,APA,C,DEO.A\G 
GO  TO  802 

800  WR  T  TE( 2  0 )  pP.aAhAYN 

ro  806  m.6A 
ro  806  19=1,28 

CO  806  K*1,4 

806  ICR(K,  I  t  I  R  )  «  0 
GO  TO  803 

802  WRITE (20) ( ( ( ICR( I « J»K> , 1*1,2) , J*1.64) ,K*1,28) 

WRITE  (10)  MAp.NDAY 

DO  807  1*1. 6« 

CO  807  IR=1.3l 

807  PP ( I , I R  )  *0  .  0 

803  PND  FILE  20- 
PACKSPACE  20 

end  file  10 

BACKSPACE  10 

GO  TO  (  11  .  250  )  SSWTCHF  ( 4  ) 

250  FND 


Fig.  A-5.  Continued. 
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PROGRAM  VENRLOl 

C  LU  24  IS  THE  INPUT  TAPE.  LU  27  IS  THE  OUTPUT  TAPc. 

C  SSW4  DOWN  TO  OMIT  PHOTOS. 

C  SSW5  DOWN  TO  CANCEL  A  SETTING  COMMAND. 

C  ARRAYS  IN  COMMON 

COMMON  IDISPLAYCIOO ) .PWRARAY (64,31) , HTFCN(2,64,28) , I NP (101*101) 
COMMON  bCATLAW<31) 

C  VARIABLES  IN  COMMON 

COMMON  EL  PRlMt,BPRlME,FCL,D0PA,N0ISEBX,A,R8SHiFT,ISUR,ISCALE>lNP 

1  MIN.  INPMAX, iPOlNT.ISTEtR. RADIAN. I D A Y NUM , D A YR 'JnS . I FL*G . P I . SMDYrN . 

2  INPMAX1.  I NP M[N1.XMAC,AAHAYN, DEC, ANG. DOPPLER  A , PL ANET A  ,  FBSH  I  F T 
DIMENSION  INPl(lOl.lOl) 

EQUIVALENCE  (  INP1,PWRARAY), (DOPA, GAMMA) 

REWIND  27 
I  FI  AG  =  3765432lE 

P I =  3 • 14 1592654  I  RADIAN  =  PI/180.0  $  C -2 . 997929E+ 0 8 

XMAC=3777777b  %  DT=5.0E-04  $  I NPH I N=377777778  I  INPMAX=-1 

1  WRITE(59,99) 

99  F0RMATH7H  LOAD  DATA  CARDS.  /10H  PRESS  GO.  ) 

GO  TO( 51 »  213 ) » SSWTCHF ( 4 ) 

213  PAUSE  4444 

51  I  ST EER  =  2  4  GO  TO  52 

52  CONTINUE 

PO  53  KRAnGE*1 *  31 
TO  53  NFREQI1,64 

53  PWRARA Y ( NFREO , KRANGE  )  =0  .  0 

PEAD(60,98)ISu6,ISCALt,IPOINT,NOISEBX,RBSHlFT,FCLX,A,ID.FBSHlFT 
9fl  FORMAT ( 4  I  5, 3T 10 . 0 , I  5, F10  .  0  ) 

I F ( ISUB)215, 214,215 

214  I SUB=3 

215  I F ( ISCALEJ217.216.217 

216  ISCALE'100 

217  IF ( I  POINT >219. 218. 219 

218  I PO I  NT  =  2 

219  IF (NOISEBX)22l, 220.221 

220  NO  I SEBX  =  2 
R8SH I  FT  =0 .975 

221  IF ( A ) 223 , 222  »  223 

222  A  =  6 . 0 55 

223  a»A*1.0E*06 

WRITE (61, 90) I  SUB, I  SCALE* I POl NT ,N01 SEBX » RBSH I  FT , FCLX ,  A,  ID.FBShIfT 
90  FORMAT ( 4I5,3FiO . 4, Is.FlO .4 ) 

19  READ< 24  )  I DAYNUM, THETA, FCL , DOPA ,B  PRIME, EL  PR  I  ME , DOPPLER  A, 

1  PLANET  A.DAYRUNS,  DEC, ANG 

GO  T0(198,199)E0FCKF(24) 

199  READ  (  24 )  (((  I WTFCN( I , J,K) , I =1,2) , Jsl,64) ,K  =  1.28) 
PEAD<24)PwRARaY,AAHAYN 
IF ( SSwT  CHF ( 5 )  ,LT.  2)21,22 
22  I F ( I D  .EO.  0  .OR.  ID  .EQ.  I D A YNUM ) 2 0 . 1 9 


Fig.  A-6.  Coded-pulse  display  and  coordinate  transformation  program. 
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20  IF(FCLX)225. 22^,225 
225  FCL=FCLX 

224  WRITE(6l, 97HDAYNUM, THETA. FCUDOPA, EL  PRIME.  B  PRIME. DAY  RUNS, 

1  AAHaYN.DOPPlERA, planet  a.dec.ang 
97  FORMAT(12hODAy  NO.  IS  16,  2X  16H  FRINGE  RADIUS  IS  EU.J.2X 

1  35H  CENTER-T0-LIM8  DOPPLER  IN  HERZ  IS  F9.4/  18H  DOPPLER  ANGLE  IS  F9.4. 

2  F9.4,  2X  29H  SUB  RADAR  LONG.  AND  LAT.  IS  F8.4,  2X  4H  AND  2X 
3F8.4/16H  NO.  OP  RUNS  IS  F6.0/22H  DAYS  AVERAGE  NOISE  IS  EH.  3/ 

451H  ANGLE  BETWEEN  DOPPLER  AXIS  AND  CELESTIAL  NORTH  IS  F9.4/ 

554H  ANGLE  BETWEEN  HESPERIAN  NORTH  AND  CELESTIAL  NORTH  IS  F9.4/ 

6  41H  MEAN  DECLINATION  AND  RIGHT  ASCENSION  IS  F9.4.2X  4H  AND  F9.4) 
KDEX=NOISEBX*l 

no  200  krange=kdex»3i 
DO  200  NFrEO  =  1 »  64 

PWRARaY(NFREQ,KRAnGE)=AMAX1(0 . 0.PWRARAY(NFREQ.KRANQE>  ) 

200  CONTINUE 

C  APPLY  SCATTERING  lAW,  ETC,  ,  AND  XFORMAT I  ON  To  GRID  OF  MAPPING  PLANE. 

CALL  MAPXFORM 
GO  T0(23.24),SSWTCHF(4) 

C  MAKE  TEST  PHOTOGRAPHS  . 

24  CALL  TESTFOTO(l) 

23  CONTINUE 

c  at  a  later  date  enter  editing  routines. 

GO  TO(208,212)  .SSWTCHF (4) 

212  WR I TE ( 59. 101 ) 

101  F ORMAT ( 35h  TYrE  1.  IF  DATA  IS  UNSATISFACTORY.) 

READ(58,102)OUESTION 

102  FORMAT(FIO.O) 

IF (QUESTION)208, 208,209 

209  DO  109  KK*1 , 3 
BACKSPACE  27 

109  BACKSPACE  24 
GO  TO  1 

208  GO  TO ( 210 , 211 ) I  STEER 

210  CALL  TAPEADD 
GO  TO  1 

211  CALL  TAPEPACK 
GO  TO  1 

198  WR I TE ( 59. 103 ) 

WR  I  TE ( 61, 103  ) 

103  FORMA  r  t lBH  EOF  ON  DATA  TAPE./ 

1  48H  LOAD  NEW  DaTa  TAPE  OR  ABORT  RUN  AND  SAVE  LU  27.  ) 

21  REWIND  24  $  GO  TO  1 

END 


Fig.  A-6.  Continued. 
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SUBROUTINE  TEsTFOTO(K) 

C  MAKES  PHOTOGRAPHS  and  PROVIDES  SAMPLE  PRINTOUTS  TO  TEST  GOODNESS  OF  RUN 
C  IN  A  ROUGH  FASHION 
C  ARRAYS  IN  COMMON 

COMMON  I  DISPLAY! 100 ) ,PWRARAY ( 6  4,31 ) . IHTFCN(2,  64  ,  28) , INPCIOl# 101 ) 
COMMON  SCATLAWOl) 

C  VARIABLES  IN  common 

COMMON  6L  PRImE,BPRIM£,FCL,DOPA,NOISEBX,A,RBSHiFT,ISUB.ISCALE<INP 

1  MIN. INPMAX, I  POINT.  I  STEER. RADI  AN. I D A YNUM . D A YRUNS .  I FL AG . P I , SMDYRN . 

2  INPMAXl,  INPmIN1,XMAC,AAHAYN,DEC,ANG. DOPPLER  A  ,  PL ANE T A , F8 SH I F T 
DIMENSION  I  NP 1(101.101) 

EQUIVALENCE  (INP1,PwRARAY),(DOPA. GAMMA) 

BP  =  BPRIME/RADIAN  $  ELP=EL  PR  I ME/R  AD  I  AN 
DO  TO (401, 402), K 

401  HR  I  TE( 61. 10  0  )  IDAYNUM.FCL. ELP. BP, DA YRUNS, AAHAYN, DEC, ANQ, DCPPLERA 

1  .PLANETA 

100  FORMA T ( 13H1DAY  NUMBER1  16, 2X  24H  CENTER-TO-L I  MB  DOPPLER5  F8.3» 

1 2 X  18H  SUBRADAR  LONG  IS  F8.3,  2X  12H  aND  LAT  IS  F8.3/ 

2  27H  NO.  of  RUNS  IN  THE  DAY  IS  F6.0,  2X 

3  46H  THt  AVERAGE  NOISE  (FIRST  TWO  RANGE  BOXES)  IS  E10.2/ 

4  22H  THE  MEAN  DEC  AND  RA  5  F9.4,  1H,  2X  F9.4,  1H./ 

5  55H  OOPPLtR  AND  HESPERIAN  NORTH  W.R.T.  CELESTIAL  NORTH  IS  F9.4, 

6  2x  4 H  AND  f 9 . 4  ,  2X  1H  .  ) 

C  DO  34  MM® 1 » 1 0 1 

C  M=102-MM 

C  wRI TE ( 6 1 . 1 0 1 ) (  INP(M.L) ,L  =  46,55) 

C101  FORMAT ( 2X  1 0  1 1 0  ) 

C34  continue 
C  IMP  =  3  77  777  7&B 

IMP®2750000 
GO  TO  403 

402  WR I Tfc( 61 . 102 ) SMDYRN , I UAYNUM 

102  F0RMAT<19h1T0TAl  NO.  OF  RUN3=F10.0,  2X 

i  33h  the  last  included  day  number  is  16) 

DO  454  MMrl.lOl 
M=102-MM 

WRITE(61.101)  (INP1(M,L)»Ls46,55) 

C454  CONTINUE 

I MP  =  I  NPMAx/2 

403  WRITE(59,103) 

103  F0RMAT124H  READY  POLAROID.  HIT  GO.) 

PAUSE  1717 

73  CALL  RESETD ( *6 . -10 > 

UL  =  V  6  ®  0  iUR  =  VT®100.0 
CALL  SCALE< UL, UR. VB.VT. 0.23, 0,23) 

DO  26KK=l,ld 

26  CALL  GR I D(  0  . -60  .  0 ,40 . 0  . -40  .  0 . 60  .  0  .  10  .  0  .  10  .  0 . 2. 2  , -2  . -2. 5. 5, l ) 

413  GO  TO( 407, 40e  )  ,  K 

40  7  CALL  INTRPLOT ( INP.101, 101. 0,2,0 ,2. 5, 5, 0 . 45, INPMIN, IMP  ) 

GO  TO  409 

408  CALL  INTRPLOT ( INP1, 101, 101, 0,2, 0,2, 5, 5, 0, 45. 0, IMP) 

40  9  CALL  ADVF i LM ( 0 . 1 ) 

WRI  TE<59. 103)  IMP 

105  FORMAT (2lH  The  PRESENT  PMAX  IS  18/ 

1  33H  TYPE  IN  A  NEW  PMAX  IN  FORMAT  17.  ) 

READ (38,106)  ImP 

106  FORMA  T (  17) 

IF ( IMP ) 71, 71 . 73 
71  RETURN 
END 


Fig.  A-6.  Continued. 
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SUBROUTINE  TApEADD 
C  ARRAYS  IN  COMMON 

COMMON  IDlSPLAY(100>.PWRARAYl64,31>,IwTFCN(2,64,28),INPtl0l,l0l) 
COMMON  SCaTLAW (31  ) 

C  VARIABLES  IN  COMMON 

COMMON  EL  PRIME.BPRIME.FCL.COPA.NOISEBX, A.RBSHIFT, ISUB. ISCAlE. INP 

1  MIN. INPMaX, I  POINT. I  STEER. RAD  I  AN. IDA YNUM. DAYRUNS. I  FLAG. PI , SMDYRN, 

2  INPMAXl,  INPMINJ.XMAC. AAHAYN, DEC. ANG. DOPPLER  A , PL ANE T A . F BSH  I  F T 

DIMENSION  INPKlOl.lOl) 

EQUIVALENCE  ( INPl.PWRARAY)  ,  (DOPA.GAMMA) 

GO  TO  15 
FNTRY  TAPEPACK 
15  CONTINUE 

nOPA  =  DOPA/RADI AN  $  EL  PR  I  ME  SEL  PR  I  ME/HAD  I  AN  SBPR I  ME  =  BPR I  ME/RAD  I  AN 
WR I T  E ( 27 ) IDAYNUM.DOPA.EL  PRIME, B  PR  I  ME . FCL . NO  I SEB  X , RB SH I F T . I NPM I N . 
1  . I NPMAX , dayruns 
WRITE (27) I NP 
WR I TE ( 2 7  > IFLAG 
I STEER=2 
RETURN 
FND 


SUBROUTINE  TEST  DATA(KRANGE,NFrEQ,X1,X150.X2,X2SQ,X1PRSU.X2PRSQ, 

1  X1PR . AREA1 , AR6a2, XMULT  ) 

COMMON  lDtSPLAY(lO0),PwflARAY(64,31>,IBTFCN(2.64,28>,INP(lOl»lOi> 

COMMON  SCATLAg(31) 

c  variables  in  common 

common  EL  PHImE.BPRIME.FCL.DOPA.NOISEBX, A.RBSHiFT.  I  sub. I  scale,  inp 

1  MIN. INPMaX, IPOINT. ISTEER.RADI AN. JDaYNUM, DAYRUNS, I  FLAG. PI. SMDYrn. 

2  I  NPM  A  XI.  INPMIM.XMAC,  AAHAYN,  DEC,  ANG,  DOPPLER  A  ,  PL  ANE  T  A  .  FBSW  1  F  T 
DIMENSION  INPlUBl.lOl) 

F  QU I  VALENCE  (  [NPl,PRRARAY),(DOPA,GAMMA) 

WRlTE(61,lOQ IkPANGE.NFREQ. Xl.XlSQ,X2.X2SQ,XlPRSQ.X2PRSQ»XlPfi, 
l  AREA1.AREA2. XMULT 

100  FORMAT ( BHOKRAmGE*  13,  3X  7M  NFREQ=  13// 

1  4H  X 1 *611 , 4 ,  3x  6H  X1SQ*  Ell, 4,  3X  4h  X2=  Ell. 4,  3X/6H  X2SQ  =  Ell. 4,  3X 

2  4,  3X  BH  X1PRSQS  Ell. 4,  3x  8H  X2PRSQ=  E11.4.3X/6H  X1PR=  Ell. 4,  3X  7H  AREA 

3  7R  AREA l *  Ell. 4,  7h  AREA  2  s  Ell. 4,  3X  7H  XMULTs  Ell. 4//) 

RETURN 

FND 


Fig.  A-6.  Continued. 
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FUNCTION  aREAfCN(X1.X2,Z1,Z2> 

COMMON  IDiSPLaY(100)»PWRARAY<64,31>. IwTFCN<2,64,28), INPUOl.lOl) 
COMMON  SCATLAW(31> 

C  VARIABLES  IN  COMMON 

COMMON  EL  PRImE,BPRIME,FCL,DOPA,NOISEBX,A,RBSHIFT.ISUB.ISCALE«INP 

1  MIN.  INPMAX,  I  POINT,  I  STEER*  RAD  I  AN.  I DAYNUM, DAYRUNS .  IFLAG.PI  , SMQYRN, 

2  INPMAXl, INPMINI.XMAC, AAHAYN, DEC, ANQ. DOPPLER  A  ,  PL ANET A , FBSH l F T 
DIMENSION  I NP 1(101,101) 

EQUIVALENCE  <lNPl,PtfRARAY),(DOPA, GAMMA) 

ISTEE  =  1 
I F( Zl-1 • 0  )1 , 2 , 1 

2  ISTEE  =2 
X  1 1  =  0 . 0 

1  Dl=l • 0/SQRT ( 1 . 0-Xl*Xl ) 

D2  =  1.0/SQRT(1.0-X2*X2)  S  B2  =  l . 0 /S QRT ( 1 . 0 - Z2 *22 > 

GO  TOO, 4),  ISTEE 

3  B1=1.0/SQRT(1.0-Z1*Z1) 

4  XI22=X2*AC0SF(Z2*D2)-Z2*ASINF(X2*B2)*ASINF(X2*Z2*B2*D2) 

XI12  =  X1*AC0SF(Z2*D1 ) - Z2* A S I NF ( X1*U2 ) * AS  I NF ( XI *Z2 * D1 *B2 ) 

XI2  =  XI22-XU2 

GO  TO ( 5 , 6 ) , ISTEE 

5  XI21=X2*AC0SF(Z1*D2>-Z1*ASINF(X2*B1)+ASINF(X2*Z1*02*B1> 
XI11=X1*AC0SF(Z1*D1)-Z1*ASINF(X1*B1)*ASINF(X1*Z1*B1*D1) 

X  1 1*X 1 21 *X 1 11 

6  AREAFCN  =  XI2-XU 
RETURN 

END 


Fig.  A-6.  Continued. 
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SUBROUTINE  MApXFORM 
C  ARRAYS  IN  COMMON 

COMMON  IDiSPLaY(100)»PWRARAY<64,31).  I WTFCN ( 2 , 64 » 28  )  , INP<10l»l0l> 
COMMON  SCaTLAW ( 31 ) 

C  VARIABLES  IN  COMMON 

COMMON  EL  PRIME, BPRIME»FCL»DOPA,NOISEBX»A,RBSHlFT» I  sub,  I  SCALE.  I  np 

1  MIN. INPMAX, I  POINT, I  STEER. RAD  I  AN. I DAYNUM. DA YRUNS , I  FLAG. PI .SMDYRN, 

2  INPMAX1, INPMINI.XMAC, AAHAYN, DEC. ANQ, DOPPLER  A , PL ANET A » FBSH I F T 
DIMENSION  INPKlOl.lOl) 

EQUIVALENCE  ( I NP1 , PWRARA Y ) , ( DOP A, GAMMA ) , (RB.SlNlNC) 

REAL  ku.kl.kupluskl 
REAL  N0ISEBX1.N0ISEBX2 
C=2.997929E*08  $  DT=5.0E-04 

TWO  A  OV  C  T  = ( 2 . 0* A ) / ( DT*C ) 

DEL  XD  =1.0/<0.992*FCL> 

GAMMA=GAMMA*RaDI AN 
EL  PR  I  ME  =  EL  PR  I ME*RAD I  AN 
B  PRIME=  B  PRIME  *  RADIAN 
COSBPRM=COS(B  PRIME) 

S  I  NBPRM  =  S  I  N ( B  PRIME) 

COSGAM=COS(U*MMA) 

SINGAM  =  SIN<GAmMA  ) 

COSLPRM=COS(EL  PRIMt) 

S I NLPRM=S I N ( El  PRIME) 

DO  91  J=l,101 
DO  91  K  =  1 , 101 
91  I NP ( K, J  )  *0 

KWTMIN  =  37777777B  $  KWTMAXs-KWTM  I  N 

DO  270  K*1 »  28 
DO  270  J=l,64 
DO  270  1=1.2 

I  WTFCN  U.  J.K)=MaX0(1.  I  WTFCN  (  I  .  J .  K  )  ) 

270  CONTINUE 
CONSTANTS 

CONSTANT  =  NOISEBX  +  RBSHIFT*0 .5 
XN0ISEBX*N0ISEBX*1.5  *RBSH  I F  T 
KDEX=NOISeBX*1 
NOISEBXl=CONSTANT+0.5 
NOISEBX2=CONSTANT-0.5 
DXD=0.5*DEL  XD 

DO  509  KRANGE=KDEX,31 


DEPTHr(KRANGE-CONSTANT)/  TWO  A  OV  C  T 
DEPTHl=(KRANGE-NOlSEBXl)/  TWO  A  OV  C  T 
DEPTH2=(KRANGe-N0ISEBX2)/  TWO  A  OV  C  T 
Z1  =  AMIN1(1 . 0. (1 . 0-DEPTW1)  ) 

Z2  =  AMAX1(0 . 0. (1. 0-DEPTH2)  ) 
X1PRSQ«1.0-Z1*Z1 
X2PRS0=1 . 0-Z2*Z2 
COS  I NC  =  1 .0-DEpTH 

SININC  =  SORT( AMAX1( 0 . 0.1 ,0-COSINC*COSINC)  ) 
M=(SININC*0 .135*COSINC)/0 .135 
SL=H*H*H/COSlNC 

NPREOM I N=1  S  NFREQM AX*64 

SL  =  SL  *EXP  (  1.0/COSINC) 
RBSQ=DEPTH*(2.0'DEPTH) 

RB=SORT(R0SQ) 

RB  =  S  I  N I NC 

TERM=64 . 0- ( RB/DEL  XD) 


Fig.  A-6.  Continued. 


87 


NFREQMIN=MAX0(1, IFIXC33. O-TERM) ) 


NFREQMAX=MlN0(64,lFIX(33.0*TERM)> 

c  ************ . ********** 

DO  610  NFrEQs1,o4 

IF(PWRARAY(NFrEQ,KRANGE >>61 0,610. 611 

611  IF (NFrEQ-nFREQMIN)612»613»613 

613  IFtNFREQM AX -NFREQ) 612.615,6 15 

612  PWRARAY ( NFRtO, KRANGE  )  =0 . 0 
GO  TO  610 

615  XD=(NFRE0-33)*DEL  XD 
XI s A8S ( XD-DXD ) 

X2  =  A  BS ( XD-DXD ) 

IF (X 1-1. 0)599, 612,612 
599  IF(X2-1. 0  )596, 612,612 

598  IF(X1-X?)614. 616,614 

616  XMULT  =  2.0 
X1*0 .0 

GO  TO  617 

614  XMULT  si . 0 
XX1=X1 

X  X2- X2 

X2= AMAX1 ( XX1» XX2 ) 

X1=AMIN1(XX1.XX2) 

617  X1SQ-X1»X1 
IF(X1S0-X2PRS0)597,612,612 

597  X2S0=X2*X2 

IF(X2PRS0-X2S0)619.616,618 

619  X2=S0RT<X2PRSq) 

GO  TO  621 

618  IF(X1PRSQ-X2SQ)621,620,620 

620  ARE A1 s 0 . 0 
AREA2sAREaFCN(X1,x2,Z1,Z2) 

GO  T0(  462,463  )  ,SSWTCHF(2) 

462  CALL  TEST  D A T A ( K R A NGE , NF REQ , X 1 , XI SQ , X2 , X2S0 , X1PR SQ . X 2PRSQ 

1  XlPR.AREAl. AREA2, XMULT  ) 

463  CONTINUE 
GO  TO  622 

621  I  Ft XlPRSO-XlSQ >623 . 623, 624 

623  AREA1=0.0 
AREA2=AREAFCN(Xl,x2,1.0,Z2) 

GO  T0i464,4o5).SSwTCHF(2) 

464  CALL  TEST  D AT A ( KR ANGE , NFRE Q , XI , X1SO , X2 , X 2SQ , X1PR SO . X2PRSQ 

1  XlPR.AREAl. ARfcA2, XMULT) 

465  CONTJNUt 
GO  TO  622 

624  X1PR=S0RT ( XIPrSQ) 

AREA1:AREaFCN(X1,X1PR»Z1,Z2) 

area2sAREaFCN(X1PR,X2,1.0,Z2) 

GO  T0(466,467) ,SSWTCHF(2) 

466  CALL  TtST  D A T A  (  KR A NGE  ,  NFREQ . XI , X1SO . X2 , X2SQ , XlPR SO » X2PR SQ 

1  XlPR, AREA1 . AREA2, XMULT  ) 

467  CONTINUE 

622  AREA  =  XMULT* ( AREA1  +  AREA2  ) 

PNRARaV<NFREQ,KRANGE)=PWRARAY(NFREQ,KRANGE)*(Sl  /AREA) 

610  CONTINUE 
509  CONTINUE 

lSU8  =  l+< (MAX0 (1, 1  SUE) )/2)*2 
LIMIT=( ( ISUb-l>/2)*l 
SURSO=ISUB*ISuB*ISCALE 


Fig.  A-6.  Continued. 
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ARG  =  50 . 0*RADIaN 
ERG  =  60.0*RADUN 
XRG=40 • 0*RADI AN 
DELYV=DELXV=RaDIAN 
DDXV=DDYV=DELYV/ISU8 
666  INPMIN-1000000  S  INPMAXa-1 

no  53  Mai, 101 
ENG  =  -XRG*(M-1)*DE|_YV 
DO  45  L  =  1 , 10 1 
BANG=-ERG+(L-1)*DELXV 
SUBPOWERa  o . 0 
SUBLEVEL=0 . 0 
DO  46  J  =  l, I  SUB 
YVANGrENG*  ( J - L I M I T  )  *DD Y V 
YV=SIN< YVANG) 

COSYVANGacOS(WANG) 

DO  46  Na 1 , ISU0 

XVANG=BANG*(N-LIMIT)*DDXV 

XVaCOSYvANG*SlN(XVANG) 

ZV=SORT(1.0-XV*XV-YV*YV) 

Zl=ZV.COSLPRM*XV*S INLPRM 
XlrXV*COSlPRM-ZV*SINLPRM 
Y 1  a  Y  V 

Z2=Zl*COSBPRM*Yi*siNBPRM 
Y2= Yl*COSBPRM-Zl*SINBPRM 
X  2  a  x  1 

XD=X2*C0SGAM*y2*S  INGAM 
YD=Y2*C0SGAM-X2*S INGAM 
ZD  =  Z2 

NFREQ=33.  *<XD/DEL  XD)*FBSMIFT 
GO  TO(630,631),SSWTCHF(3> 

630  WRITE(61,632)M,l, J.N.XV, YV.ZV.XD. YD.ZD.NFREQ 

632  FORMAT (  32H 0 M , L , J . N , X V , Y V , Z V . XD . YD » ZD , NFREO  4  I  5 . 6E14 . 6/ I  5 > 

631  CONTINUE 

IF  (NFRE(J)48»4g,47 
47  I  F  (  NFREQ-64  )  300.300.48 
300  KRANGE=(1.0-ZD)*TWO  A  OV  C  T  *XNO I SE8X 
GO  T0(629, 628) »  SS WTCHF (3) 

629  WRITE (61, 627  )  KRANGE . PWRA RAY ( NFREQ » KR A NGE > 

627  FORMAT ( 8H  KRAnGE*  15,  9H  PwRARA Y:  E14.6) 

628  CONTINUE 
IF(KRANGE-3D302,302.48 

302  KU=FLOAT( IWTFCN(1, NFREO, KRANGE-3)) 

KL=FLOAT(IWTFcN(2, NFREO, KRANGE-3)) 

KUPLUSKL=KU*KL 

IF(YD)303,304,304 

304  SUBP0WER=SUBP0WER*PWRARAY(NFRE0,KRANQE)*(KU/KUPLUSKL) 

SUBLEVEL  a  SUBLEVEL *PWRARAY( NFREQ i KRANGE) 

GO  TO  46 

303  SUBPOWER=SUBPoWER*PWRARAY<NFREQ,KRANGE)*(KL/KUPLUSKL> 

SUBLEVEL =SUBLEVEL*PWRARAY( NFREO. KRANGE) 

46  CONTINUE 

GO  TO(305,3Q6)» ISTEER 

306  SOVERSaSUBLEVEL/SuBSO 
GO  TO  307 

305  SOVERSaSUBPOWER/SUBSQ 

307  IF(SOVERS-XMAc)49,49,680 

49  CONTINUE 

INP(M,L)=SOVErS*0  .5 
INPMAXaMAxO (  INPMAX.  INP(M,  L>> 

I  NPM  I  N  =  M I  NO ( INPMIN. INP(M.L) ) 

GO  TO  45 


Fig.  A-6.  Continued. 
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4  8 


IMP(M,L)=0 


45  CONTINUE 
53  CONTINUE 

ZD  I FF  = I NPMAX- I NPM I N 
ARRAYMAX=JMAX=37777776B 
DO  682  L*l>10l 
DO  682  MS1 » 1 0 1 

68?  INP(M,l>  =  (  ( INP(M.L>-INPMIN)/ZDIFF)*ARRAYMAX 
GO  TO ( 681 , 308 ) , ISTEER 

308  DP=DOPA/RADI an  S  ElP=bL  PRIME/RADIAN  $  BP=0  PRlME/RADI AN 

WRITE (27)  I  DAYnUM. DP. ELP.BP.FCL.NO I SEBX.R0SH I  FT, INPMIN. I NPMAX, 
1  DAYRUNS 
WR  I  TE ( 27  )  I NP 
WR  I  T E ( 2 7  )  I  FLAG 
ISTEER=1 

GO  T0(666,683)  .SSWTCHF (4) 

683  CALL  TESTFOTO(l) 

GO  TO  666 
681  RETURN 

680  SUBSO=SUBSO*10 . 0 
GO  To  666 
END 


Fig.  A-6.  Continued. 
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PROGRAM  HAYFORO 

0  I MFNS  I  ON  RA<50  >.OECS<50  >.NCAY  <50  >  .c  I-ST<50  >»  IDATA<2050  > .  1M102A  > 
1 , IW <1024) .RE <5121 ,RI <5121 »HF<5l2> ,WP< 512) »PH<200  >, INFO <200  ) 
2,CRQT<5l2) .SR0T<512). ARE <51 2 )»  ARI <512) . AHP<5l2  >,  AWP<512  ) 

FQU I V.ALENCF  <IH,IOATA(3)>.<IW,I[:aTa<1027)),<FM!N.ICATA<4)) 

P 1*3.1415026536 
TO  2  I=l»50 

RE  AO  1«  MRS.  AmIN.SECS.  OEGS.CM  INS.  CSEC.OIST<  I  >  .N'OAY<  I  > 

1  FORMAT  <  7  r  1  0  .  7  »  1  3  ) 

RA(I)=HRS*PI/12.0  *AMIN*PI/720.0+SECS»Pl/43200.0 

2  CECS<I>=DFGS*Pl/10O.O*SIGN<(OMINS*PI/lOeOO.O*CSEC*PI/640COO.O 
1  )  .PEGS  > 

3  CONTINUE 
I  P  =  U 

ICnUNT= ICMFCK=U 
PO  4  1=1.512 

A  H  P  (  I  ) = ARP  <  I  )  =0  .  0 
A  R  F  <  I  )  =  A  R  I  <  I  >  =  0.0 

4  RF< I >*Rl < I )=Hp< I )=WP< I  )s0.0 

5  RtFFFR  in  (30.1)  ( IDATA(I), IO*TA(2050)  ) 

ID  r.o  TO  <  10 , 20 . 30 , 40  )  ,  UM  TSTF  <  30  ) 

20  l 0  =  LENG T HF  <  3 0 ) 

I  P  =  0 

IF<LU.EQ.?050)  50,60 

30  PRINT  3 1 

31  FORMAT  < 12H  EOF  UNIT  30) 

GO  TO  60 

40  PRINT  41.jp 

41  FCRMAT<1X.12MPARITY  ERROR. I l 0  ) 

IP  = IP+1 

I F ( IP-5142, 5. 5 

42  PACKSPACF  30 
no  TO  5 

50  ro  200  1=1.512 

II  =  1*1-1  5  1 4  *  I ♦ I 

Nl *  I H ( I1)Sh2=IH( 1 2 ) SRI  *  I W ( 11 ) $ W2  = 1 W  <  12  ) 

HP (  I ) =HP ( I  )  +H1 *H1  +H2*H2 

WP(I>=WP(l)+Wi*Hl+W2*k2 

A=lW<il)*rR0r(I)-iW<l2)*SR0T(I) 

P  =  I  W  < I1)*SR0T ( I )*[W< I2)*CR0T(  I  ) 

RE< I >  *RE  <  I )*IH<  I1  ) *  A  *  I H ( I2)*B 
200  RI < I )*Rl < I >  ♦  I M  < I1)*B*IH< 1 2 )  *  A 
ICOUNT=ICOUNT*1 
IF< IOOUNT-15)  5,210,210 
210  I  COUNT  =  0 
A  A  =  BB  =  0  ,  0 
CO  220  1=92,110 

AA  =  AA*RE  <  I  > 

220  PB  =  BR*R  I  (  I  ) 

PHASE *<180. 0/pI)*(ATAN<BB/AA)*<0. 5-SI GN<0.5«AA))*SIGMPI» BP >) 
PRINT  221,  PHASE 

221  FORMAT(20X»9H  PHASE  =  .F7.2) 

CO  222  1=1.512 

A H P (  I )  =  AHP (  I  )  +HP  (  I  ) 

AWP (  I )  =  AWP (  I  )  +WP (  I  ) 

HP(  I  ) =WP (  I  )  =  0  ,  0 

AL  =  SQRT  <  AA*AA  +  BB*BB  ) 

AC  =  RE(I)*AA/al*RMI)*BB/AL 
AO=RI(I)*AA/AL-RF(D*BB/AL 
RE  <  I  >=RI <  I  )  =  0  .  0 


Fig.  A-7.  CW  phase  calibration  program. 
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A  RF ( I )  =  ARE  < I >+AC 
22?  ARI  (  I  ) 1 AR I ( I )+AD 
I  CHECK  =  1 
RO  TO  5 

100  IF(LD.E0.23)  300.110 

110  PRINT  HI 

111  FORMAT ( 28H  NO  MORE  DATA  LOAC  NEW  TAPE  > 

PAUSE 

RO  TO  3 

300  NDA YT= IDATA ( 1 ) *100+ IDATA ( 2) 

TIME=IDATa<3>+FMIN/60.0 
MOIS=TDATa(1 ) 

IDAY=IDATa(2> 

T  R  =  0 

DO  304  1=1.50 

I F ( NDAYT-NDAY ( I ) )  304,305.304 
305  T  B  =  I 
304  CONTINUE 

PRINT  3  0  6.NDAYT,NDAY< IB)»RA< I B  > »  DECS  < IB) 

7  0  6  FORMAT  <2I20»2p20.5) 

IC=I0+1SID=IB-1 

CRA=RA(I0)*(RA(IC)-RA(ID))*TIMF/48.O 
CDFC=DECS( I B ) + ( DECS ( IC)-DEcS(ID))*TIME/48.0 
ANR=CRA-179.7*PI/180 .0 
NDAYN=(M0IS-8)*31  +  I C  A  Y  *212 

ST  I  ME »  < (TIME/24,0 )♦ <NDA YN-212) )*2.0*P 1*1 . 002737909 
1  +(15.0/24.0+45,0/1440.0  +45 . 868 / < 24 . 0 *36 0 0 . 0 > ) +2 . O+P I 
PN=3245O.OS0A=44.O+PI/18O.O+19.O*P 1/(180. 0+6 0.0) 

BHA=3l.0+PI/l80.0  ♦  22 . 0+P 1/ ( 180 . 0+60 . 0 ) 

PLN=RN+<CnS<CDEC>+SIN<BA)+CCS(BA)+CCS<8HA-STIMF+CRA)+SlN<CrEC>) 

BLW=BN+<COS(Ba>+SIN(BHA-STIME+CRA>) 

BLA=ATAN(PLW/BLN) 

PLI=S0RT(RLN+BLN+BLW+BLW) 

FLN=94.7+(C0S(23.3+PI/180.0)+COS(CDEC)“SIN(CDEr)+SIN(23.3*Pl/lp0.0 

l)+SIN(ANO))+1.068+(RA(IC)-RA(ID))+180.0+30.0/Pl-0.72+CCS(CRA-ST 

2IMF)/DIST( IB) 

FLW=94.7+SIN<23.3+PI/180.0)+CCS(ANG)-1.06e+<DErS(IC)-DECS(lD)) 

1  +  180 • 0  +  30 . 0/P  I 
DOPA  =  ATaN<FLW/F|_N) 

FL=SQRT(FlN*FiN+FlW+FLW) 

ROT=2.0+Pl+SlN<D0PA-BLA)+BLL+6.055/(FL*DIST(lB)+149600.0) 
THFTA=360.0+CoS(DOPA-BLA)*BLL+6. 055/ (CIST(lB)* 149600.0) 

PRINT  307,NDAYN,CRA,CCEC. A NG. STI ME, BLA.BLL.DOP A, FL.RCT.TIF eta 

307  FORMATdX,  I3»1(JF12.4) 

ENCODE<800, 308.  INFO)  NDAYN, TIME. FL. THETA 

308  F0RMAT<4HDAY  ,  I7.1X.5HTIME  . F5 . 2 , IX . 2HFL , F5 . 2 , 1 X , 6H THE T A  ,F7.?> 

DO  309  1=1 .512 

CROT  (  I  )=COS<  (  I-IOD+RCT) 

SROT  (  I  )  =  S  I  N  <  (  I  ■  1 0 1 )  +RC  T  ) 

7  0  9  CONTINUE 
RO  TO  3 

60  IF(ICHECK)  100.100.90 

90  TCHECK= ICOUNT=0 
HANO=WESN=0 . 0 
DO  91  1=257.512 
HANO  =  HANC  +  AHP (  I  ) 

91  WESN=WESN+AWP  (  I  ) 

DO  92  1=1.256 

AHP< I )=AHP( I >+256.0/HAN0  *  1.0 
A WP  (  I )  =  A Wp ( I )+256.0/WESN  -  1.0 
ARF(  I  )  =  ARp( I  )+256.0/SQRT(HANO  +  WESN  ) 


Fig.  A-7.  Continued. 
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ARI  < 1 )  =  AR  f < I >*256.0/SCRT<MANO*WESN  > 

ARI  < I >  =  ARl  < I >/SORT<ABS(AHP< I >*AUP< I > ) ) 

arf<j>=arp<i>/sqrt<abs<ahp<i>*awp(j))) 

52  CONTINUE 

C A I  L  LlMITS<5o. 0.150. 0,-180. 0,180.0) 
no  93  I *5l » 15l 
A  =  I  -  1 

PH (  I  )=ATAN(ARl  (  I  ) / ARE (  I  )  ) 

PH< I > =PH( I >♦ <  0 • 5-SI GN<0 .5. ARE < I >  >  >*S IGN< P 1 , AR I ( 1 >  > 

C  A I  L  POINTS  ( A . PHI  I )*180  .  0/P  I , 1R* . 1  ) 

53  CAI  L  POINTS  <  A  <  1  8  0  .  0  *  S  QR  T  <  A  R  E  <  I  )  *  A  RE  <  I  >  ♦  A  R  [  <  I  )  *  A  R  I  <  I  >  )  ,  1 R  C  .  1  > 

C  A I  L  LABELS ( AhFRPQ. 1 . 3HPHA , 1  ) 

CAI  L  GRIDS(50. 0,10. 0.-180. 0.60.0) 

CAl  L  GRAPHS ( I NfO. 15. 200 . 1 ) 

CAI.L  LIMITS(5o. 0.150. 0,0. 0.180.0) 

CO  94  I *51 . 15i 
A  =  I  -  1 

CALL  PO  INTS ( A , AHP( I  )  . 1RH. 1  ) 

Y  =  10  .  0*AWP(  I  ) 

54  CAI  L  POINTSIA, Y.lRW.l) 

CAI  L  LABELS<4HFREQ»1.3HPOW,i) 

CAI  L  GRICS(  50. 0,  10  .  0.  0. 0.20.0  ) 

CAI  L  GRAPHS! INFO, 15, 200,1) 

HANO  =  HANC/1000  0  .  0 
WESN  =  WESN/1000  0  .  0 

WRITE (2 0.95)  nDAYn.CRA.CDEC.ANG.TIME.ELA.fi  L.DPPA.FL.RCT.TPfcTA. 
1 HANO. WESN 

55  FORMAT  < I  3 , 1 2 F i 0 . 4  ) 

WR!TE<20.96)<AHP<I>,AWP<I),ARE<I>.ARI<I).  [  =  1,200) 

56  FORMAT ( 4E1 8 . 11  ) 

FND  FILE  2" 

BACKSPACE  20 
IF(LD.E0.23)30U,98 

58  CONTINUE 

CAI  L  UNLOAD  ( 30  ) 

PAUSE 

GO  TO  <97,3)  SSWTCHF(l) 

57  FND 


Fig.  A -7.  Continued. 
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PROGRAM  HFPI.OT 
CHARACTER  MAP 

COMMON  ANH(200  ),ANM200  ),RR(  200  ),RI  (200  )#  INTO  (  200  ) 

1, ACSf 200) 

2, ACH(200),ACW(200),AR(200),AI(?G0),PHC<200) 

3, MA(101,101)»MAP<121»121> 

4 » P ( 10  01 ) 

Pl=3.l4l5g26536 
?5l  CONTINUE 
ice*o 

DO  7  L=1'101 
TO  7  N= 1 ' 1 0 1 

7  MA(L'N)=0 

DO  8  K*1 » 121 
CO  8  M*1 » 121 

8  mAP(K,M)=0 

1  C  =  0  .  0 

ndaya=acra=acdec*adopa=afl=atmfta=o . 0 
r-  DO  9  1*1.200 
ACM (  n»o.o 
ACWU>  =  0,0 
aru>=o.o 

9  A  I ( I ) =0  .  0 

10  READ (30,20)  NDAYN, CRA , CDEC, ANG , T IME , BLA ,BLl  » DOPA , FL » ROT, T  FTA 
20  FORMAT! J3.10F10.4) 

R£AD<30.30)  UnH<  I  ),  ANM  I  >»RM! )  <RI  <  1  ) »  1*1.200  > 

3D  FORMATMEie.il) 

ENCODE  <  800  ,68 ,  INFO)  NDAYN,  T  IME  .  FL  .Tt-ET  A 

60  F0RMAK4HDAY  ,  I7.1X.5HTIME  .  F  5  .  ?  » 1 X  .  2  NFL .  F5 . 2 , 1 X  .  6H  THE  T  A  ,F7.?> 
GO  TO  (62,61)  SSWTCHF(2) 

61  CALL  L5MJTS!5o.O»150.0, -180,0,180.0) 

DO  70  1*1,200 

A*  1  "1 

- *ATAn<R!(^)/RR(1)^- 

PH  =  FH  ♦(  0 . 5 -S  I GN  C  0 .5,  RR  (  I  )  )  )«SIGN  (F  I  ,  RI  (  I  )  ) 

PN  =  PH  *180. 0/PI 

0  A  t  L  PG  I N  T  S  !  A  ,  A  tv  H  {  I  ) -180 , 0,1  RH  ,  1 ) 

c  A I  L  PO  I  NTS ( A , ANW (  I)*10. 0-160. 0.1fiU,l) 

CALL  POINTS  (A,. PH  »1R*»1) 

70  CALL  POINTS-! A . 1 80 . 0*SQRT ! RR  ! I )**2*RI ( I )**?) ,1RC»1 ) 

GAIL  LABELS  (4HFREQ.1.3HPHA,1) 

TAIL  GRIDS  (50.0,10.0,-180.0,60.0) 

CALL  GRAPHS-  ( INFO. 15,200, 1 ) 

62  CONTINUE 
PAUSE 

— Gtf  TO!  1*k714-=SSWTcHF (54 

71  DO  80  1=2,200 
K=  I  -1 

-ACHjK)*ACH<K)*ANH( I)— 

ACW(K)=ACW<K)+ANW( I ) 

AR(K)  =  AR(K)*RR( I )*SQRT ( ABS ( ANN ( I ) *  ANH ( 1 ) ) ) 

80  A  I (K)  =  A I (K>*RI ( I )*SQRT( ABS ( ANW ( I >*ANH( I ) ) ) 

C  =  C  +  1  .  o 

nda ya=nda YA*NDA yn 

_ ACRAv=AC^LA*C.flA— — — 

ACDEC*ACCEC+CdEC 

ADOPA=ADCPA+DoPA 

AFL«AFL*FL 

ATHETA=ATHFTA*TMFTA 

GO  TO  (90,10)  S SWTCHF ( 3 ) 


Fig.  A-8.  CW  averaging  and  transform  program. 
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90  CALL  LIMITS  < 5 0  .  0 , 15 0  .  0 . 0 . 0 » 12 0  .  0 > 


DO  100  1*50.150 
P*  I 

CALL  POINTS  <B»ACH< I )/C,lfiH,l) 

100  CAI  L  POI NTS(B, 10 . 0*ACM I )/C«lRW,l > 

C A l  L  LABELS  (4HFRE0.1.3HPCW.1) 

CALL  GRIDS  (50,(1,18 .0.0 < 0.20 . 0 > 
thet=atheta/c 

FNCODF! 800,101, JNFOJNCAYN.TIME.C.THET 

101  FORM AT<7HAVER AGE, 10X.3HDAY, I  7, 4X , 4ht  IPE.FIO . 2. 4X1HC ,F7 . 1,  X. 
15HTHET  A  »  F7 . 1 > 

CALL  GRAPHS! Info. 15, 200,1) 

CALL  LIMITS  <50.0,150,0,-180.0,180.0) 

DO  110  1=50,150 
P»I 

— ARC  H-AR<  l  >/<  ABSOACBI J  )-*SQRT  <  ABS  <  ACK  <  100  >/ACH<  100  )  )  )  )  > 

A  I  <  I  >  =  AI < I )/( aBS< ACH< I )*S0RT<ABS<  ACM100  >/ACH<lOO  ) ) )  ) ) 

PHC  (  I  >  =  A  T  AN ( A  I ( I  )  / A R  (  I  >  ) 

PHC ( I )*PHGti >*<0 , 5-SIGN <0.5, AR( I ) ) >*SIGN  <PL»AI <  J ) ) 

PHC ( I )=PHC< I )*180 . 0/PI 
CAt  L  POINTSiB.PPCU  ),1R*»1> 

AD*SQRT4AfM=f  l*Aft<  f  >*  A  I  <  I  >*  A  I  <  I  >  4 

110  CAI  L  PC  I  NTS < B, 360 . 0*AC-180 . 0  ,  IRC,  1 ) 

CAI  L  lABElS(4hFREQ,1,3HPHA,1) 

CALL  <3Rf  DS<50  ,OT10 .8.-180 ,0,60.0  ) 

PRINT  116  1 

116  FORMAT <1P1.12X,4HhPOW,6X,4hNPCW,6X,4HCORR.5X,5hPNASE» 

-14X,4HCHAN,4X.4HpPOW.6X,4HWPCW,6Xr4HCORR,5X,5HPHASE,4X,4PC  AN, 
DO  112  1=50,148,2 
AB=ACH< I )/C 
AC»ACw( I )/C 

AD  =  SORT ( Ar( I )*AR( I )  +  AI < I )*AI  <  !  ) ) 

AE=PHC< I ) 

J=  T+l 

PB  =  ACH<  j)/c 
PC  =  ACW <  J  > /C 

PD=SORT ( ARC J)*AR( j)*A I ( J)*AI ( v) > 

PE=PHC< J) 

PRINT  111,  AB,  AC,  AD.AE,  I  »  BB ,  BC ,  6D ,  PE 

111  FORMAT! 1  OX, 4F10 .3, I6.4F10 .3,  16 ) 

112  CONTINUE 

CAI  L  GRAPHS ( I  NFC , 15, 200 , 1 ) 

nday=neaya/c 

DEC*ACCEC*180 ,0/(PI*C) 

ANG=ACRA*180 ,0/<PI*C) 

THPTAs ATHFTA/C 

fl=afl/c 

DOPA=ACOPA*180.0/(PI*C) 

ACD  =  ACH  (100  1 
ALPHA  =  0 .11 
DO  520  1*1,1001 
- P  =  (4~-l^y 

520  p(  t ) rExp<-l .0/SQRT <1.0-R*R> >/<R+ ALPHA* SORT <1.0-R*R> )**3 
PXX  =  0  .  0 

DO  5 ?2  K  =  0»1000 

TR  =  SORT(FLOAT(  K*K )  )  +1 . 5 

522  PXX=PXX*P( IR) 

DO  600  1=50,150 
PX  =  cx  =  0  .  0 

I  X  =  A  RS  < I_99.00>*1000.0/FL+0.5 


Fig.  A-8.  Continued. 
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!V=SORT(1000000.0-IX*IX)*0.5 
ro  550  K - 0  »  IX 

i r= sort (Float ( ix*ix*k*k>  >+i  .5 

PX=pX*p ( Ir ) 

550  CX=CX*COS(PI*THETA*K/180000.0)*F(IR) 

ACF=ACH(  I  ) 

ACM  I )=ACH< I ) -PX*ACD/FXX 
CR=CX/PX 

PRINT  560,  ACH<1> ,PX.CR 
560  rORMAT<10X*3F3.0,A) 

ACS(  I  )  =  ACM  I  >/<  ACP  +  ACE/100 . 0  ) 

AR ( I >  =  AR(  I  )  -CR*0 , 9 

AR( I )»AR( I )*AcP/< ACP+ACD/100 .0 ) 

600  A  I ( I )  =  AI ( I )*AcP/( ACP+ACD/10  0 .0 ) 

r  smoothing 

rC  601  1=52.146 
Tl=I-2*I2=I-i$I3=I*l$I4=I+2 

RR(  I  )=AR(  II  )  *2  •  0  *  AR ( 12  )  *2 . 0  *  A  fi ( I  )  *2  . 0*AR( I3)*AR( 14) 
RI<I>=AI(ll)*2-0*AlU2)+2.0*A!(I)+2.0*AI<I3)+AT(I4) 

601  A  N  H  (  I  >  =  ACS(11  >*2 . 0*ACS( I  2 )  +  2 . 0* ACS( I )*2 . 0*ACS< I  3 )*ACS( M 
CC  602  1=52.148 

A R ( I  )=RR(  1)64.0 
A  I  ( I )=RI <  I )/4. 0 

602  A  C5  ( I >  =  A  N  H  < I )/4,0 
AR(50)=AR(51)=AR(149)=AR(150)=0.0 
AI ( 5  0 )  =  A  I ( 51 >  =  A  I ( 1 49  )  =  A  I  (15C)=0.0 
ACS(50 ) SACS (51  Abacs (149)  =  ACS (150)  =  0.0 

c  plank  out  the  SLR  RADAR  PCINT 

AC':(95)=ACS(96)=ACS(97)=ACS(103)=ACS(104)=ACS(105)=0.0 
— 4*8--«0  3  1=08.102 

603  A  R ( I  )  =  A  I ( T )  =  AcS( I)  =  0 . 0 

TAIL  LIMITS  <50.0,150.0,-180.0.180.0) 

CO  620  1=50,150 
PM 

PMC (  I  )  =  ATA N ( A  I  ( I  )/AR(  I  )  ) 

-PH044)  =  PH€(-LJ+(0.5-SIGN(0 .5.  Afi(  I  )  )  )  * S  I  GN  (  P  I  .  A  I  (  I  )  ) 
p  h  r ( i )=pHn< I )*ieo:o/Pl 
TAI  L  PC  INTS (B, PRC ( I ) . IR* , 1 ) 

CAOt  PC  I  NTS (R, *CS( I ) *360 ,0-180 .0,  IBM  .1 ) 

620  CAI  L  POINTSCB, (360.0  ) *SCRT ( A R ( I ) * « 2* A  I ( I ) **2 ) - 1 8 0  .  0 , 1 RC  1) 

TAIL  LA0ELS(4hFRFQ,1,3hPHA,1) 

CAi L  6RJ OS (50 . 0, 10 . 0.-180.0,60 . 0 ) 

C Al  L  GRAPHS  (  0  »  0 . 0 . 1 ) 

C  FILL  MA  ARRAY 

CO  130  1=50,150 
CfisAR (  I  >*10000  .  0 
C  I  = A  I (I >*10000.0 
C A  =  ACS ( I ) *10000 ,0 
TC  130  N= 1 » 1 0 1 
L  =  T -  49 
KF=N-51 

130  MA(L,N)=MA(L»N)*CR*COS(THETA*FI*XF/SOOO.O)  *CI*SIN (THETA*  I* 
1XF/9000.0  )  -CA*  ( 1- ICC  ) /2 . 0 
!CC  =  1 

r Al  L  LIMITS (1.0, 101. 0,51. 0,111.0) 

ro  i4o  n= 51 » i o i 

r.n  i40-  i-=i  ,iOi 
Ai  =  N 
A2  =  L 

MO*MA(L.K)/10oO.O  *0.5 
140  CAI L  P0INTS(A2.A1,MQ,1) 

CAI  L  GRAPHS  (  0  «  0  .  0 , 1 ) 


Fig.  A-8.  Continued. 
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(“Ail  LIMITS  (1.0.101.0,-9.0.51.0) 

TO  15>0  N  =  1 »  51 
TC  150  L=l»10l 
Al  =  N 
A  2  =  L 

pQ  =  MA(I.»N)/1000,0  +0  .5 
150  r  A  |  L  PCIMS(A2,A1,MQ,1) 

CAI  L  GRAPHS  ( 0  .  0 , 0 . 1 > 

PA1.SE 

RC  TO  ( 155, 1  >  SSWTCHf ( 6  > 

155  ANG=A|\|G-179.7 

rDFC  =  CCS (OFC*P  1/180.0) 

SCFC  =  S  I  N  < DFC*P  1/180 . 0 ) 

CANG=COS< AMG*P  I /ISO  .  0  ) 

SANG=S I N( ANG*P  1/180.0) 
rTN=OoS(23.3*PI/180.0) 

?Th  =  SIM23.3»PI/180.0) 

ALATR=(18P.0/pl)«ASINF(-CTH*SCEC-STH*SANG*rDFC) 

ALONR=( 180.0/p  I )*aTAN(<SANG*CCEC*CTF-SDEC*STH)/CANG*CDEC) 

ALONR  =  ALCnR  +  16.0+(NDAY-241)*360. 0/245.0 
PLA=ATAN(CANG*STH/(CTH*CDEC-STH*SCFC*SANG)) 

APa  =  PlA*180  .O/PI 

rrop=cos (nopA+p  i  /  iso .  o ) 

RDrP=SlN(POPA*PI/l80.0) 

NRUN=NnAY*10000+TjME*100 

PRINT  160, NR Un.NOaY.DEC.ANG. THETA, FL-EOPA.aLaTr.ALCNR.AFA 

160  FORMAT ( 21 10 , 8F10  .  2  ) 

PE=PI/180.0«(l6.0+(NnAY-24l)*360. 0/245.0) 
ro  180  L *  1  *  1 0 1 

rc  180  N=1,101 

Y=(L-51)*CPCP/FL+(N-51)*SCOF/50.0 
Z=(L~51)«snOP/FL  -  (N-51 > *CCOF/50 . 0 
IF(1-Y*Y-Z*Z>160,161,161 

161  Y=SQRT (1- Y*Y-Z*Z )  ’ 

ya=x*cang»odec-y*sang+z*scec*cang 

YA=X+(CTH*SANG*CDEC-STH+SCEC)+Y*rANG+CTH+Z+(CTH+SCEC+SANG  STH+oDFC 
1  ) 

7A=X*( -CTN*SDeC-STH*SANG*CDEC)-Y*CAKG*STH+Z*(CTH*CCEC-STH  cDFC*SAN 
IP) 

xe=XA+COS(PE)-TA*SIN(FF) 

'YB  =  XA*SIMPE>  +  Ya*COS(FE> 
k=yB*60 .0*61.5 
p=7A*60.0  +6a . 5 
IF(XB)  170,175,175 
170  MAP(K»M)=lR$ 

GO  TO  180 

175  8AP(K,M)=MAP<K,M>+(MA(L»N)*XB/1000.0) 

180  CONTINUE 

TO  190  J=l,19 

aa=(j-io)*pi/i8.o 

CO  190  1=1,361 

BB=( I-101)*PI/36O.O 
K=COST  AA)*S1N<PBM60t0*61i5 
n  =  S IN( AA ) *60 . 0*61  .5 
IF (MAp(K-M) )  190,181,190 

181  NAP (K,M ) »1R* 

ISO  CONTINUE 

TO  200  J=1 , 19 
AA  =  ( J-10 ) *P  I  /18 . 0 

rc  200  1=1,361 

PB= ( I-101)*PI/36O.O 


Fig.  A-8.  Continued. 
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K=rOS(BB)*SIN(AA)*60.0+61.5 
m=SIN(BB)*60.0  *61.5 
IF<MAp(K,M))  200,191,200 
151  MAP(K,M)*1R* 

200  CONTINUE  I 

CAI  L  UMITSd.O  ,61. 0#  61. 0/121.0) 
rc  210  K  =  1 »  61 
ro  210  M=61.121 
PAPK=MAP<KfM) 

AK  =  K 
AM  sM 

f-P  f  M  A  P  K  >2  0  9 , 2 1 C  >'  2  0  9 

209  CAI  L  PCIMS(AKfAM.MAPK.l) 

210  CONTINUE 

CALL  GRAPHS  ( Q»0 / 0> 1 ) 

CALL  LlMlTSifOl  .0*121 .0/61 .0/121 .0  ) 
CO  220  K*61/I2l 
TO  220  M=6l » 1 2 1 
AK  =  K 
AM  =  M 

MAPK  =  M Ap ( «  »  M ) 
tF{MApK)2l9. 220,219 
215  CALL  POINTS(AK/AH,MAPK,l) 

220  CONTINUE 

CAI  L  GRAPHS  (  0  /  0  /  0  >  1 ) 

CAI  L  LIMITS  (1.0  ,61.0/1.0/61.0) 

DO  230  K=l,61 
CO  230  M*1 »  61 
AK  =  K 
AM  =  M 

MAPK  =  MAP ( K  >  M  ) 

TF(MApK)229,230,229 
229-OAfL  POINTS!  A«#  AM, MAPK/l) 

230  CONTINUE 

CAI L  GRAPHS  <0» 0/0/1) 

MAIL  t  IttI  TS<  6i  i  0x"l21  tO,  1  tO/-<1  •  G  ) 

CO  2 A 0  K  =  61fl21 
CO  240  MS1 , 61 
AH«K  _  — 

AMSM 

HAPK*MAP(K,M> 

IF(HAPK)239,240,239 

239  CAIL  POINTS(A«.AM,MaPK,1) 

240  CONTINUE 

CALL  GRAPHS  10,0,0,1) 

PAUSE 

SO  10(242,241}  SSWTCHF  < 1 ) 

241  CALL  SPFTTJOi 
QACKSPACE  20 

242  CONTINUE 


pro  nte  20 

BACKSPACE  20 

SO  TO  (251/250)  SSHTCHF(4) 
250  FND 


Fig.  A~8.  Continued. 
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PROGRAM  HFMAP 
CHARACTER  HP 

DIMENSION  MAP (121, 121) ,MP(l2l,121> 

PI*3,l4l5v2653 

IO«0 

IC»0 

DO  1  K*l»l2l 
DO  1  M*l»  l2l 

1  m AP ( K , M ) *  o 

2  RE AD ( 30 )  MP.NDAY 

GO  TO(30»3)»£OfCKF(30) 

3  GO  TO  ( 2 » 1 4 ) » SSWTCHF ( 2 ) 

14  PAUSE 

DO  10  K«l,121 
DO  10  M*l,121 
IF(MP(K,M)  ,tQ,lR*)4,5 

4  MP(K,M)=0 

5  IF(MP(K,M) ,EQ.1RS)6,7 

6  HP ( K »  M ) s  0 

7  I  A*MP ( K, M ) 

J  F  (  I  A  -  3 1 )  9 »  9 , 8 

8  J  A  *  I  A  -  64 

9  MAP(K,M)»MAP(K*M)*IA 
10  CONTINUE 

PRINT  20 , no a V 

20  FORMAKlOX,  JHDAY.4X,  13) 

I C= I C*1 

GO  TO  (30,2),  SSWTCHF ( 3 ) 

30  DO  40  K*l,  121 
DO  40  M*l, 121 
MD*MAP(K,M) 

40  MAP(K,M)imD 

41  DO  190  J*l,l9 
AA=( J-10)*PI/18,0 
DO  190  I  *  1 , 36 1 

69  = ( I *181 )  *P I /360 . 0 
K»COS(AA)*SIN(BB)*60  .0*61.5 
M»SIN(AA)*60.0*61.5 
IF(MAP(K,M))  190,181,190 
181  HAP(K,M)»1R* 

190  CONTINUE 

DO  200  J*1 *  19 

AA«(J-10)*PI/l8.0 

DO  200  1*1, 36i 

69* ( I-181)*PI/360 .0 

K*COS(BB)*SIN(AA) *6  0 .0*61.8 

M*SIN(8B)*60.0  *61.5 

I F ( MAP ( K . M  )  )  200,191,200 

191  MAP(K,M)*lR* 

200  CONTINUE 

CALL  LlMlTSd.O  ,61.0.61.0,121,0) 
DO  210  K *  1 , 61 
DO  210  M«61,121 
MAPK=MAP(«,M) 

AK*  K 
AM*M 

IF(M*PK)209, 210,209 
209  CALL  POINTS! AK, AM, MAPK.l) 


Fig.  A-9.  CW  display  programs. 
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210  CONTINUE 

CALL  GRAPHS  <0#0#0.1) 

CALL  LIMITS (6l. 0.121. 0,61.0 .121,0) 

DO  220  K*61 »  121 

DO  220  M*61  *  121 

AK*K 

AM»M 

MAPK*MAP(K»M) 

IF(MAPK)2l9,220,2l9 

219  CALL  POINTS! Ak. AM, MAP*, 1) 

220  CONTINUE 

CALL  GRAPHS  <0*0, 0»D 

CALL  LIMITS  (1.0  ,61,0,1,0,61.0) 

DO  230  K* 1  *  61 

DO  230  M*l, 61 

AKxK 

AM»M 

MAPKSMAP(K.M) 

IF(MAPK)229,230,229 

229  CALL  POINTS! A«, AM, MAPK,1) 

230  CONTINUE 

CALL  GRAPHS  <0*0»0*1> 

CALL  LIMITS(61. 0,121.0, 1.0, 61,0) 

DO  240  K*  61  *  121 
DO  240  M*1 » 61 
AK*K 
AM*M 

MAPK*MAP<K»M) 

IF(MAPK)239,240,239 

239  CALL  POINTS! A«, AM, MAPK.l) 

240  CONTINUE 

CALL  GRAPHS  <  0  *  0  »  0  *  1 > 

WR I TE ( 20 ) M AP  *  NDA Y 
end  file  20 

DO  245  K«1.12l 
DO  245  M"1 *  121 
245  MAP(K,M)»0 
IO*IO*l 
REWIND  30 

250  READ<30)  MP»ND*Y 

GO  TO(4l,300),EOFCKF(30) 

300  DO  310  K*l,l2l 
DO  310  M*l,l2l 
I F ( MP ( K , M ) .EQ.1R*)304,305 

304  MP ( K , M ) * 0 

305  I F  <  MP ( K , M  )  ,EQ.1RS)306,307 

306  MP ( K  *  M  )  *  0 

307  IF(MP(K,M)-I0)3l0,308,306 

308  IF(MP(K,M)-30)309,309,310 

309  MAP(K,M)«MAP(k#M)*1 

310  CONTINUE 
GO  TO  250 
END 


Fig.  A-9.  Continued. 
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PROGRAM  HFDISP 

DIMENSION  MAP(121»121),X(361)»Y(361) 

REWIND  30 

READ  <  30 )  (<MAP(K»M)»M*1,121)»K=1»121  ).NDAY 
P I  =  3 • 14159265 

1  CALL  PESETDI -6, -10  ) 

CALL  SCALE<-1.,1..-1.,1.,32,32.0.64) 

DO  10  J= 1 » 1 9 

AA=< J-10)*PJ/l8.0 

DO  2  1=1.361 

BB= ( 1-181) *P 1/36 0,0 

X< I)=C0S(AA)*SIN<BB) 

2  Y< I ) =S IN<  A  A) 

DO  3  K  =  1 »  5 

3  CALL  PL0K2. 360, X.Y, 1,1,0) 

10  CONTINUE 

DO  20  J  =  1 , 19 

AA= ( J-10 )*P 1/10 . 0 

DO  12  1=1,361 

BB=  < 1-181) *P 1/360,0 

X(  I)=CCS(BB)*SIN(AA) 

12  Y ( I )=SIN<Bfi) 

DO  13  K  =  1 , 2 0 

13  CALL  PL0T(2, 360, X.Y, 1,1,0) 

20  CONTINUE 

DO  24  K  =  1 , 121 
DO  24  M=1 , 121 
IP (MAP(K.M)-lR*)  24,22,24 
22  MAP(K,M)=0 
24  CONTINUE 
IMAX  =  -1  .E10 
I  M I N  =  -  I M  A  X 
DO  50  1=1,1*641 
I M  I N  =  M I  NO ( I M I N « M A P ( I ,1 ) ) 

50  IMAX  =  MAX0 ( IMAX,MAP( I , 1 ) )  $  W R  I  TE ( 5 9 , 6 0  )  I  M  I  N .  I  M A X 

60  FORMAT (10HIMJN, IMAX=215) 

WRITE(59.40) 

40  PORMAT ( 21hWR  I  TE  I  M  I  N . I M A X , NGR A Y ) 

CALL  GE  TDaT  A ( 50 , A , B , C  ) 

IMIN=AIIMAX=HJNGRAY=C 

call  INTRPLOT(MaP,121,121,32,2,0,2,4,4,0,NGRAY, IMIN, IMAX) 
CALL  ADVPILM(O.I) 

PAUSE 
GO  TO  1 
END 


Fig.  A-9.  Continued. 
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